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ABSTRACT 


This report presents the results of work on the modeling of com- 
munication receivers at both the circuit detail level and at the block 
level. The largest effort was devoted to developing new models at the 
block modeling level. The available effort did not permit full develop- 
ment of all of the block modeling concepts envisioned, but idealized 
blocks were developed for signal sources, a variety of filters, limiters, 
amplifiers, mixers, and demodulators. These blocks were organized into 
an operational computer simulation of communications receiver circuits 
identified as the Frequency And Time Circuit Analysis Technique (FATCAT) . 
The simulation operates in both the time and frequency domains, and per- 
mits output plots or listings of either frequency spectra or time wave- 
forms from any model block. Transfer between domains is handled with 

a fast Fourier transform algorithm. 

A separate block model effort was devoted to developing block models 

for use with the circuic simulation MIMIC. 

Two efforts in modeling at the circuit detail level were also car- 
ried out. One of these demonstrated the feasibility of interfacing the 
time-domain analysis program CIRCUS with a frequency-domain analysis pro- 
gram to provide a more powerful model. The other effort demonstrated 
the use of ECAP to determine input parameters for CIRCUS which are other- 
wise difficult to determine. 

Computer listings of the software developed and examples of use of 
the programs are included. 


pTj * 




-■ 17 ) 


FOREWORD 


This report was prepared at the Engineering Experiment Station at 
the Georgia Institute of Technology for the Astrionics Laboratory of 
Marshall Space Flight Center under Contract NAS8-28148. The work was 
carried out under the direct supervision of Mr. J. R. Walsh, Project 
Director, and under the general supervision of Mr. D. W. Robertson, 
Chief of the Coranunications Division. The report describes the results 
of a one-year effort on the modeling of communications systems. 


PAGE 


* 


TABLE OF CONTENTS 


I . INTRODUCTION 

II. MATHEMATICAL MODELING OF COMMUNICATIONS CIRCUITS ON A 
CIRCUIT DETAIL BASIS 

A. Interfacing CIRCUS and the FFT 

1. Time Domain vs Frequency Domain Analysis . . 

2. Description of Interface and Software . . . 

3. Example of Using the CIRCUS-Frequency Domain 

Interface Program 

B. Use of ECAP for Transfer Function Evaluation . . 

C. Conclusions and Recommendations 

III. MATHEMATICAL MODELING OF CCK2JNICATI0N SYSTEMS 

ON A BLOCK BASIS 

A. General Description 

B. Program Control 

I. Block Input Commands 



a. 

Sources . . 


b. 

Filters . . 


c. 

Demodulators 


d. 

Other Blocks 

2. 

Control Commands 

3. 

Special Commands 


C. Simulation Examples using FATCAT 

1. AM Receiver 

2. FM Receiver 

D. Conclusions 

IV. TIME DOMAIN SIMULATION 

A. Introduction 

B. MIMIC: A Continuous System Simulation Language . . . 

C. TIMSIM: A User-Oriented Communications System Block 

Diagram Simulation Program 

1. Subprogram Derivation 

2. TIMSIM Subprogram Library 

3. Applications 

D. Conclusions 


1 

5 

5 

5 

6 

8 

19 

25 

29 

29 

30 

30 

31 

32 

33 

33 

34 
37 
39 
39 
52 
57 
71 

71 

72 

79 

79 

81 

81 

95 


v 


TABLE OF CONTENTS (Continued) 


V. REFERENCES 


PAGE 

96 


APPENDICES: 
APPENDIX A 
APPENDIX B 
APPENDIX C 
APPENDIX D 


97 

101 

117 

179 





LIST OF FIGURES 


Test Circuit for the CIRCUS Frequency Domain Interface 

Program 

Execution of CIRCUS for Test Circuit 


Execution of Interface Program 


Spectrum of Data Set 3 

Time Waveform of Data Set 3 . . 

Request for Calcomp Plot ..... 

Calcomp Plot of the Frequency Spectrum of the Time 


Waveform Shown in Figure 5 

Circuit Diagram of CCS Telemetry Filter 

Circuit Description in ECAP Format 

Comparison of Response Calculation for CCS Telemetry 


Filter 


Sample of ECAP Output 

Block Diagram of Amplitude Modulation Receiver and 
FATCAT Model of Receiver 

Start of FATCAT Run for AM Receiver, Showing Spectrum 
at the Output of the Signal Generator 

Listing of the Frequency Function at the Output of 
the Signal Generator 

Listing of Frequency Function at the Output of the 
First Bandpass Filter 

Spectrum at the Output of the First Bandpass Filter . 

Spectrum at the Output of the Ideal Multiplier . . . 

Spectrum at the Output of the Second Bandpass Filter 

Spectrum at the Output of the Amplifier 

Spectrum at the Output of the AM Demodulator .... 

Time Waveform at the Output of the AM Demodulator . . 

Listing of Time Function at the Output of the AM 

Demodulator 

Spectrum at the Output of the AM Demodulator .... 

Spectrum at the Output of the Low Pass Filter .... 

Block Diagram for FATCAT Analysis of an FM Receiver . 


vii 


* 


LIST OF FIGURES (Continued) 


PAGE 




26. Start of FATCAT Run for FM Receiver, Showing Spectrum 

at the Output of the Signal Generator 

27. Spectrum at the Output of the First Bandpass Filter 

28. Spectrum at the Output of the Ideal Multiplier 

29. Spectrum at the Output of the Second Bandpass Filter .... 

30. Spectrum at the Output of the Amplifier 

31. Spectr um at the Output of the FM Demodulator 

32. Time Waveform at the Output of the FM Demodulator 

33. Listing of the Frequency Function at the Output of 

the FM Demodulator 

34. Spectrum at the Output of the Low Pass Filter 

35. Time Waveform at the Output of the Low Pass Filter 

36. Listing of the Frequency Function at the Output of 

the Low Pass Filter 

37. Detailed MIMIC Block Diagram for a System Described by 

x + x + x = 0 

38. MIMIC Instructions for Simulation of the System of 

Figure 37 

39. Block Diagram of Feedback Example 

40. MIMIC Simulation of System of Figure 39 

41. Subprogram Example 

42. Single-Tuned Filter Flow Diagram 

43. Hypothetical AM Communication System 

44. TIMSIM Listing, Hypothetical Communication System 

45. TIMSIM Results, Hypothetical AM Communication System . . . . 

46. Plot of TIMSIM Results for AM System 

47. Hypothetical Automatic Gain Control System 

48. TIMSIM Listing, Automatic Gain Control System 

49. TIMSIM Results, AGC System 


58 

59 

60 
61 
62 

63 

64 

65 

66 

67 

68 

74 

74 

76 

76 

78 

80 

84 

85 
88 

90 

91 

92 
94 


| 



viii 


>■* 


I. INTRODUCTION 

This report discusses the results of a one-year effort directed 
toward the modeling of communications receivers. This effort followed 
a program carried out during the previous year under contract NAS8-20054 
in which transmitter modeling was considered. Most of the modeling work 
carried out on the former program was directed toward modeling at the 
circuit detail level, and the simulation program CIRCUS was used as the 
primary analysis tool. Accomplishments under that program have been 

documented in previous reports [1,2]. 

On the current program, modeling of communications receiver circuits 
has been carried out on both the circuit detail level and the block model- 
ing level, with emphasis on the latter. Several existing circuit simula- 
tion programs have been investigated and three of these (CIRCUS, ECAP, 
and MIMIC) were used in the modeling efforts reported herein. In addi- 
tion, using a block modeling approach, work on a new circuit simulation 
program was undertaken and the program was developed to an operational 
state that provides a highly useful circuit analysis tool. However, the 
present form of the program is far from realizing the potential that the 
approach offers. 

Two lines of effort were pursued in circuit detail modeling. One of 
these was an investigation of the feasibility of interfacing CIRCUS with 
some other circuit analysis program which operates in the frequency do- 
main, thereby creating a more powerful analysis tool. The other effort 
was concerned with using ECAP as a basic tool for providing certain in- 
put parameters for CIRCUS which are difficult to determine otherwise. 

The results of both investigations were promising. 

The need for interfacing CIRCUS with some other program arises from 
th fact that CIRCUS operates entirely in the time domain. Processing 
of signals requires step-by-step integration of the time function. Such 
a procedure is necessary when dealing with nonlinear portions of a cir- 
cuit; when linear circuits with known transfer functions are involved, 
the effect of the circuit is much more readily computed in the frequency 
domain. Thus a model with the capability of processing nonlinear circuits 
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in the time domain, and linear portions in the frequency domain would be 
very attractive. Since CIRCUS is a very powerful time domain simulation, 
it would provide a good starting point for creating a two-domain model 
if it could be successfully interfaced with a frequency domain model. 

The investigations conducted were really in the nature of feasibil- 
ity investigations. Methods were developed for extracting the computed 
time functions from CIRCUS, converting them to the frequency domain, and 
displaying either time or frequency functions. This work is discussed 
in Section II and all software developed in this effort is listed in 
Appendices A and B. 

The work with ECAP was directed toward using the a-c analysis cap- 
abilities of ECAP to investigate circuit parameter effects on the trans- 
fer functions of certain circuits, and thereby determine suitable para- 
meter values for entry into CIRCUS. There are certain parameters requir- 
ed in specifying circuits to CIRCUS which are not easily estimated. One 
example is the values of coupling coefficients in coupled circuits. Un- 
less reasonably good values for these parameters can be determined prior 
to entry into CIRCUS, interpreting the program's output may be very near- 
ly impossible if it runs; sometimes the program will not even operate if 
the values are unrealistic. 

The investigations reported herein show that ECAP is a useful tool 
for determining parameter values in many cases. Details of the investi- 
gation are also discussed in Section II. 

The major effort in this program was devoted to circuit models based 
on a block modeling (as opposed to circuit detail modeling) approach. 

The block modeling work can also be divided into two lines of effort. 

The primary effort was devoted to developing models for circuit blocks 
and combining them to create an entirely new simulation program. The 
other effort consisted of investigating the use of the simulation pro- 
gram MIMIC in modeling circuits with feedback loops on a block modeling 
basis . 

Development of the new simulation program was directed toward model- 
ing each subsystem (mixers, filters, amplifiers, etc.) as a single block 
in order to obtain an efficient program for rapid analysis of communica- 
tions receivers. The simulation was planned from the beginning to operate 
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in both the time and frequency domains. A representation of the signal 
is stored in a complex data array, and all processing blocks operate on 
this array. The array contents at any given time may be either a dis- 
crete spectrum representing the frequency function, or data points repre- 
senting samples of a time waveform. The array contents are transformed 
from one domain to another as needed by a fast Fourier transform (FFT) 
algorithm. 

The block model for each subsystem can be modeled in either the time 
domain or the frequency domain. Blocks have been created for signal 
sources, a variety of filters, mixers, limiters, amplifiers, and demodu- 
lators. Other needed software was developed to create an operational 
simulation identified as the Frequency and Time Circuit Analysis Tech- 
nique (FATCAT). Control software was designed for conversational opera- 
tion from a real time computer terminal; it is readily adaptable to 
batch processing, however. 

The program was developed and completely implemented on a Univac- 
1108 computer at Georgia Tech. Required modifications were then made to 
adapt the program to a SIGMA-5 computer at Marshall Space Flight Center. 
The program is discussed in Section III and descriptions and listings of 
all software for the Univac-1108 version are given in Appendix C. List- 
ings of the SIGMA-5 versions of those routines requiring modification 
are given in Appendix D. 

Although the version of FATCAT presented here was brought to an opera 
tional state and has proved to be a highly useful simulation, it is by 
no means the ultimate simulation of its type. Many ideas exist for im- 
proving the program which could not be investigated during the time avail- 
able. Ideas for improvement include changes to upgrade some of the exist- 
ing model blocks, creation of new model blocks, and additional features 
that could be added to the control software which would make the simula- 
tion even easier to use. Some specific examples of improvements that 
should be developed are pointed out in Section HID. 

The other block modeling effort carried out under this program made 
use of the simulation program MIMIC, a dynamic system simulation program. 
MIMIC accepts user inserted "models" in the form of equations, and this 
feature was investigated as a method of modeling circuit blocks. The in- 
vestigation is discussed in Section IV. 
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II. MATHEMATICAL MODELING OF COMMUNICATIONS CIRCUITS 
ON A CIRCUIT DETAIL BASIS 

A. Interfacing CIRCUS and the FFT 

1. Time Domain vs Frequency Domain Analysis 

One of the major efforts related to circuit detail modeling 
was devoted to interfacing CIRCUS with other programs so that frequency 
domain representations of signals could be obtained. CIRCUS itself 
operates entirely in the time domain. Processing the signal requires 
step-by-step integration, often with very small steps. Although the 
process can be time consuming, it is a good method for performing analy- 
ses that must be performed in the time domain, such as analyzing non- 
linear circuits or calculating the transient response of any circuit. 

If only the steady state solution is desired for linear circuits, 
frequency domain analysis can be applied. Since frequency domain analy- 
sis is generally much less time consuming, and also since the steady 
state response is of prime interest in most circuit investigations, it 
would be preferable to analyze linear circuits in the frequency domain. 
However, entire equipments (such as a communications receiver) cannot 
usually be analyzed entirely in the frequency domain since they usually 
contain nonlinearities. The ideal solution would be to have an analysis 
program that operates in both domains and to have the ability to trans- 
fer from one domain to another. 

The investigation reported here was directed toward investigating 
the possibility of combining circuit detail models which operate in the 
two domains, and thus create a single more powerful model with more effi- 
cient operation. The job of actually creating such a model was beyond 
the scope of this investigation, and only a feasibility study has been 
conducted . 

Specifically, the task undertaken was to develop an interface with 
CIRCUS which would permit extraction of computed time waveforms. These 
signals could then be operated on with a fast Fourier transform (FFT) to 
obtain frequency domain representations of the signals. Such an inter- 
face has been constructed, along with some supporting software that 

TX ( r W ^ “v ' rrrp 

5 ' 




permits recovery of the CIRCUS-generated signals, transforming between 
domains, and output of results. Details of the interface and software 
are described below, followed by sample runs. 

2. Description of Interface and Software 

For purposes of this feasibility study, the modifications made 
internally in CIRCUS were kept to a minimum. The immediate objective 
was to gain access to the computed time waveform data points, and to 
write these into a file external to CIRCUS so that the data would be 
stored when the run terminated. Other software was then constructed to 
access the stored data file and process it. 

After examining several of the CIRCUS subroutines, it appeared that 
the desired data could be obtained easily by constructing a subroutine 
to be called by the PLOT statement. As a working procedure, subroutine 
PLOTER was disabled (by commenting out its call in subroutine LINK6A) , 
and a new subroutine, PLTDTA, was constructed to extract the data. 

PLTDTA writes the data out on one of the FORTRAN output units (currently 
set at 19) and then terminates. This form of modification permitted 
testing without the necessity of revising the CIRCUS command structure. 
The PLOT statement could be used to specify the outputs wanted, and out- 
puts from as many nodes as desired could be obtained. Thus the changes 
in CIRCUS were limited to construction of PLTDTA plus minor modifications 
to LINK6A. The modifications made in LINK6A were (1) deleting the calls 
to PLOTER by inserting comments on line, and (2) inserting calls to 
PLTDTA at the appropriate point. Listings of PLTDTA and the modified 

version of LINK6A are given in Appendix A. 

The output generated by PLTDTA consists of two parameters and three 
arrays. In order these are: (1) the parameter NPNT specifying the num- 

ber of points per data set; (2) the parameter NPLOTS specifying the num- 
ber of data sets to be outputed; (3) the array TIME containing NPNT 
entries of the times corresponding to output data points; (4) the array 
PLOT containing NPLOTS sets of output data, each consisting of NPNT 
values of a node voltage, an element current, or whatever variable was 
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called for in the PLOT statement;* and (5) the array TITLE containing 
the title of the run (not used by the current version of the interface 
program) . 

In using the program, the modified version of CIRCUS is run with a 
PLOT statement included that has a list of the variables to be written 
out. When the computations are finished the CIRCUS run is terminated. 

The interface program can then be activated to call the data file and 
operate on it; the call is non-destructive and the same data file can 
be processed many times by the interface program. 

There is one specific requirement that the user must be aware of 
and provide for. The FFT is designed to work only with data sets of N 
points where N is a power of two (N = 2 IGAM , where IGAM is an integer). 
Since the theory underlying the FFT requires that it be applied only 
to periodic signals, the user must pre-determine the period of the out- 
put waveform to be generated by CIRCUS. The period must then be divided 
by the chosen value of N (with N = 2 IGAM ) to determine the time step 
at which outputs will be generated. This procedure insures that exactly 
N consecutive outputs will represent exactly one period of the output 
waveform. 

While the time step size must be carefully determined, the overall 
time interval for which outputs are generated is not critical. More 
than N points can be outputed so long as exactly N consecutive points 
represent one period. In general it will be desirable to output more 
than N points so that if the first few contain transient effects they 
can be discarded. 

Ihe interface program, PLT, is used to access the data file and dis- 
play the data in either the time domain or the frequency domain. Display 
capabilities include printed listings and plotted time waveforms or fre- 
quency spectra. Two types of plot routines are available, printer type 
plots (scaled not to exceed 72 columns so that they can be displayed by 
a teletype terminal), and Calcomp plots (frequency spectra only). 


*This list can contain any of the output variables specified in the CIRCUS 
users manual [3]. 
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In the next section a specific example of using the program is 
shown, and the overall operation is further exaplined. A listing of all 
software used in the program is given in Appendix B. 

3. Example of Using the CIRCUS -Frequency Domain 
Interface Program 

Consider the simple transistor feedback amplifier shown in 
Figure 1. This amplifier, when operating in its linear region, would 
have a gain approximately equal to the ratio of the collector to emitter 
resistances. The drive has been set to a level such that the amplifier 
operates in a nonlinear region, theyby producing distortion in the out- 
put signal. A sample execution of the CIRCUS -frequency domain interface 
program for this circuit is presented to demonstrate the capabilities of 
the program. 

Before executing CIRCUS, the file to accept the time waveform data 
from CIRCUS must be declared and linked to logical unit 19. The setup 
on the Univac-1108 is as follows: 

(?ASG,A PFILE. 

(?USE 19, PFILE 

PFILE could be either a tape file or a mass storage file; as used here 
it is a Fastrand (mass storage) file. Execution of CIRCUS with a plot 
statement containing the data to be plotted will write these data on 
PFILE. An execution of CIRCUS for the simple circuit shown in Figure 1 
is shown in Figure 2. 

To use the transform and plotting program PLT, the file containing 
the CIRCUS data must be assigned to the run by the statements: 

^ASG,A PFILE. 
i^USE 19, PFILE 

The command to execute the plotting program produces several messages 
relative to the number and size of the data sets generated in CIRCUS and 
to the number of points in a transform. Figure 3 shows a sample execu- 
tion of the transform and plotting program. The first messages received 
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Figure 1. Test Circuit for the CIRCUS Frequency 
Domain Interface Program. 
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©CAT PFILE.,F2 
HEADY 

©ASG,A PFILE. 
rtEADY 

©USE 1 9,PFILE 
READY 

©XQT CAA3 • 


START 51 

©ADD CD1.DATA17 
'SINGLE STAGE TEST CIRCUIT' 


HI , 

1 , 

2, 

100. 


R2, 

3 , 

7, 

8.2E3 


K3 , 

0, 

3, 

2.2E3 


RA, 

A , 

7, 

1 *E3 


H5 , 

5, 

6 , 

68. 


H6 , 

0, 

6, 

270 . 


Cl , 

2, 

3, 

1 .E-8 


C2 , 

0, 

6, 

1 .E-6 


VS1 

, 0 

, 7 

, 12. 


T 1 , 

3 , 

A, 

S, 2N697 


SV1 

, 0 

, 1 

, 0 » , 0.5 

, 0 • , 

DEV 

ICE 

PARAMETERS 


TRANSISTOR, 2N697, 

NPN, 

RB , 

70 

. , 

HC , 1*5, 

RE, • 

PHI 

1 , 

.9, 

Nl, .37, 

A2 , 


N2 , .36, IES, A.A7E-1A, 

THETA I, 35.5 


2.E-6 


001, Al, 5A.E-12, 
30.3E-12, PH 12 , .9, 

ICS, 2.2AE-13, THETAN, 


BN, 0., 1 .E-5, l«, 68* 


BI, 0 * , 1 .E-5, 1 *, 4.7 
TCN, .001, 9.93E-10 
TCI, .001, .91 E-7 

END 

INTERVALS, .125E-6, 8.E-6 
PRINT, VN1 , VN2, VN4* IBT1, 
PLOT, VN1, VN2 , VN4, IBT1, 
HOLD FINAL CONDITIONS 
EXECUTE 
EXECUTE 

LINKS 50 


1CT1 

ICT1 


38*8, 


Figure 2. Execution of CIRCUS for Test Circuit. 
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SINGLE STAGE TEST CIKCUIT 


» 



TIME 
(NSEC ) 

VN1 

0N2 

UN4 

IBT1 

ICT1 


.0 

0 .000 

1 .676-08 

6.91 

7.491 -05 

5.094-03 


125.0 

.191 

• 148 

5.59 

4.245-04 

6.407-03 


250 *0 

.354 

.296 

3.74 

4.858-04 

8*264-03 


375.0 

.462 

.399 

2.38 

4.811 -04 

9.618-03 


500 .0 

.500 

.430 

2.19 

5.375-04 

9.806-03 


625 .0 

.462 

.414 

2.17 

3.399-04 

9.833-03 


750 .0 

.354 

.354 

2.15 

-1 .160 -04 

9.853-03 


875.0 
CUSEC ) 

.191 

.241 

2.84 

-5.401 -04 

9.156-03 


1 .000 

1 .589-08 

4 .633-02 

5.76 

-3.996-04 

6.241 -03 


1 .125 

-.191 

-.144 

8 .46 

-2.984-04 

3.538-03 


1 .250 

-.354 

-.306 

10.6 

-2.144 -04 

1 .443-03 

* 

1 .375 

-.462 

-.420 

1 1 .7 

-1 .022-04 

2.884-04 


1 .500 

-.500 

-.466 

1 1 .9 

3.461 -06 

5.021-05 


1.625 

-.462 

-.436 

1 1 .9 

6.266-05 

9.485-05 


1 .750 

-.354 

-.343 

1 1 .4 

1 .624-04 

6.221-04 


1 .875 

-.191 

-.201 

9.95 

2 .823-04 

2 .054-03 


2 .000 

-3 .1 79-08 

-2.918-02 

7 .91 

3.838-04 

4.088-03 


2 .125 

.191 

.145 

5.72 

4.551 -04 

6.282-03 


2 .250 

.354 

.296 

3.77 

4.886-04 

6.228-03 


2 .375 

.462 

.400 

2 .40 

4.794-04 

9.598-03 


2 .500 

.500 

.431 

2.19 

5.303-04 

9.806-03 


2.623 

• • 1 r 2 

• i’Ui 

2.17 

3.332-04 

9.832-03 


2 .750 

.354 

• 355 

2.15 

-1 .222-04 

9.852-03 


2 .875 

.191 

.240 

2.90 

-5.305-04 

9.102-03 


3 .000 

2 .265-07 

4 .604-02 

5.79 

-3 .955-04 

6.206-03 


3.125 

-.191 

-.144 

8 .49 

-2.976-04 

3.514-03 


3 .2 50 

-.354 

-.306 

10.6 

-2.141 -04 

1 .424-03 

* 

3 .375 

— .462 

-.419 

1 1 .7 

-1 .013-04 

2.773-04 

jt ' 

3 .500 

-.500 

- .466 

12.0 

2.678-06 

4.784-05 

* t 
iP 

3 .625 

-.462 

-.435 

1 1 .9 

6.238-05 

8.995-05 

f 

3 .750 

-.354 

-.342 

1 1 .4 

1 .616-04 

6.091 -04 

,fe 

3 .875 

-.191 

-.200 

9.96 

2.817-04 

2 .036-03 

§:■ 

4 .000 

-4.212-07 

-2.903-02 

7.93 

3.832-04 

4.070-03 

k- ; ■ 

4.125 

.191 

.145 

5.74 

4.544-04 

6.263-03 

'• 

4.250 

.354 

.296 

3.79 

4.876-04 

8.21 1 -03 


4 .375 

.462 

.400 

2.42 

4.774-04 

9.584-03 

■*$, 

4.500 

.500 

• 432 

2.19 

5.247-04 

9.805-03 


4 .625 

.462 

.415 

2.17 

3.278-04 

9.832-03 

& 

4.750 

.354 

.356 

2.15 

-1 .271 -04 

9.852-03 


4 .875 

.191 

.239 

2.94 

-5.232-04 

9.061 -03 

: ,T . 

5.000 

5.563-07 

4.583-02 

5.82 

-3.925-04 

6.178-03 

& 

5.125 

-.191 

-.144 

8.50 

-2.970-04 

3.495-03 

’k: 

5.250 

-.354 

-.306 

10.6 

-2.139-04 

1 .409-03 


5 .375 

-.462 

-.419 

11 .7 

-1 .007-04 

2.683-04 

i 1 

5.500 

-.500 

- .465 

12*0 

2.095-06 

4.604-05 


5.625 

-.462 

-.435 

11 *9 

6.217-05 

8.619-05 


5.750 

-.354 

- .342 

11 .4 

1 .609-04 

5.989-04 


5.875 

-.191 

-.200 

9.98 

2.812-04 

2.022-03 

V. 

6.000 

-8.106-07 

-2.892-02 

7.95 

3.827-04 

4.055-03 

■jl, .. 

6.125 

• 191 

• 146 

5.75 

4.538-04 

6*249-03 



Figure 2. 


(Continued) . 
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END L5 43 


'SINGLE STAGE 


TINE V/Nl 


(USED ) 


6.2 50 

.354 

6.375 

.462 

6.500 

.500 

6.625 

.462 

6.750 

.354 

6.875 

.191 

7 .000 

1 .005-06 

7.125 

-.191 

7.2 50 

-.354 

7 .375 

-.4 62 

7 .500 

-.500 

7.625 

-.462 

7.750 

-.354 

7 .875 

-.191 

8 .000 

-1 .319-06 


vN2 VN4 


.296 

3 .80 

.400 

2 .43 

.432 

2.19 

.415 

2.17 

.356 

2.15 

.239 

2.97 

4.567-02 

5.84 

-.144 

8.52 

-.306 

10.6 

-.419 

1 1 .7 

-.465 

12.0 

-.435 

1 1 .9 

-.342 

1 1 .4 

-.200 

9.99 

-2.882-02 

7 .96 


IBT1 

4.867-04 
4.759-04 
5.202-04 
3 .23S-0 A 
-l .310 -04 
-5. 176-04 
-3.901 -04 
-2.965-04 
-2.137-04 
-l .002-04 
1 .661 -06 
6.201-05 
1 .604-04 
2.808-04 
3.823-04 


END OF JOB 
END OF JOB 
ENDJOB 43 

END 8049 MLS EC 


Figure 2. (Continued). 


TEST CIUCUIT' 


ICT 1 

8.197-03 
9.573-03 
9.805-03 
9.832-03 
9.851 -03 
9.028-03 
6.157-03 
3.481 -03 
1 .397-03 
2.61 1 -04 
4.468-05 
8.331 -05 
5.908-04 
2.011-03 
4.043-03 
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<sxPT 

PL T. 



THE 

NUMPEP 

OF 

PO I N T 5 

THE 

NIMFEP 

OF 

PO I N T S 

THE 

NUMPEP 

OF 

DATA S 


PH P DATA SET * 65 

DPEI IN A TFANSFOPM = 
TS = 5 


ENTEP O FOP PPINT OF DATA SET EPOM Cl PCI'S 
ENTEP PLUS TATA SET NUMPEP EOP EFEP E CN 
ENTEP MINDS DATA SET NIWPEP FOP TIME F CN 
V 

ENTEF V FOP TIME LISTING OP DATA SET NDMPEP 


TIME 


0. OOOO 
7. 5808-07 
1 . 5000-86 
p. P500-86 
3. 8000-86 
3. 7 580-86 
A. 5080-86 

5. P500-86 
6. 8880- 06 

6.7 588- 06 
7. 5000-06 


1 . P600- 07 

8 . 7 500- 07 
1.6P50-06 
?. 37 50- 06 
3. 1 258- 06 

3.87 58- 06 
A. 6P50- 06 

5. 3750- 06 

6 . IP50-06 

6.87 50- 06 
7.6 250-06 


P. 5000- 07 

1 . 0000 - 06 

1 . 7 500-06 
p. 5000-06 
3. P500-06 
A. 0000- 06 
A. 7 500— 06 
5. 5000- 06 

6 . P500-06 

7 . 0000 - 06 
7.7500-06 


3. 7500-07 
1 . 1P50- 06 
1 .8750-06 
?. 6P50- 06 
3. 3750- 06 
A. 1P50- 06 
A. 8750- 06 
5. 6P50- 06 
6. 3750- 06 

7. 1 P 50- 06 

7. R750-06 


ENTEP 0 FOP PPINT OF DATA SET EPOM Cl PCI'S 
ENTEP PLUS DATA SET NDMPEP EOF FFEP E CN 
ENTEP MINUS DATA SET NIWPEP FOP TIME F CN 


0 

ENTEP 0 FOP TIME LISTING OP 


DATA SET NIMFEP 


3 


data set niwpep 3 


6.9060+00 
?. 1A67+88 
1 . 19 50+ 0 1 
3.77 17+00 

5. 79 A3+00 
1 . 1391 + 01 
p. 19 A5+ 00 
1 . 059 1 + 0 1 
7 .9 A52+88 
P. 1 A86+00 
1 . 19 55+01 


5. 59 3A+00 
p. 8 A A A+ 00 
1 . 19 0 5+0! 
P. AO 1 7+00 
8 . A 8 59+ 00 
9.96AP+00 
p. 1677+00 
1 . 17 32+01 
5.751 A+00 
?.97 1 7+00 
1.19 17+01 


3 . 7 36 ?+ 00 
5.7588+00 
1 . 1 378+0 1 
?. 19A1 + 00 
1 . 0576+01 
7.930 A+00 
?. 1 A 8 1 + 00 
1 . 19 5A+0 1 
3.8028 + 00 
5.8A35+00 
1 . 1 A09+0 1 


?. 38 16+00 
8 . A6P0+00 
9.9 A65+00 
?. 1 675+00 
1 . 1723+01 
5.7368+00 
p. 9 39 A+ 00 
1 . 19 1 A+01 
?. A27P+00 
8 . 519 5+00 
9.989 1+00 

Cl PCI'S 
E CN 
F CN 


ENTEP 0 FOF PPINT 
ENTEP PLUS DATA 
FNTEP MINUS DATA 


OF DATA SET FPOM 
PET NUMPEP FOP EPEP 
PET NUMPEP FOP TIME 


3 


3P 


5. 0000-07 
1 . P500-06 
p. 0000- 06 
P. 7500-06 
3. 5000-06 
A. P500-06 

5. 0000- 06 

5. 7 500-06 

6 . 5000-06 

7. ? 500- 06 

8. 0000-06 


?. 1936+00 
1.0557+01 
7.911 7+00 
P. 1A75+00 
1 . 19 52+01 
3. 7890+00 
5. 8 P 19+00 
1 . 1 AO 1 + 0 1 
?. 19A8+00 
1 . 060 3+0 1 
7.9569+00 


6 . P 500- 07 
1 . 37 50- 06 
P. 1 P50- 06 
p. 87 50- 06 
3.6P50-06 
A. 3750- 06 
5. 1 P 50- 06 
5 . 8750- 06 
6.6P50- 06 
7. 3750- 06 


?. I67P+00 
1 . 1 7 1 ?+ 0 1 

5.7 18 3+00 

p. 8980+00 
1.19 10+01 
P. A 1 59+ 00 
8. 50A7+00 
9.9782+00 
?. 1679+00 
1. 17 39+01 


Figure 3. Execution of Interface Program. 
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EN1EF IFTAFT 
39 

THE PE.FI 0 V OF THE TIME FUNCTION = 4 . <*0(*<*- (*6 SEC 


fcNTfcP 100 kOP PHINT/ 010 POP CALCOMP PLOT/ OP 00 1 FOP *IT V PLOT 
1 00 


LIMP 


kNTkP 

kNTFP 

3 


PkAL 


IMAG 


DP 


1 

-P. 7697-0P 

0. 0000 

-31.15 

? 

-4.9 144-04 

6.6766- 04 

-61.63 

3 

P.009 1-0P 

7.9110-03 

- 33. 31 

4 

-Q . 1 36 5-05 

-6.3436-04 

-63.86 

b 

- P. 3060- 0P 

9. 5P55- 03 

- 3?. 06 

6 

-6.6337-04 

P. 3? 14- 04 

-63. 06 

7 

p.9 1 55- 0? 

- 4. 58 56- 0P 

-P5. 30 

8 

- 1 . Q 304- 04 

- 1 . 16P9-03 

- 58. 57 

9 

-7.8478-0P 

- 1 . 64 30- 03 

-PP. 10 

1 <A 

- 5. 3355-04 

3.96 14-04 

-63. 55 

1 i 

P. 309 1-01 

1 . PI 46- 0? 

- IP. 7? 

IP 

1 . 09 4 3- 0 3 

- 1. 6PP9-03 

- 54. 17 

1 3 

- 1 . 307 0-01 

3. 1 161-0? 

-17.43 

1 4 

- 1.7313-03 

-8.69 59 - 05 

- 55. PP 

1 b 

l . 47 38+00 

-P. P6P5+00 

«. 63 

! 6 

-K.8Qt>3-04 

- 3. 4519-03 

-48.96 

17 

6.9555+00 

0. 0000 

16.85 

IP 

-8.8963- 04 

3. 4519-03 

-48.96 

19 

1 . 47 38+00 

P. P6P5+00 

8.63 

9V 

- 1. 7313-03 

8. 69 59 - 05 

- 55. PP 

PI 

- 1 . 307 0-01 

-3. 1 161-0P 

- 17.43 

92 

1. 0943-03 

1 . 6PP9-03 

- 54. 17 

23 

?. 309 1-01 

- 1. PI 46-0? 

- IP. 7P 

9 A 

-5. 3355-04 

- 3.9614-04 

-63. 55 

2 5 

-7.8478-0P 

1 . 6430-03 

-PP. 10 

P6 

- 1 .9 304- 04 

1. 16P9-03 

- 58. 57 

PI 

P.9 1 55- 0P 

4. 58 56-0P 

-P5. 30 

2P 

-6. 6337-04 

-P. 3014-04 

-63. 06 

?9 

-P. 3060- 0P 

-9. 5P55-03 

- 3P. 06 

30 

-9. 1 365- 05 

6. 3436-04 

-63.«6 

3 1 

P. 009 1 — 0P 

-7.9110-03 

- 33. 31 

3? 

-4.9 144-04 

-6. 6766-04 

-61.63 

riD 

n h.'»d DPTMT Ok LATA SkT P POM CIFCOS 


PLUS L'ATA SkT NUMPkP POP P PkO P CN 
MINUS L'ATA SkT NlMPkP POP TIMfc P CN 


Figure 3. (Continued). 
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ENTEP ISTAPT 

3 ENTEP 100 FOP PPI N T» 010 FOP CALCOMP PLOT* OP 00 I FOP TTv PLOT 
1PP 



LINE 

PEAL 

I MAG 

L'P 


1 

9.964?+ PP 

0. 0000 

19.97 


? 

7.9 304+ PP 


17.99 


3 

5. 7368+00 

O' • 

1 5. 17 


4 

3. 789P+0P 

04 0000 

1 1 . 57 


5 

P. 4159+00 

0 4 

7.66 


6 

?. 194 5+ 00 

0 • 00 00 

6.8 3 


7 

P. 1 677+ PP 

0. 0000 

6. 7P 


8 

?. 1 48 l+PP 

0 • 0 0 00 

6. 64 


9 

P.9 39 4+ PI PI 

0. 0000 

9. 37 


If 

5.8? 19+ Pif 

0.0 000 

1 5. 30 


1 1 

8. b047+PP 

0 • 0000 

18. 59 

* 

1? 

1 . P59 1+P1 

0. 0000 

P0. 50 


13 

1. 17 3P+P1 

0. 0000 

21. 39 


14 

1. 19 54+01 

0 • 0000 

PI. 55 


1 5 

1.19 14+01 

0 • 0000 

PI. 5P 


16 

1. 14P1+01 

0* 0000 

PI. 14 


17 

9.978 P+ 0 0 

0 • 0000 

19.98 

*;• 

1« 

7.945 P+ 0 0 

0 • 00 00 

1R. 00 


IP 

5.7 514+P0 

0#00 00 

15.P0 


77 

3.RPPR+00 

04 0000 

1 1.60 

H 

PI 

?. 4P7P+ 00 

04 0000 

7.70 


PP 

P. 194R+PP 

04 0000 

6.R3 

w 

?3 

P. 1679+00 

04 0000 

6.72 

iv-' 

?4 

P. 1 4R6+00 

04 0000 

6.64 

1* 

?5 

P. 97 1 7+0P 

0. 0000 

9.46 

'§{ 

?6 

5.R435+00 

0 4 0000 

15.33 

jfc 

> i ■ 

27 

R. 5195+00 

0# 0000 

1R.61 

£ 

PH 

1.0603+01 

04 0000 

?P. 51 


?9 

1. 17 39+P1 

04 0 000 

PI. 39 

M;- 

30 

1. 1955+01 

0 4 0000 

PI . 55 


31 

1. 1917+01 

04 0000 

PI . 5P 


3P 

1. 1409+01 

0. 0000 

PI. 15 


ENTEP 

0 FOP PPINT OF 

L AI A SLT FFOM Cl PCI'S 


fi 

ENTER 

PLUS PAT A SET 

NtWFfcP FOP FPFP FCN 


■wii< 

ENTEP 

MINUS DATA SET 

NIWPEP FOP TIME F CN 



Figure 3. (Continued) 
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after execution begins indicate a data set size of 65, the number of 
points in a transform to be 32, and the number of data sets to be 5 
(these correspond to the 5 variables in the PLOT statement, see Figure 
2). A message set then follows which requests the entry of an integer, 
a zero for print of a data set, a positive data set number for output 
of the frequency function, or a negative data set number for output of 
the time function. 

The example shows the entry of a zero followed by another zero in 
answer to the next question, thus indicating data set number 0 (the time 
listing). This is followed by the 65 entries from CIRCUS making up data 
set 0. After the printing of these data is complete, the original ques- 
tions pertaining to the data desired are repeated as shown in the example 
of Figure 3. This time a print of data set 3, the node 4 output volt- 
age, from CIRCUS was requested. 

The next response was a +3 requesting the frequency domain repre- 
sentation of data set number 3. The program then requests the value of 
ISTART, an integer specifying the starting index for the 32 values to be 
used in the transform. Care should be exercised to insure that ISTART 
is never greater than the number of points in the data set minus the num- 
ber of points required for a transform; otherwise, the transform data 
will be taken from two adjacent data sets (or partially from an adjacent 
storage area) and the transform will be meaningless. 

After entry of the value of ISTART, the period of the time function 
is displayed. This display occurs only on the first call for either a 

time function or a frequency function. 

Next a statement telative to the type output desired is displayed. 
The code for response to this statement is 100 for print, 010 for Cal- 
comp plot, and 001 for teletype plot. The example shows a request for 
print, followed by a printing of the 32 values of the frequency function 
obtained from the original data set number 3 starting at array element 
number 32. The three output types may be called individually, in pairs, 
or all at one time. For example, responding with 111 would generate 
printed output, a Calcoop plot, and a teletype plot of the specified data 
set. The response 101 would produce printed output and a teletype plot. 
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After execution of the print command in the example of Figure 3, 
the program cycles back to the first set of questions and asks again for 
the entry of a 0, or plus or minus a data set number. The response 
in the example of -3 indicates a request for the time waveform of 
data set number 3. Next, in response to the questions generated by the 

program, ISTART is specified as 32, and an output type of 100 (printed 

output) is entered by entry of the number 100. This is followed by a 
print of the 32 data points which would be input to the FFT had a fre- 
quency function been called. These 32 points are from data set 3, start- 
ing at element number 32 and ending at element number 63. This provision 
f 0; - selecting any 32 consecutive points desired from a 65 point data set 
allows the first part of the data to be skipped when it contains a tran- 
sient response. 

After execution of the given command che program again responds with 
a request for output data type desired. The example of Figure 3 shows a 
request for the frequency function of data set 3 starting with data set 
element number 32 and for a teletype plot of these data by entry of 001 

to the request for output type. The program next asks for entry of FLO, ? 

the lowest frequency desired in the plot. This can be any frequency com- | 

patible with the size of the stored frequency function. Entry of this | 

frequency is followed by a request for entry of FH1 , the highest frequency | 

desired in the teletype plot. Entry of these two frequencies satisfies j 

the requirements of the program and a plot is produced on the teletype \ 

as shown in the example of Figure 4. The ordinate of the plot is dis- 
played horizontally and is calibrated in decibels with automatic scaling. 

The abscissa of the plot appears vertically (along the teletype paper) 

{ 

and can be any length. This axis along the paper is the frequency axis 

extending from the FLO to FHI specified. 

Completion of the teletype plot again generates the request for the 
output data desired which in the example is answered by a -3, a request 
for the time function of data set number 3. Again the request for ISTART 
is responded to with a 32 and a teletype plot is requested by entry of 
001. Since a time plot has been specified, the program requests that the 
starting and stopping value of the index of the data array be entered. 


s 
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ENTFF I STAFT 
3 ? 

THE PEPI 0 E> OF THE TIME FINCTION * 4 . 0000-06 SEC 

ENTEP 100 FOP PPINT# 010 FOP CALCOMP PLOT/ OP 001 POP TTV PLOT 
P01 

ENTEP FLO 
- 3 . 7 5 E 6 
ENTEP FHI 
3 . 7 5 F 6 


i 


NSI 7 F = 31 


FPEOUENCY (MHZ) 

DECIPELS 




-30 - ?0 

- 10 



10 

I + I — 

— + I 

--I 

» - + ~ - 

-- I-- 


-3.750 ! 

-3. 500 I 
-3.250 I 
-3.000 I 
-P.750 I 

-P.500 I 

-P.P50 I 
-P.000 I-- 

- 1 • 7 50 I 

- 1 . 500 I 

- 1 . PSP I 

- 1 . 000 I 

-.750 I 

-.500 1 

- • P50 I 

• OOP 1"**“*"*““*"""*"""*""*""" 

• P50 I 

. 5 PiP> I 

• 7 5fl I 

1.000 1 

I.P50 I 

1.500 I 

1.7 50 I 

P.000 I 

P.P50 I 

P. 500 1 

P.750 I 
3.000 I 
3. P50 I 

3. 500 I 
3.750 I 

FPEP 

ENTEP 0 FOP PFINT OF DATA SET FPOM Cl PCI'S 
ENTEP PLUS DATA SET NUMBEF FOP FPEG F CN 
ENTEP MINUS DATA SET NUME*EH FOP TIME FCN 
-3 


P0 

- I 


Figure 4. Spectrum of Data Set 3. 
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The data array available for the time plot is contained in an array 
having the length of the transform array specified (in this case 32). 

Any stopping value up to the maximum length of the array may be speci- 
fied (in this case 32). Ihe entry shown is 32. The next request by 
the program is that for NJUMP which specifies the number of points 
skipped in the array between points plotted. This allows skipping 
points between plotted points and is a very useful feature when large 
data arrays are processed. In the example shown the transform data 
array size is small and the value of NJUMP is entered as 1. The out- 
put is shown in Figure 5 which displays the time waveform of the ampli- 
fier output. 

The last request illustrated is that for a Calcomp plot. The 
Calcomp plotting routines and methods used for setting up the plot 
files are probably unique with the Georgia Tech 1108 computer. Simi- 
lar routines and specific control directives should allow the use of 
the plotting routines with a minimum of changes at other installations. 
The Calcomp frequency plot routine contains as its last instruction a 
write statement to write a message "Plot Complete" indicating that the 
program has processed the plot routine. The command to generate a Cal- 
comp plot of the frequency function is a plus data set number (3 in this 
case) followed by ISTART (32) followed by 010. These instructions are 
shown in Figure 6, and the resulting plot in Figure 7. 

To exit the program, the command @E0F is given. The remaining in- 
formation relative to the Calcomp plot and the time to execute the run 
is generated by the computer. 


B. Use of ECAP for Transfer Function Evaluation 

One of the difficulties frequently encountered in the use of CIRCUS 
was that of obtaining an accurate representation of a circuit for analy- 
sis. If the circuit description presented to CIRCUS is not an accurate 
representation of the actual circuit, then the results produced by CIRCUS 
will not accurately represent the circuit response. An example of this 
type of difficulty is that of entering an interstage bandpass filter 
which contains mutual inductance and is situated between two active 
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ENTFP I STAPT 

3 P 

ENTEP IAA FOP PPINT* AIA FOP CALCOMP PLOT* OP AA 1 FOP TTY PLOT 
AA 1 

APPAY SI 7 E = 3P N STOP MUST BE EQUAL TO UP LESS THAN THIS VALUE 


ENTEP N ST APT 

1 

ENTEP N STOP 

3? 

ENTEP N.JIMP 

1 


+ 

AMPLITUDES MIN . AAAA * MAX . 1196+AP VOLTS 


+ 


0 

I- 

1 I 

2 I 

3 I 

4 I 

5 I 

6 I 

7 I 

8 ! 

Q I 

1 0 I 
1 1 I 
1? I 
1 3 I 

14 I 

15 I 

16 I 

17 I 

18 I 

19 I 

20 I 

21 I 

22 I 

23 I 

24 I 

25 I 

26 I 

27 I 

28 I 

29 I 

30 I 

31 I 

32 I 
N 


1 • 2 *3 • 4 *5 

I i i i I 


* 


5k 


5k 

★ 

* 

* 


* 


* 


sk 


>k 

Jk 

* 


sk 


>k 


6 

I 


5k 


5k 


ENTEP A FOP PPINT OF DATA SET FPOM Cl PLUS 
ENTEP PLUS DATA SET NIWPEP FOP FPEQ F CN 
ENTEP MINI'S TATA SET NIMBFP FOP TIME F CN 


7 

I 


* 


* 


3 


8 .9 1 . A 

I I I 

5k 


5k 

5k 

5k 

♦ 

5k 


5k 


5k 


ik 


5k 

m 


ik 


Figure 5. Time Waveform of Data Set 3. 
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ENTER ISTAPT 

3 ihe PERIOD OF THE TIME FUNCTION = 4.O00O-06 rec 

ENTER lOO FOP PRINT.* 010 FOP CALCOMP PLOT* OP 001 FOP TTV PLOT 
^EN TFP THE HIGHEST DESIRED FPEO IN THE SPECTPlM* FMAX 
3.7 5E6 


PLOT COMPLETED c 

ENTEP <* FO P PPINT OF LATA SET FPOM Cl PCI'. 

ENTEP PLUS LATA SET NLWBEP FOP FPEP F CN 
ENTEP MINUS DATA SET NIWEEF FOP TIME F CN 
#E0 F 

PLOT ?• 7 MIN 1*3 FT O^F 11 ?<*7 ? IQ 1 I 04 

tN p 1 34P ML SE C 



I 

't' 

| 


Figure 6. Requeet for Calconp Plot. 
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(FREQUENCY - FCENTER) OlViDED BY FSCflLE . (HZ) 

FCENTER - 0. FSCfiLE - 1000000* 

Figure 7. Calcomp Plot of the Frequency Spectrum of the Time Waveform 
Shown in Figure 5. 


(I 


devices. Such filters usually have bandwidths wide enough to pass the 
significant sidebands of a signal of interest after allowing for such 
factors as frequency instability and temperature effects. For CIRCUS 
to provide an accurate representation of the circuit response, the trans- 
fer function of the filter must meet these bandwidth requirements. Since 
CIRCUS determines the circuit response by solving for its time response 
(which yields only a time waveform), the transfer characteristics of the 
bandpass filter cannot be readily determined from CIRCUS calculations. 

In a search for better methods of coping with such problems, other 
analysis programs were investigated. Those having an a-c analysis cap- 
ability in the frequency domain were of particular interest since the 
circuit description could be entered in a manner similar to that used 
in CIRCUS. If the computed circuit response was not as desired, adjust- 
ments in circuit parameters could be made and the transfer function eval- 
uated again. Thus the variation of the transfer function with variation 
in the value of circuit parameters could be obtained. 

With such an a-c analysis capability, much of the guesswork of 
entering a circuit into CIRCUS could be removed. This is not to say 
that the response of a filter, for example, is independent of the devices 
preceding and following it, but that values of circuit parameters could 
be determined to yield the desired filter response. These could then 
be adjusted for drive and load impedances, possibly by using CIRCUS to 
calculate the impedances presented by the circuit external to the filter. 

The analysis program ECAP was selected to be used for the a-c analy- 
sis calculations. This selection was based largely on the availability 
of ECAP on the Georgia Tech Univac-1108 computer. 

To illustrate the use of ECAP, a fairly complex circuit with a known 
response was selected. This circuit was the CCS telemetry bandpass fil- 
ter for which calculated and experimentally determined responses have 
been previously obtained [1]. A circuit diagram of this filter is shown 
in Figure 8 with the input transformer replaced by its "T" equivalent 
circuit. The underlined numbers in the figure show the node numbers 
used for entry of the circuit into ECAP. The branches are shown in the 
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figure as numbers preceded by a "B" and are from left to right and 
from top to bottom. The circuit description in the format required by 
ECAP is shown in Figure 9. The circuit response calculated by ECAP 
agreed very closely with that obtained earlier, and is shown in Figure 10 
by the "x" points superimposed on the response calculated previously [1]. 

A sample of the output produced by ECAP is shown in Figure 11. The 
output of ECAP is not easily interpreted. ECAP essentially re-executes 
for each new frequency and therefore produces an excessive amount of 
superfluous output information. A much easier to interpret output for- 
mat would result if the output data were displayed in tabular form. 

The ability of ECAP to solve for a transfer function of a circuit 
has been demonstrated and a familiarity with the user program interface 
was obtained. 


C. Conclusions and Recommendations 

An interface with the time domain program CIRCUS has been constructed 
which allows the frequency domain representation of time waveforms gene- 
rated with CIRCUS to be produced. These frequency domain representations 
of signals are generated by use of the FFT algorithm. In addition, the 
capability has been provided for easily obtaining listings of the data 
sets obtained from CIRCUS, or for generating Calcomp or teletype, time 
or frequency plots of portions of these data sets* 

The a-c analysis capabilities of ECAP have been investigated and 
found to be a useful tool for accurate preparation of data for input to 
programs such as CIRCUS. This a-c analysis capability is particularly 
useful in the evaluation of transfer functions. 

It is recommended that the capabilities of the two programs, CIRCUS 
and ECAP (preferably the latest versions of these programs) be interfaced 
such that nonlinear circul ts could be handled with CIRCUS and linear cir- 
cuits with ECAP. The link between the frequency and time domains would 
be the fast Fourier transform algorithm. Part of this interface has al- 
ready been established with techniques discussed in this section. There 
remains a need to interface the a-c analysis capabilities of ECAP with 
the frequency domain data obtained from CIRCUS. A study should be under- 
taken to determine an efficient method of providing this interface. 
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Figure 9. Circuit Description in ECAP Format. 
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III. MATHEMATICAL MODELING OF COMMUNICATION SYSTEMS 
ON A BLOCK BASIS 

A. General Description 

The Frequency and Time Circuit Analysis Technique (FATCAT) is a 
computer implemented program for analyzing communications circuits and 
is designed on a block modeling approach. Circuits are represented as 
a linear collection of sub-assemblies and FATCAT provides a model block 
for a variety of sub-assemblies including signal sources, filters, demodu- 
lators, amplifiers, limiters, mixers, etc. The program is designed for 
use with either remote terminal or batch processing mode. For remote 
terminal operation the program is conversational in nature and provides 
considerable flexibility to the user. Input statements to specify the 
circuit configuration, to direct processing, and to direct output are 
given in alphanumeric codes along with numeric specification of para- 
meters. The input formats were designed to be relatively simple and easy 

to use . 

FATCAT is designed for steady-state analysis of circuits in which * 

the signal flow is sequential; no provision is presently included for 
feedback loops. The model is designed to operate in both the time do- j 

main and the frequency domain. The signal being processed is stored in jj 

a complex array as either a frequency domain representation of the sig- 
nal or one complete period of a time wave form. Transition from one j 

state to another is made using a subroutine which performs a fast Fourier ; 

transform. Each model block was developed in whichever domain was most 
convenient for modeling that block. When calling for signal processing 
through any block, the domain representation of the signal which current- 
ly exists is checked and, if necessary, the conversion is automatically 
made. The same automatic conversion is available on output calls; this 
permits examining either the time waveform or the frequency spectrum at 
the output of any block. Outputs include both printed and plotted values 
of the time waveform or the frequency spectrum. 

The entire program was developed and implemented on a Univac-1108 at 
Georgia Tech. Operation and command structure of the program along with 


29 


4 


I! 


examples of using the program are given in the following sections. A 
brief description of each software unit along with a program listing is 
given in Appendix C. 

After development of the program was completed to its present form, 
suitable modifications were made to implement it on a SIGMA-5 computer 
at Marshall Space Flight Center. Appendix D gives program listings of 
those routines that were changed for use on the SIGMA-5. 

B. Program Control 

Operation of the program is accomplished by giving a sequence of 
defined commands, each command being given basically with an alphanumeric 
string of 6 or less characters. In most cases other data, usually nu- 
meric will follow the basic command. 

The set of commands can be broken into several categories (1) input 
(block specification) commands, (2) control commands, and (3) special 
commands . 

After starting execution of the program, the first input command 
will either be specification of the first block of the circuit, or one 
of two special conmands. These two special conmands are (1) the ability j 

to list the command structure of the program, and (2) the ability to ! 

list the input formats for all input commands. These two commands are 

provided as a convenience to the remote terminal user and are described j 

in detail in Section 1IIB3. 

r 

1. Block Input Commands X 

Each circuit block is specified with a line entry composed of 
the alphanumeric code for that block, followed by the parameters of the j 

block. The entry items on a line are separated by commas; no comma is 
used after the last item. Each block has a specific entry format and 
these are described below. 

The blocks of a circuit are entered sequentially in the order that 
they occur in the circuit. Thereafter they are identified by number, 
i.e., the first block entered is block number 1, the second is block 

f 

r * 
j. 
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number 2, etc. The block numbers are important in directing processing 
and outputs. 

Flexibility is provided for the remote terminal user in entering 
the blocks. All blocks comprising a circuit may be entered before any 
processing takes place, or only part of the blocks may be entered and 
processed through before the remaining blocks are entered. The user has 
the option of adding new blocks to the end of the input string at any 
time (provided the block count does not exceed the maximum allowable num- 
ber of blocks, currently set at 20). One additional restriction is that 
if a block containing a source frequency (a mixer with a local oscilla- 
tor, for example) is to be added after processing has begun, the source 
frequency must be chosen to be periodic in the time interval represented 
by the time function stored in the data array. This time interval (PERIOD) 
is printed out at the beginning of processing so that the choice of a 
suitable source frequency should be relatively easy. 


a. Sources 

Two signal source blocks are provided. The first simulates a 
signal generator and is called with 


where 


SIGGEN, FO, FMOD, AM, PM, FM, A 

FO = carrier frequency (Hz), 

FMOD = modulation frequency (Hz), 

AM = percent AM modulation, 

PM = peak phase deviation (radians), 

FM = peak frequency deviation (Hz), and 
A = peak amplitude (volts). 


The signal generator block provides only sine wave modulation. Com 
binations of AM-PM or AM-FM modulation may be specified, but not PM-FM 
(either PM or FM must be zero). As an example, a signal generator with 
a carrier frequency of 1.0 MHz, a modulation frequency of 10.0 kHz, 507. 
AM modulation plus FM modulation with 0 - 2, and a peak amplitude of one 
volt would be entered as SIGGEN, 1.E6, 1.E4, 50., 0., 2.E4, 1. 
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A call for processing this block will generate the appropriate time 
waveform and store it in the data array. 

The second signal source is a flat spectrum (impulse time function) 

generator specified by 

FLATSP, AMP, DELF, N 




i 
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where 

AMP = amplitude of spectrum lines, 

DELF = frequency separation of spectral lines (Hz), and 
N = array size. 

The user specified values of DELF and N determine the overall frequency 
spread of the spectrum being produced and hence determines the period of 
the impulse function being represented. When using FLATSP, the usual 
computation of the array size and period is not used. The actual value 
of N that is used will be a power of 2; if the input N is not a power of 
2, it will be automatically raised to meet this requirement. (Example: 
if N = 500 is specified, N = 512 will be used.) 


b. Filters 

Eleven different filter blocks are included which simulate the 
action of Buttorworth, Tchebysheff, and synchronously tuned filters in 
low pass, high pass, bandpass, and (except for sync tuned) band stop con 
figurations. Since the input statements are similar, the commands will 
be grouped by types and the inputs for Butterworth, Tchebysheff, and 
synchronously tuned filters of each type given in that order. The com- 
mands for low pass filters are: 

BWLOWP, FC, NR 
CHLOWP, FC, NR, EPSDB 
SYNLP, FC, NR 

where 

FC * corner frequency (Hz), 

NR » number of filter sections, and 
EPSDB * Tchebysheff ripple factor in decibels. 





32 


•* 








i 

f 


The inputs for high pass filters are: 

BWHIP. FC, NR 
CHH1P, FC, NR, EPSDB 
SYNHP, FC, NR 

where the parameters are identical to those defined for low pass filters. 
The inputs for bandpass filters are: 

BWBNDP , FO, BW, NR 
CHBNDP , FO, BW, NR, EPSDB 
SYNBP, FO, BW, NR 

where 

FO = center frequency (Hz), 

BW = total bandwidth (Hz), (3 dB bandwidth for synchronous and 
Butterworth filters and ripple amplitude bandwidth for 
Tchebyscheff filters), 

NR » number of filter sections, and 

EPSDB = Tchebysheff ripple factor in decibels. 

The inputs for band stop filters are: 

BWBSTP, FO, BW, NR 
CHBSTP, FO, BW, NR, EPSDB 

and the parameters are identical to those for band pass filters. 


with 


where 


c. Demodulators 

Demodulators for AM, FM, and PM are specified, respectively 

AMDEMO, FO 
FMDEMO, FO 
PHDEMO , FO 

FO = center frequency (Hz). 


d. Other Blocks 

In addition to the blocks categorized above, several other 
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blocks are provided. An amplifier is specified with 
AMP, Gain 

where 

Gain = amplifier gain in dB. 

A limiter is specified with 

LIM, CL, CH, GL 

where 

CL = low clipping level (volts), 

CH = high clipping level (volts), and 
GL = limiter gain (volts/volt). 

A frequency multiplier (a wide band harmonic generator) is specified 

with 

FRQMUL 

(no parameters necessary). 

An ideal multiplier can be specified with 
IDIMUL, ALO, FLO 

where 

ALO = amplitude of the LO signal (volts), and 
FLO = frequency of LO signal (Hz). 

2. Control Commands 

A number of control commands are used to direct processing and 
to generate output. Processing of the signal to the output of any block 
is initiated with the command 

BLOCK, N 

where N is the number of the block. Note that processing occurs only 
when this command is given; entry of a block specification does not 

cause processing through that block. 

Signal processing is non-reversible; if processing to the output 
of a given block has been completed, the outputs of earlier blocks are 
no longer available. Error detecting features are included in the pro- 
gram so that a BLOCK call with an N that is wrong will not upset the 
computations . 
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The command* 

PRIME FACTORS 

will list the prime factors of all source frequencies. These can be 
helpful when it is necessary to adjust one or more frequencies slightly 
so that all source frequencies will be periodic cn an interval of rea- 
sonable size. A periodic interval that is too long will require an array 
of excessive size. 

The command 

END OF JOB 

will terminate a run. 

Other control conmands are used to generate output. Outputs include 
printed listing, printer type plots, and plotting equipment outputs of 
both the time function and the frequency function. When any of the out- 
put commands are given, the output data is from the block output where 
processing currently stands. For example, in a six-block circuit, if 
processing through block 3 has been effected by the command BLOCK 3, then 
any output generated will be for the output of block 3. To examine the 
signal at the output of block 5, it will first be necessary to process 
the signal with BLOCK, 5, and then call for the output. 

For a printed listing of the time function, the command is 

PRINTT, NSTART, NSTOP 

where 

NSTART = the starting array index, and 

NSTOP = the final array index (net to exceed the array size 
being used) . 

In conversational mode, an alternative command is 
PRINTT 

following which the machine will ask for the starting and stopping indices. 


*Dnly the first six characters of any command are important, and spaces 
(blanks) are ignored. This command can be abbreviated to PRIMEF. 
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For a printed listing of the frequency function the command is 
PRINTF, FLO, FHI 

where 

FLO = low frequency limit (Hz), and 
FHI = high frequency limit (Hz) . 

Alternatively, the command 

PRINTF 

will produce questions asking for FLO and FHI, following which the same 
output will be generated. The use of these parameters permits only the 
portion of t e spectrum of interest to be printed. 

The commands for plotted output utilize the same parameter forms 
described above, and when operating from a remote terminal the commands 
can be given with or without the parameters. When the parameters are 
omitted from the command, the computer will ask for them. 

Printer plots (line printer for batch processing, teletype plots for 
operation from a remote teletype terminal) of the time function are ob- 
tained with 

TPLOTT, NSTART, NSTOP, NJUMP 

where 

NJUMP = interval between plotted points and the other para- 
meters are the same as for PRINTT. 

Printer plots of the frequency spectrum are generated by 
TPLOTF, FLO, FHI 

where FLO and FHI are the same as for PRINTF. 

Plotter outputs are similarly generated with 
CP LOTT, NSTART, NSTOP 
for the time function, and 
CPLOTF, FLO, FHI 
for the frequency function. 
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3. Special Commands 

Three special commands are included for use in connection with 
establishing the size of the data array. Initial processing of any cir- 
cuit containing source frequencies (such as those originating in a signal 
generator or the local oscillator of a mixer) will start with a determina- 
tion of the smallest time increment on which all of the source frequencies 
are periodic, and a calculation of an array (sample) size that meets the 
Nyquist criterion for all frequencies. The calculated period and array 
size will then be printed out and the user asked if the size is satisfac- 
tory. At this point the program is positioned at a special input posi- 
tion, at which only four commands will be recognized: 

YES 

NO 

N, NSIZE 

END OF JOB 

where NSIZE is an integer which specifies the size of the array. If the 
answer is NO, the program will ask for input of an integer for NSIZE. 

For either method of entering NSIZE, the input integer will be adjusted 

to a power of 2 by increasing it if necessary, and the new values of N, 

IGAM, and DELTA-T will be printed out. The question of whether these 
values are satisfactory will then be repeated. Thus this section of the 

program is a loop and will be exited only when a YES response (or an END 

OF JOB) is given. Processing of the signal to the block output designated 
will continue after a YES is entered provided N lies within acceptable 
limits. To be acceptable, N must be at least as large as the value ini- 
tially computed in order to meet the Nyquist criterion; at the same time, 
it must be no larger than the size of th s main data array declared in the 
main program. Since cases can occur where both these conditions cannot 
be met, recognition of END OF JOB has been included at this point to per- 
mit a normal termination of the run. Any attempt to continue a run with 
N outside the required limits will cause an error termination. 

Two other special commands are included as an aid to the remote 
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terminal user. The command* 


» 


LIST COMMANDS 

followed by a comma and one of several second words will produce a list- 
ing of part or all of the FATCAT commands. Permissible second words, and 
their effect are: 

LIST COMMANDS - print instructions for using list commands, 

ALL - list all operating commands, 

SOURCES - list commands for inserting source blocks, 

FILTERS - list comnands for inserting filters, 

DEMODULATORS - list commands for inserting demodulators, 

MISC - list other block input commands, and 
CONTROL COMMANDS - list control commands. 

The above instruction will only list the command word for each in- 
struction along with identifying information. Input data formats will 
not be listed. If the user wants information on the input data format 
for any block specification command, the command 

INPUT FORMAT 

followed by a comma and the name of any block specification command will 
produce a listing of the complete command with the required parameters, 
and will identify the meaning of each parameter. For example, the com- 
mand 

INPUT FORMAT, BWBNDP 

w iH list the complete input instructions for a Butterworth band pass 

filter. 

The cotmnand 

INPUT FORMAT, ALL 

will list the formats for all of the block specification commands. 

*May be abbreviated to L1STC0. 
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C. Simulation Examples using FATCAT 

As examples of the use of the FATCAT program, two sample problems 
will be presented. Both examples represent receiving systems, the tirst 
being an amplitude modulation receiver and the second a frequency modu- 
lated receiver. 

1. AM Receiver 

A block diagram of an amplitude modulated superheterodyne re- 
ceiver is shown in Figure 12. Figure 12a shows the actual configuration 
of the receiver, while Figure 12b shows the way that the receiver would 
presently be modeled using FATCAT. The signal source is presently repre- 
sented by a signal generator which produces single tone modulated radio 
frequency signals. In the final version of FATCAT as envisioned, this 
signal source could be modeled as an actual transmitter with modulators, 
frequency multipliers and amplifiers. FATCAT does not presently include 
a propagation model or antenna models. Therefore the output of the sig- 
nal generator represents the input signal to the receiver. The control 
program is structured so that additional model blocks such as those for 
transmitters, propagation, antennas, and many more, can be easily added 
to the model blocks already included in the program. 

The first step in the analysis of the amplitude modulated receiver 
is a command for the execution of the FATCAT program. Figure 13 shows 
the command for execution of the program contained in file "R" as abso- 
lute element "U". Following this command the program responds with the 
word "START", which indicates that the program is ready for input com- 
mands. Following the sequence of blocks shown in Figure 12b the first 
input instruction given is that for an AM signal generator with a carrier 
frequency of 10 MHz, a modulation frequency of 10 kHz, 50 percent ampli- 
tude modulation, no phase or frequency modulation, and a carrier peak 
amplitude of 1 microvolt. This is the "SIGGEN" command shown in Figure 
13. Next the conmand, "BLOCK, 1" is given which requests the program 
to process the submitted instruction and arrive at the output of block 1. 
After the first conmand to process to the output of a block for which 
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Figure 13. Start of FATCAT Run for AM Receiver, 
Showing Spectrum at the Output of 
the Signal Generator. 
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input data has bean submitted (which can be a single block or multiple 
blocks) the output information shown in Figure 13 is given. This infor- 
mation contains PERIOD which gives the period of the time function, 

DELTA-F, which is the spacing of possible spectral component in the fre- 
quency domain, N which is the array size needed to meet the minimum 
Nyquist requirements, IGAM the exponent to which 2 is raised to give N, 
and DELTA-T the sampling interval in the time domain. After the quanti- 
ties have been printed the program asks if these values are satisfactory. 
Response to this question can be "YES" or "NO". If a "YES" is given pro- 
cessing proceeds to the output of the block specified. In the example of 
Figure 2 processing is completed through block 1. If "NO" had been given 
in response to the above question the program would have responded with 
"ENTER N, VALUE" which allows a new value of N to be inserted. This 
provision allows the number of samples per cycle of the highest frequency 
in the input signals to be increased. This is especially useful for plot- 
ting since just over two samples per cycle, which meets the Nyquist cri- 
terion, is not adequate for detailed plotting of the time function. 

Next it is desirable to observe the output of the signal generator, 
and an observation of the frequency spectrum was selected. This is pro- 
duced by the command "TPLOTF" (teletype plot of the frequency function) 
which is followed by a request for the low and high frequency limits of 
the spectrum to be observed. These were specified to be 9.9 MHz to 10.1 
MHz as shown in Figure 13 and the program produced the spectrum plot shown 
in that figure. Note that the ordinate of the plot is automatically 
scaled. The spectrum of the AM signal is displayed with a carrier ampli- 
tude of 0.5 microvolts or -126 dB (one side of a two-sided spectrum); 
the amplitudes of each sideband component is 0.125 microvolt or -138 dB. 

A print of the frequency function is shown in Figure 14 which shows the 
amplitude and phase of these spectral components. 

The next block in the diagram of Figure 12b is that of a single sec- 
tion Butterworth bandpass filter with a total 3 dB bandwidth of 20 kHz. 
This block command along with the conmand to process to the output of 
block 2 is shown at the top of Figure 15. The 3 dB point of the filter 
was placed at the sideband frequencies of the AM signal so that the effect 
of the filter could be observed. The print of Figure 15 when compared 
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Figure 14. Listing of the Frequency Function at the Output 
of the Signal Generator. 


43 


* 


FVPNL'P, 1 . E7, 2. E4, l 
BWPN DP, I • E7, P. F4, 1 
PLO CK, P 
PLO CK, ? 

PROCESSING COMPLETE THRU PLOCK P 
PPINTE 
PPINTF 

ENTRP LOW, HIGH FFEPL'ENCIES 
9.9E5-6, 10. !E6 


LINE 

F PE P 

PEAL 

?01 5 

9.9000+Ob 

- . 000 

PO 1 6 

9 . 9100+06 

- . ooo 

?0 17 

g . gpoo+ofe 

- . 000 

PO 18 

g. 9 300+06 

- . ooo 

PO 19 

g. g400+06 

— . 00 0 

POP 0 

9 • g 500+ 06 

- . ooo 

P021 

g. g6oo+o6 

- . ooo 

POPP 

9.9700+06 

- . ooo 

POP3 

g.g« 00+06 

- . ooo 

P0P4 

9. 9900+06 

• ooo 

P0P5 

1 . 0000+07 

-. 000 

P0P6 

1. 0010+07 

- . ooo 

P0P7 

1. 0020+07 

- . 000 

pops 

|. 0030+07 

-.000 

POPP 

1. 0040+07 

- • OOO 

PO 30 

1. 0050+07 

- . OOO 

2031 

1. 0060+07 

- . 000 

P03P 

1 . 0070+07 

- • 000 

2033 

1. 00R0+07 

- . 000 

P03 4 

). 009 0+07 

- . 000 

2035 

1.0100+07 

- . 000 


IM AG 

MAG 

VB 

• 000 

• 000 

-P56. P? 

. 000 

• 000 

- 254. 06 

. 000 

. 000 

-PSP. 76 

a 000 

• 000 

-P50. PI 

a 000 

. 000 

- 248.41 

• 000 

• 000 

-?44. 51 

• 000 

• 000 

- p/iO. 53 

. 000 

. 000 

-P37. 58 

• 000 

• 000 

-PP8.95 

- • 000 

• 0 0 0 

- 141.07 

— • 000 

a 000 

- 1P6. OP 

— • 000 

a 000 

- 141. 07 

- . 000 

a 000 

-PP9. 48 

- . 000 

a 000 

-238.8? 

- . 000 

a 000 

-P4 1.69 

• a 000 

a 000 

- ?46. 1 0 

• • 000 

a 000 

- P 50. 78 

— • 000 

a 00 0 

-P5P.83 

— • 000 

a 0 00 

-P55. 64 

. 000 

a 000 

-P57. 49 

- . 000 

a 000 

-P59.P 1 


Figure 15. Listing of Frequency Function at the Output 
of the First Bandpass Filter. 
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with Figure 14 verifies this. The sidebands at 9.99 MHz and at 10.01 MHz 
have an amplitude of -141.07 dB out of the filter and had an ampli- 
tude of -138.06 at the input to the filter — a difference of 3.01 dB. 

Note that the carrier amplitude is unchanged. A plot of the spectrum at 
the output of the filter is shown in Figure 16. 

Next in the block diagram of the receiver is a mixer. This is mod- 
eled in FATCAT as an ideal mixer which produces sum and difference fre- 
quencies only. Another useful feature of FATCAT is the input formet re- 
quest command, shown at the top of Figure 17, which can be used when the 
input format of a block is not known. The request "INPUT FORMAT", fol- 
lowed by a conma, and then followed by the name of a block in the model 
library, yields the input parameters and their definition as shown m 
Figure 17 (use of this command has no effect on the circuit being pro- 
cessed). The ideal multiplier IDLMUL requires the peak amplitude of the 
local oscillator be entered as well as the frequency of the local oscilla- 
tor. This amplitude in the example was given as 1 volt and the frequency 
as 9 MHz. This should result in the spectrum of the AM signal being 
shifted to 1 MHz and 19 MHz. The spectrum around the difference fre- 
quency is shown in Figure 17. 

The next two blocks of the receiver are made up of a four section 
Butterworth bandpass filter with a center frequency of 1 MHz and a band- 
width of 20 kHz followed by an amplifier having a gain of 120 dB. The 
input commands for these blocks and the spectrum at the output of each 
block are shown in Figures 18 and 19, respectively. 

Proceeding along the block diagram of Figure 12b, the next block is 
an amplitude demodulator. The input conmands for this block and the out- 
put spectrum are shown in Figure 20. Note that the only two components 
in the spectrum displayed are the d-c component and the demodulated sig- 
nal at 10 kHz. 

Figure 21 shows the result of a call on the teletype time plot rou- 
tine. This routine was called to plot the time function at the output of 
the demodulator. The request for the entry of NSTART, NSTOP, and NJUMP 
allows the starting point in the array containing the time function to 
be specified, the stopping point to be specified, and the number of points 
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Figure 16. Spectrum at the Output of the First 
Bandpass Filter. 
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Figure 17. Spectrum at the Output of the Ideal Multiplier. 
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Figure 18. Spectrum at the Output of the Second Bandpass Filter. 
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Figure 19. Spectrum at the Output of the Amplifier. 
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Figure 20. Spectrum at the Output of the AM Demodulator. 
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Figure 21. Time Waveform at the Output of the AM Demodulator. 
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to be Skipped between plotted points to be specified. For the plot of 
Figure 21 the starting value was 1, the stopping value was 2048, and 64 

points were skipped between points plotted. 

Figure 22 illustrates the capability of printing the time function 
by use of the command "PR1NTT" . Following this command a request for 
the entry of the low and high array indices is printed. Figure 22 shows 
a print of the time function from array index 1 to index 30. 

A plot of the baseband spectrum at the output of the demodulator 
is shown in Figure 23 to demonstrate the flexibility of the plotting 
program. The spectrum can be seen to contain a positive and a negative 
frequency component as well as a d-c term. 

Input of the final block in the AM receiver is shown in Figure 24. 
This is a five section Tchebysheff lowpass filter with a corner frequency 
of 10 kHz and an inband ripple of 1 dB. Observation of the 10 kHz com- 
ponent in Figures 23 and 24 indicates that it has been reduced on the 
order of 1 dB as it should since this component lies at the filter cor- 

ner frequency. 

Exit from the program is effected by entering the command END OF 
JOB, which is followed in Figure 24 by the time required to simulate the 

receiver . 


2. FM Receiver 

The second example is that of the analysis of a proposed fre 
quency modulated receiver. The block diagram used for the FATCAT analy 

sis of the receiver is shown in Figure 25. 

The characteristics of the blocks making up the FM receiver are: 


( 1 ) 


A signal generator with a carrier frequency of 10 MHz, 
a modulation frequency of 10 kHz, no AM or PM modulation, 
FM modulation producing a peak frequency deviation of 
30 kHz, and a carrier peak amplitude of 1 microvolt. 


(2) A single section Butterworth bandpass filter with a band 
width of 80 kHz. 


(3) An ideal mixer with a local oscillator amplitude of 1 volt 
peak and a frequency of 9 MHz. 
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Figure 22. Listing of Time Function at the Output of the 
AM Demodulator. 
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Figure 23. Spectrum et the Output of the AM Demodulator . 
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Figure 24. Spectrum at the Output of the Low Pass Filter, 
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(4) A six-section intermediate frequency Butterworth band- 
pass filter with a center frequency of 1 MHz and a band- 
width of 80 kHz. 

(5) An amplifier with a voltage gain of 140 dB. 

(6) A frequency demodulator with a center frequency of 1 MHz. 

(7) A ten-section Butterworth low-pass filter with a cutoff 
frequency of 15 kHz. 

The frequency spectrum at the output of each block of the receiver, 
from block 1 (the output of the signal generator) to block 6 (the output 
of the frequency demodulator), is shown in Figures 26 through 31. Figure 
26 shows the output from the signal generator for a modulating frequency 
of 10 kHz and a peak frequency deviation of 30 kHz which gives a modula- 

tion index of 3. 

The output spectrum of the demodulator (Figure 31) shows distortion 
products produced by the bandpass filters in the receiver. Figure 32 is 
a time plot of the output of the frequency demodulator. A print of the 
frequency function at the frequency demodulator output is shown in Fig- 
ure 33. This print shows the amplitudes and phases of the distortion 
products over a greater frequency range than the plot because of the 
limited 50 dB range of the plot. 

Figure 34 shows the baseband spectrum at the output of a 10 section 
Butterworth low-pass filter with a cutoff frequency of 15 kHz. Figure 35 
shows a plot of the time waveform at the output of this filter, and Fig- 
ure 36 is a print of the frequency function at the filter output. It 
can be seen from Figure 36 that the ten-section Butterworth filter eli- 
minates virtually all distortion products except the 20 kHz component, 
and it reduces this component's amplitude approximately 60 dB. 

The cotmnand END OF JOB terminates the program and is followed by 
the time required to analyze the FM receiver. 

D. Conclusions 

A new circuit simulation program, FATCAT, for analyzing communica- 
tions receiver circuits has been developed to an operational state. 
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Figure 26. Start of FATCAT Run for FM Receiver, Showing 

Spectrum et the Output of the Signal Generator. 
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Figure 27. Spectrum at the Output of the Firet Bandpasa Filter. 
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Figure 28. Spectrum et the Output of the Ideel Multiplier. 
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Figure 29, Spectrum at the Output of the Second Bandpass Filter, 
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Figure 30. Spectrum at the Output of the Amplifier. 
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Figure 31. Spectrum at the Output of the FM Demodulator. 
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Figure 32. Time Waveform at the Output of the FM Demodulator. 
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Figure 33. Lieting of the Frequency Function at the Output 
of the FM Demodulator. 
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Figure 34. Spec t run et the Output of the Low Peee Filter. 
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Figure 35 . Time Waveform at the Output of the Low Paea Filter 
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Figure 36. Listing of tho Frequency Function et the Output 
of the Low Pete Filter. 
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although there are many ways in which the system can be improved. The 
simulation consists of block models of signal sources, filters, ampli- 
fiers, mixers, demodulators, etc., along with an operating framework. 

The program has been fully implemented on both Univac-1108 and SIGMA-5 
computers. Use of the present form of the program has shown it to be 
a useful tool for rapid analysis of communications circuits. 

The time and effort available for this work did not permit investi- 
gation of numerous ideas for improvement. Additional development work 
could produce a far more powerful and flexible program. Some of the 
areas in which additional work should be undertaken are. 

(1) The present signal generator block permits simulation of 
only simple modulation frequencies. Development of a modu- 
lator block in which an arbitrary modulating signal could 
be stored in a data array is needed. 

(2) Models for coupled circuits are needed. 

(3) Present form of the FM and PM demodulators operate satis- 
factorily only when the demodulator frequency is perfectly 
aligned with the carrier frequency, hence no d-c level is 
generated. Better demodulation models would be desirable. 

(4) Control software could be modified to permit insertion, 
deletion, or replacement of a given block. 

(5) A restart capability without the necessity of re-entering 
the circuit blocks would be helpful. 

(6) A routine to print out the circuit description is needed. 
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IV. TIME DOMAIN SIMULATION 


A. Introduction 

Digital simulation of a system governed by differential equations 
can be based upon either a time domain or a frequency domain description 
of the system. For time domain simulation, the describing differential 
equations are usually expressed in state variable form and the system re- 
sponse is obtained by numerical integration of the state equations. The 
time domain simulation method is applicable to both linear and nonlinear 
systems. For frequency domain simulation, the describing differential 
equations are expressed in transfer function form, and simulation is ac- 
complished by application of a fast Fourier transform (FFT) algorithm. 
Frequency domain simulation normally requires less computer time than 
does time domain simulation, but the transform theory upon which the fre- 
quency domain approach is based is limited to linear systems. According- 
ly, the simulation of linear systems is usually based upon the frequency 
domain while that of nonlinear systems is based upon the time domain. 

Cascade systems whose only nonlinearities are algebraic are excep- 
tions to the above classification and can be effectively simulated using 
the FFT algorithm. Since such systems have no feedback, the complete 
response (for the time range of interest) of the first block or subsystem 
can be obtained before considering the response of the following blocks. 
Similarly, the response of the second block can be found after the first 
block's response has been found. The output of a block represented by 
an algebraic nonlinearity (e.g., multiplication, saturation, polynomial 
nonlinearities, etc.) depends only upou the instantaneous input to the 
block and can be readily determined from the algebraic description of the 
nonlinearity without recourse to differential equations or frequency do- 
main transfer functions. A cascade system composed of linear dynamic 
blocks and nonlinear algebraic blocks can be simulated in a serial fash- 
ion using the FFT algorithm to transform back and forth between the time 
and frequency domains as required. This concept has been utilized at 
the Engineering Experiment Station in the development of the FATCAT simu- 
lation described in Section III of this report. 
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The digital simulation of communication systems with feedback and 
with nonlinear subsystems requires the use of the complete time domain 
formulation. The remainder of this section describes the development of 
a user-oriented communication system block diagram simulation program 
applicable to cascade or feedback, linear or nonlinear systems. Although 
the resulting time domain program has general applicability, its use is 
recommended only for those systems for which the more efficient FFT simu- 
lation program (FATCAT) is not applicable. 

A number of computer languages have been developed for digital simu- 
lation of dynamic systems in recent years. These programming languages, 
such as MIMIC, DSL- 90, CSMP, and CSSL [4], are intended to make it easy 
for an engineer experienced in analog computer simulation to use a dig- 
ital computer. The languages have simple input/output instructions and 
generally have efficient numerical integration algorithms. They can be 
considered block diagram simulators only in the limiting case where each 
block is no more complex than a single integration or addition; i.e., 
only when the block diagram is in reality an analog computer flow diagram. 
However, at least one of these languages, MIMIC, has a subprogram cap- 
ability which allows one to construct elementary blocks better suited to 
communication systems simulation than simple integrators and summers. 

With these block subprograms the communications system engineer can eas- 
ily prepare the necessary simulation input data directly from his block 
diagram. Using the MIMIC simulation language as a basis, a modular time 
domain simulation program, TIMSIM, has been developed for communication 
system block diagram simulation. In the following, the basic language 
MIMIC will be described briefly and its use with TIMSIM will be developed. 
Application of TIMSIM to the simulation of an AM receiver and an auto- 
matic gain control system will also be given. 

B. MIMIC: A Continuous System Simulation Language 

The simulation language MIMIC was developed at Wright-Patterson AFB 
in the mid- I960' s for the digital simulation of dynamic systems. De- 
tailed instructions for its use are given in the original MIMIC report [5] 
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and in the previously referenced textbook by Stephenson [4]. Reference 
should be made to these sources for a more complete description of MIMIC 
than appears in the following. 

The MIMIC language provides a set of functions (including integra- 
tion) specifically chosen to perform the operations necessary to solve 
systems of ordinary differential equations. A function is used by list- 
ing the name of the output (beginning in column 10) and the name and 
arguments of the function (beginning in column 19) on standard computer 
cards. For example, the equation x = log(y) would be programmed by 
punching "X" in column 10 and "LOG(Y)" in columns 19-24. A block-oriented 
program for simulating a dynamic system is obtained by first drawing a 
detailed block diagram of the system and then listing the interconnec- 
tions of the blocks. The block diagram and MIMIC program for a system 
described by x + x + x = 0 for zero initial conditions is given in Fig- 
ures 37 and 38. The first four lines in Figure 38 correspond directly 
to the connections at the four blocks in Figure 37. The name NEG2DX was 
arbitrarily selected to represent the negative of the second derivative 
of x. The first box follows from -X = x + x, and the second box and in- 
struction reverse the sign of the variable to give x. The third box and 
instruction correspond to the integration of x to give x while the fourth 
ones further integrate x to obtain the variable of interest, x. In the 
integration instructions, the "0” corresponds to the initial conditions. 
The detailed block diagram and resulting MIMIC instructions closely fol- 
low analog computer programming techniques, but do not require the ampli- 
tude scaling and time scaling required for analog computation. The last 
four MIMIC instructions in Figure 38 are bookkeeping instructions. The 
FIN(T,10) instruction causes the simulation to stop when time, T, reaches 
ten seconds. The HDR(TIME.X) instruction establishes the headings on the 
computer output as TIME and X. The 0UT(T,X) instruction means that the 
variables T and X are to be tabulated as simulation output. The last 
instruction, END, simply indicates the end of instructions. 

The CON and PAR instructions in MIMIC are used to load numerical 
values for constants and parameters (constants which change from run to 
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Figure 37. Detailed MIMIC Block Diagram for a System Described 
by x + x + x = 0. 


10 

19 

NEG2DX 

ADD(X.IDX) 

2DX 

NEG (NEG2DX) 

1DX 

INT(2DX,0.) 

X 

INT(1DX,0.) 
FIN (T, 10.) 
HDR(TIME,X) 
0UT(T,X) 

END 


Figure 38. MIMIC Instructions for Simulation 
of the System of Figure 37. 
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run). For example, if the system of Figure 37 had involved a constant, 
K2, in the differential equation, the instruction ”C0N(K2)", beginning 
in column 19, would precede the instructions in Figure 38 and the "END" 
instruction would be followed by a data card with the value of K2 entered 
in columns 1-12. 

As another MIMIC example, consider unity negative feedback around 
a lowpass filter as shown in Figure 39. The closed loop transfer func- 
tion is 


X _ K/(s + A) _ K ( 1 ) 

U “ 1 + K/ (s + A) s + (K + A) ’ 

and the associated differential equation is 

x = - (K + A) X + Ku . ( 2 > 

The corresponding MIMIC program and simulation results for u(t) = 1, 

K = 1, and A = 0.1 are given in Figure 40. 

An alternate method for MIMIC simulation of systems containing fre- 
quently occurring subsystems utilizes the subprogram feature of MIMIC. 
Rather than repeatedly deriving and programming descriptive differential 
equations for these subsystems, a library of subprograms can be estab- 
lished so that simulation of total systems is essentially reduced to 
describing the interconnection of the subsystems. Consider the develop- 
ment of a subprogram for the low pass filter previously considered. For 
subprogram purposes, rename the filter input and output variables as uj_ 
and xj, respectively. The describing differential equation can be 
shown to be iq = K-uj - A-x^, so that the describing MIMIC eqr rion is 

XI INT(K * U1 - A * XI, 0.) . (3) 

The associated subprogram is constructed by adding opening and closing 
instructions which name the subprogram and provide for arbitrary naming 
of the variables and constants within the subprogram. The resulting sub 
program (named FIRST) for the current example is as follows: 
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^x 

-v 


s + A 




Figure 39. Block Diagram of Feedback Example. 


@GT* LIB. MIMIC, IS TEST 
MIMIC-03.2-06/09/72-10: 32:17 


@ADL LARRY. 1 

1 * X 

2 * 

3* 

4* 

5* 


INT(-1.1*X+1.0,0.) 
FIN (T,5 . ) 
HDR(TIME,X) 
OUT(T,X) 

END 


**FURTHER DIAGNOSTICS AND EXECUTION FOLLOW*** 


TIME 

X 

0.00000 

0.00000 

1. 00000-01 

9.46962-' 

.20000 

.17953 

.30000 

.25552 

.40000 

.32360 

.50000 

.38459 

.60000 

.43923 

.70000 

.48817 

.80000 

.53202 

.90000 

.57129 

1.00000 

.60648 

1 . 1000 

.63800 

1 2000 

.66624 

1 . 3000 

.69154 

1.4000 

.71420 

1.5000 

.73450 

1.6000 

.75269 

1.7000 

.76898 

1.8000 

.78357 

1.9000 

.79665 

2.0000 

.80836 

2 . 1000 

.81885 

2 . 2000 

.82825 

2.3000 

.83667 

2.4000 

.84422 

2.5000 

.85097 


2.6000 

.85703 

2.7000 

.86245 

2.8000 

.86731 

2.9000 

.87166 

3.0000 

.87556 

3.1000 

.87905 

3.2000 

.88218 

3.3000 

.88499 

3.4000 

.88750 

3.5000 

.88975 

3.6000 

.89176 

3.7000 

.89357 

3.8000 

.89518 

3.9000 

.89663 

4.0000 

.89793 

4.1000 

.89909 

4.2000 

.90013 

4.3000 

.90107 

4.4000 

.90190 

4.5000 

.90265 

4.6000 

.90332 

4.7000 

.90392 

4.8000 

.90446 

4.9000 

.90494 

5.0000 

.90538 

5.1000 

.90576 


END 568 MLSEC 


Figure 40. MIMIC Simulation of System of Figure 39 
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10 19 

FIRST BSP(U1 , K, A) 

XI INT (K*U 1 - A*X1, 0.) 

FIRST ESP (XI) < 4 > 

The letters BSP and ESP are abbreviations lor "begin subprogram" and "end 
subprogram," respectively, and the arguments of BSP name the input vari- 
ables and constants while that of ESP names the output variable. In call- 
ing such a subprogram from the main MIMIC program, the variable names need 
not be the same as those used in the subprogram but must be in the same 
order in che calling statement. 

A MIMIC program using the low pass filter subprogram in the simula- 
tion of the system of Figure 39 is given in Figure 41. Note that lines 
2, 3, and 4 are a listing of the subprogram which has been obtained from 
a library as a package. Instruction 5 describes the feedback connection 
for u(t) = 1, while instructions 6 and 7 are the instructions required to 
"call" the desired subprogram. The letters CSP and ESP represent "call 
subprogram" and "end subprogram", respectively. Comparison of the argu- 
ments of CSP and BSP shows that the input to the low pass filter is 
called "E" in the main program and "Ul" within the subprogram, and that 
the constants "K" and "A" are to have the values 1.0 and 0.1, respectively. 
Note that the name of the subprogram being called is included on the CSP 
card. Instructions 7, 8, and 9 are not affected by the use of subpro- 
grams. Instruction 1 is a parameter statement which must be used tc ini- 
tialize the subprogram output variable for feedback system simulations, 
the initial value, 0., is entered in columns 1-12 of a data card. 

The use of library subprograms to represent subsystems allows the 
design engineer to digitally simulate his system without a detailed know- 
ledge of programing or of state variable techniques. In the following, 
a set of library subprograms for communication subsystems will be devel- 
oped and presented. The composite system is entitled TIMSIM, for time 
domain simulation. 
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@GT*LIB. MIMIC, IS TEST 
MIMIC-03 . 2-06/ 13/ 72-08 : 33 : 42 


C a ADD LARRY , 4 

1* 


PAR(X) 

2* 

FIRST 

BSP(U1,K,A) 

3* 

XI 

INT (K*U1-A*X1 , 0 . ) 

4* 

FIRST 

ESP (XI) 

5* 

E 

l.-X 

6* 

FIRST 

CSP(E, 1. , . 1) 

7* 


RSP(X) 

8* 


FIN(T,5.) 

A 


HDR(TIME,E,:.) 

10* 


0UT(T,E,X) 

11* 


END 


★♦FURTHER DIAGNOSTICS AND EXECUTION FOLLOW*** 

ENTER DATA NOW 
X 

0.00000 


Figure 41. Subprogram Example. 
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C. TIMSIM: A User-Oriented Communications System Block Diagram 

S imulation Program 

TIMSIM is not a new digital simulation language to compete with 
MIMIC, nor is it really a computer program. Rather, it is an operating 
philosphy and an expandable library of subprograms in the MIMIC language 
which represent the dynamics of communication subsystems. The communica- 
tion subsystems modeled for TIMSIM herein by no means exhaust the possi- 
bilities but do provide a representative group of subsystems. Included 
are subprograms which model AM signal generators, bandpass single-tuned 
filters, bandpass RC filters, mixers, and lowpass RC filters. Ihe mixer 
subprogram can also be used as a product detector. The method of develop 
ment of subprograms is described in sufficient detail to allow users to 
expand the subprogram library to meet their changing needs. 


1. Subprogram Derivation 

The single-tuned, bandpass filter subprogram STFIL is based 
upon the transfer function 



o 


as given in Pettit and McWhorter [6], where B is the filter 3 dB band- 
width (Hz), f Q is the filter center frequency (Hz), K is the center fre 
quency gain, and io 0 = 2rrf 0 . The transfer function can be rewritten as 


X sK 

U s 2 (f 0 /iu 0 B ) + s + (u) 0 f 0 /B) 


( 6 ) 


from which the corresponding differential equation could be obtained by 
inspection after cross multiplication with the interpretation of s r fac 
tors as n th derivatives. The resulting differential equation, however, 
would involve an undesirable derivative of the input, U. This undesir- 
able input signal differentiation can be avoided in the same manner as 
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with analog computer programming [7] by re-expressing the transfer func 
tion as 


X = i[2nB(KU - X) - £ (u> 0 2 X)] (7) 

where the complex variables appear only in nested multiplications by 1/s. 
The desired form of the describing differential aquation is obtained by 
drawing an elementary block diagram for this last expression, selecting 
the integrator outputs as the state variables, and writing the associated 
state equations (simultaneous first order differential equations). From 
Figure 42 the state equations are written as 

x = 2ttB (Ku - x) + y 



Figure 42. Single-Tuned Filter Flow Diagram. 
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Prior to writing the associated MIMIM subprogram, program names 
must be chosen for the variables. One of the requirements of MIMIC is 
that variable names not be repeated in other subprograms unless they 
represent exactly the same quantity. Accordingly, the TIMSIM subpro- 
grams are written with variable and constant names which end in a three 
digit number unique to that subprogram. Further, if a subprogram is re- 
peated in a simulation (e.g., two single-tuned filters in one receiver 
system), the third digit is increased by one in each additional copy of 
the subprogram. The variable names selected for the single-tuned fil- 
ter subprogram STFIL and their relationships to the algebraic variables 
are as follows: x - OUT110; y - Y110; u - XIN110; K - K110; B - B110; 

f Q - FRQ110. 

2. TIMSIM Subprogram Library . 

Table I contains the subprograms used to model AM signal gene- 
rators (AMSIG, AMSIG5) , single-tuned bandpass filters (STFIL), RC band- 
pass filters (BPFIL) , mixers (MIXER), and RC low-pass filters (LPFIL) . 
Included in the table are the definitions of the subprogram input, and 
output variables. 

3. Applications 

The operation of TIMSIM is best described by examples. The 
block diagram of Figure 43 describes a hypothetical communication sys- 
tem consisting of an ideal AM signal generator, single-tuned rf filter, 
mixer (frequency converter), single-tuned i-f filter, product detector, 
and low pass filter. The corresponding TIMSIM simulation program is 
given in Figure 44. The first five instructions establish the names 
and sequence of appearance of constants whose numerical values are to 
be entered on the data cards. Instructions 6 through 38 present the 
library subprograms which model the communication subsystem. Note that 
both the single- tuned filter and the mixer subsystems appear twice in 
the AM system diagram so that the corresponding subprograms must also 
appear in duplicate. As mentioned earlier, the repeated use of a 
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TABLE I 


TIMSIM SUBPROGRAM LIBRARY 


****BEGIN SUBPROGRAM FOR AMPLITUDE MODULATED SIGNAL GENERATOR.*** 
AMSIG BSP (XINIOO , FRQ 100 ,M100 , K100) 

THT .00 6 . 2831*FRQ100*T 

OUTIOO K100* ( 1 . -+M100*XIN100 ) *COS (THT100 ) 

AMSIG ESP (OUTIOO) 

****end of subprogram**** 

Definitions 

XINIOO = Modulating Voltage 
FRQ 100 = Carrier Frequency 
M100 = Modulation Index 

K100 = Gain 


****BEGIN subprogram for sin amplitude modulated signal generator. 

AMSIG5 BSP(XIN105,FRQ105,M105,K105) 

THT105 6 . 2831*FRQ105*T 

0UT105 K105*(l .4M105*XIN105)*SIN (THT103) 

AMSIG5 ESP(0UT105) 

****END OF SUBPROGRAM**** 


Definitions 

Variables defined similarly to AMSIG above. 


****BEGIN SUBPROGRAM FOR SINGLE-TUNED FILTER**** 

STFIL BSP(XIN110,B110,FRQ110,K110) 

OUT110 INT(YllO+(6. 2831*B110)*(K110*XIN110-OUT110) ,0.) 

Y110 INT(- (6 . 2831*FRQ110)* (6 . 283i*FRQ110 )*0UT110 ,0 . ) 

STFIL ESP(OUTUO) 

****END OF SUBPROGRAM**** 

Definitions 

XINllO = Input Voltage 
B110 = Filter Bandwidth, Hz 

FRQ110 = Filter Center Frequency, Hz 
K110 = Center Frequency Gain 

OUT110 = Output Voltage 


(Continued) 
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TABLE I (Continued) 


****BEGIN SUBPROGRAM FOR MIXER**** 

MIXER BSP(XIN120 ,FRQ120 ,KL20) 

OUT120 XIN120*K120*SIN (6 . 2831*FRQ120*T) 
MIXER ESP(OUT120) 

****END OF SUBPROGRAM**** 

Definitions 

XIN120 = Input Voltage 

FRQ120 = Mixer Product Frequency 

K120 = Gain 

OUT120 = Output Voltage 


****BEGIN SUBPROGRAM FOR LOW-PASS FILTER**** 

LPFIL BSP(XINI30,FCN130,K130) 

OUT130 INT(6 . 2831*FCN130* (XIN13O-OUT130) ,0.) 
LPFIL ESP(OUT130) 

****END OF SUBPROGRAM**** 


Definitions 

XIN130 = Input Voltage 
FCN130 = Corner Frequency, Hz 
K130 a DC Gain 
OUT 130 = Output Voltage 


****BEGIN SUBPROGRAM FOR RC BANDPASS FILTER**** 

BPFIL BSP (XINI40 , FL140 , FH140 , K140 ) 

WL140 6.2831*FL140 

WH140 6.2831*FH140 

Y140 INT(WL140*WH140*0UT140,0 . ) 

0UT140 INT(KI40*WH140*XIN140- (WL1404WH140)*(XJT140-Y140 ,0 . ) 

BPFIL ESP(OUTI40) 

****END OF SUBPROGRAM**** 


\ 


Definitions 

XIN140 = Input Voltage 
FL140 = Lower Corner Frequency, Hz 

FH140 = Higher Cornet Frequency, Hz 

K140 = Center Frequency Gain 

OUT140 = Output Voltage 



M-i 

Vi 



Figure 43. Hypothetical AM Communicat 









1* CON(FM,FC,M,K,FO,B) 

2* CON (FH1X ,KMIX) 

3* CON (F01 , BOl , KOI ) 

4* CON(FMIXl.KMIXl) 

5* CON (FLP.KLP) 

6*****BEGIN SUBPROGRAM FOR AMPLITUDE MODULATED SIGNAL GENERATOR.*** 

7 * AMSIG BSP(XIN100,FRQ100,M100,K100) 

8* THT100 6 . 2831*FRQ100*T 

9 * OOTIOO KlOO* ( 1 . 4M100*XIN100 )*COS (THT100 ) 

10* AMSIG ESP (OUTIOO) 

11*****END OF SUBPROGRAM AMSIG********************* 

12*****BEGIN SUBPROGRAM FOR SINGLE-TUNED FILTER******** 

13* STFIL BSP(XINL10 , BLIO, FRQ110,K110) 

14* OUT110 INT(Y110+(6.2831*B110)*(K110*XIN110-OUT1IO),0.) 

15 * Y110 INT(-(6.2831*FRQ110)*(6.2831*FRQ110)*OUT110,0.) 

16* STFIL ESP(OUTllO) 

17*****END OF SUBPROGRAM STFIL********* 

18*****BEGIN SUBPROGRAM FOR MIXER****************** 

19* MIXER BSP(XIN120,FRQ120,K120) 

20* 0UT120 XIN120*K120*ST.N (6 . 2831*FRQ120*T) 

21* MIXER ESP(OUT120) 

2 2***** END OF SUBPROGRAM MIXER********************* 

23*****BEGIN SUBPROGRAM FOR SINGLE-TUNED FILTER********* 

24* STFIL1 BSP(XIN111, Bill, FRQ111, Kill) 

25 * OUT111 INT(Y111+(6.2831*B111)*(K111*XIN111-OUT11i),0. ) 

26 * Ylll INT(- (6 . 2831*FRQ111)* (6 . 2831*FRQ111)*OUT111,0 . ) 

27* STFIL1 ESP(OUTlll) 

28*****END OF SUBPROGRAM STFIL1********* 

2 9*****BEGIN SUBPROGRAM FOR MIXER*********************** 

30* MIXER1 BSP(XIN121,FRQ121,K121) 

31* OUT121 XIN121*K121*SIN(6.2831*FRQ121*T) 

32* MIXER1 ESP(OUT121) 

33*****END OF SUBPROGRAM MIXER1** ************* **** ****** 

34*****BEGIN SUBPROGRAM FOR LOW-PASS FILTER********* 

35* LPFIL BSP(XIN130 ,FCN130,K130) 

36 * OUT130 INT(6.2831*FCN130* (XIN130-OUT130) ,0.) 

37* LPFIL ESP(0UT130) 

38*****END OF SUBPROGRAM LPFIL******************** 

39* X SIN (6 . 2831*FM*T) 

40* AMSIG CSP(X,FC,M,K) 

41* RSP(Y) 

42* STFIL CSP(Y,B,FO,l. ) 

43* RSP(Z) 

44* MIXER CSP (Z ,FMIX ,KMIX) 

45* RSP(ZZ) 

46* STFIL1 CSP (ZZ ,B01 , FOl , KOI ) 

47* RSP(ZZZ) 

48* MIXER1 CSP(ZZZ,FMIX1 ,KMIX1) 

49* RSP(OUTPUT) 

50* LPFIL CSP(OUTPUT, FLP.KLP) 

51* RSP(V) 


Figure 44. TIMSIM Lifting, Hypothetical Communication System. 
(Continued) 
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-*4 


A 


52* 


DTI 

53* 


DT2 

54* 


T1 

55* 


TF 

56* 


TTEST 

57* 


TLOGIC 

58* 


NTLOG 

59* 

TLOGIC 

DT 

60* 

NTLOG 

DT 

61* 



62* 



63* 



64* 




l./(50.*FM) 

1. /(20.*FC) 

(20 . / (4 .*FM) ) -DT2 

2. /FM 
T-Tl 

FSW(TTEST,TRUE .FALSE .FALSE) 
COM(TLOGIC) 

EQL(DTl) 

EQL(DT2) 

FIN(T.TF) 

HDR(TIME,X,V) 

OUT(T.X.V) 

END 


* -'FURTHER DIAGNOSTICS AND EXECUTION FOLIXM*** 


ENTER DATA NOW 

FM FC M 

10000. 1.00000406 .50000 


X FO 

1.00000-03 1.00000406 


ENTER DATA NOW 
FMIX 

5.45000405 


ENTER DATA NOW 
FOI 

4.55000405 


ENTER DATA NOW 
FMIX1 

4.55000405 


ENTER DATA NOW 

FLP 

10000 . 


KMIX 

10.000 


B01 KOI 

20000. 10.000 


KMIX l 
100.00 


KLF 

1.0000 


Figure 44. End. 


B 

20000. 
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subprogram within a TIMSIM simulation requires renumbering the variable 
names within each repetition. Accordingly, the second single-tuned 
filter subprogram is named STF1L1 and the internal variable names termi- 
nate in 111 rather than in 110. Instructions 39 through 51 constitute 
the heart of the TIMSIM input; these are the instructions which describe 
the system block diagram and specify the subsystem parameters through 
their calling of the subprograms. For instance, instructions 42 and 43 
specify that the subsystem is a single-tuned filter whose input is named 
y (the output of the AM signal generator) and that the filter bandwidth 
is the input quantity B, the center frequency is the input quantity FO, 
and the center frequency gain is unity. Instruction 43 states that the 
oucput of the filter is designated Z. Instructions 52-60 result in two 
separate time spacings between lines of printout; spacing DTI is used 
f or t < T1 while DT2 is used thereafter. This feature is used only when 
the simulation results of interest are preceded by a transient of no 
interest. Instructions 61-64 are as in the MIMIC example described 
earlier. Figures 45 and 46 present the simulation results in tabular 
and graphical form. 

The block diagram of Figure 47 describes a hypothetical Automatic 
Gain Control (AGC) system used to demonstrate the application of TIMSIM 
to feedback systems. The corresponding TIMSIM simulation program is 
given in Figure 48. The first three instructions describe the constants 
of the simulation; repeated simulations for new values of the constants 
in the PAR statements can be made without repeating the program listing. 
Following the library subprograms (instructions 4-28) , the interconnec- 
tion of the AGC system blocks is described by a series of subprogram call 
statements. Within this group, instruction 35 gives the effect of the 
feedback signal upon the forward gain, and instruction 44 demonstrates 
the use of the internal MIMIC function for a limiter. The remainder of 
the simulation is analogous to that of the AM system discussed earlier. 
The simulation results are tabulated in Figure 49. 
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TIME 

0.00000 

2.00000-06 

4.00000- 06 

6.00000- 06 
8.00000-06 

1.00000- 05 

1.20000- 05 

1.40000- 05 

1.60000- 05 

1.80000- 05 

2.00000- 05 

2.20000- 05 

2.40000- 05 

2.60000- 05 

2.80000- 05 

3.00000- 05 

3.20000- 05 

3.40000- 05 

3.60000- 05 

3.80000- 05 

4.00000- 05 

4.20000- 05 

4.40000- 05 

4.60000- 05 

4.80000- 05 

5.00000- 05 

5.20000- 05 

5.40000- 05 

5.60000- 05 

5.80000- 05 

6.00000- 05 

6.20000- 05 

6.40000- 05 

6.60000- 05 

6.80000- 05 

7.00000- 05 

7.20000- 05 

7.40000- 05 

7.60000- 05 

7.80000- 05 

8.00000- 05 

8.20000- 05 

8.40000- 05 

8.60000- 05 

8.80000- 05 

9.00000- 05 

9.20000- 05 

9.40000- 05 

9.60000- 05 

9.80000- 05 

1.00000- 04 


X 

0.00000 
.12533 
.24869 
.36812 
.48175 
.58778 
.68454 
.77051 
.84432 
.90482 
.95105 
.98228 
.99803 
.99803 
.98229 
.95106 
. 90484 
.84434 
.77053 
.68457 
.58781 
.48179 
.36816 
.24873 
.12537 
4.29471-05 
-.12529 
-.24864 
-.36808 
-.48171 
-.58774 
-.68451 
-.77048 
-.84430 
-.90480 
-.95104 
-.98228 
-.99802 
-.99803 
-.98230 
-.95108 
- . 90486 
-.84437 
-.77056 
-.68460 
-.58785 
-.48182 
-.36820 
-.24877 
-.12542 
-8.59539-05 


V 

0.00000 

-9.69228-04 

-6.16795-03 

-1.71452-02 

-3.61954-02 

-6.83026-02 

-.11594 

-.17394 

-.23469 

-.29958 

-.38015 

-.48357 

-.59961 

-.70867 

-.80481 

-.90508 

-1.0299 

-1.1756 

-1.3140 

-1.4232 

-1.5143 

-1.6194 

-1.7525 

-1.8900 

-1.9960 

-2.0634 

-2.1229 

-2.2064 

-2.3090 

-2.3936 

-2.4340 

-2.4449 

-2.4639 

-2.5082 

-2.5549 

-2.5694 

-2.5453 

-2.5109 

-2.4963 

-2.4999 

-2.4928 

-2.4537 

-2.3932 

-2.3409 

-2.3125 

-2.2943 

-2.2608 

-2.2047 

-2.1453 

-2.1069 

-2.0929 


Figure 45. TIMSIM Results, Hypothetical AM Communication System. 
(Contln ed) 
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1.02000-04 

1.04000-04 

1.06000-04 

1.08000-04 

1.10000-04 

1.12000-04 

1.14000-04 

1 . 16000-04 

1 . 16000-04 

1.20000-04 

1.22000-04 

1.24000-04 

1.26000-04 

1.28000-04 

1.30000-04 

1.32000-04 

1.34000-04 

1.36000-04 

1.38000-04 

1.40000-04 

1.42000-04 

1.44000-04 

1.46000-04 

1.48000-04 

1.50000-04 

1.52000-04 

1.54000-04 

1.56000-04 

1.58000-04 

j.. 60000-04 

1.62000-04 

1.64000-04 

1.66000-04 

1.68000-04 

1.70000-04 

1.72000-04 

1.74000-0': 

1.76000-04 

1.78000-04 

1.80000-04 

1.82000-04 

1.84000-04 

1.86000-04 

1.88000-04 

1.90000-04 

1.92000-04 

1.94000-04 

1.96000-04 

1.98000-04 

2.00000-04 

2.02000-04 

END 39256 


.12525 
.24860 
.36804 
.48167 
.58771 
.68448 
.77045 
.84427 
.90478 
.95102 
.98227 
.99802 
.99803 
.98231 
.95109 
. 90488 
.84439 
.77059 
.68463 
.58788 
.48186 
.36824 
.24881 
. 12546 
1.29795-04 
-.12520 
-.2/' "56 
-.36800 
-.48163 
-.58767 
-.68444 
-.77042 
-.84425 
-.90476 
-.95101 
-.98226 
-.99802 
-.99804 
-.98232 
-.95110 
-.90489 
-.84441 
-.77062 
-.68467 
-.58792 
-.48190 
-.36828 
-.24886 
-.12550 
-1.74173-04 
.12516 


-2.0830 

-2.0566 

-2.C173 

-1.9894 

-1.9916 

-2.0153 

-2.0339 

-2.0326 

-2.0267 

-2.0456 

-2.0982 

-2.1612 

-2.2036 

-2.2228 

-2.2486 

-2.3094 

-2.3969 

-2.4735 

-2.5141 

-2.5358 

-2.5782 

-2.6564 

-2.7415 

-2.7927 

-2.8040 

-2.8107 

-2.8478 

-2.9085 

-2.9521 

-2.9509 

-2.9216 

-2.9045 

-2.9169 

-2.9340 

-2.9193 

-2.8665 

-2.8058 

-2.7679 

-2.7506 

-2.7236 

-2.6652 

-2.5864 

-2.5178 

-2.4750 

-2.4436 

-2.3973 

-2.3291 

-2.2585 

-2.2105 

-2.1878 

-2.1696 


MLSEC 


Figure 45 . End . 
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Figure 47. Hypothetical Automatic Gain Control System. 



















C0N(FC,KP,BW,K1,VSET) 

PAR (X4 , FCR , KF , TPOINT , TF , KFORWD ) 

PAR (FM) 

SUBPROGRAM FOR LOW-PASS FILTER**** 

LPFIL BSP (XIN130,FCN130,K130) 

0UT130 INT(6.2831*FCN130*(XINI30-OUT130) ,0. ) 

LPFIL ESP(OUT130) 

F SUBPROGRAM*** 

SUBPROGRAM FOR SINGLE- TUNED FILTER**** 

STFIL BSP(XIN110,B110,FRQU0,K110) 

OUTIIO INT (Y110+(6 . 2831*B110)* (K110*XIN110-OUT110) ,0. ) 

YUO INT(- (6 . 2831*FRQ110)* (6 . 2831*FRQ110)*OUTI10 , 0 .) 

STFIL ESP (OUTIIO) 


? SUBPROGRAM**** 

SUBPROGRAM FOR SIN AMPLITUDE MODULATED SIGNAL GENERATOR. 


16* 

AMSIG5 

BSP(XIN105,FRQ105 ,M105 ,K105) 

17* 

THT105 

6 . 2831*FRQ105*T 

18* 

0UT105 

K105*(l .+M105*XIN105)*SIN(THT105) 

19* 

AMSIG5 

ESP(0UT105) 


20*****END OF SUBPROGRAM**** 

21 * 

i k 

***BEGIN SUBPROGRAM FOR RC BANDPASS FILTER**** 

22* BPFIL BSP (XIN140 , FL140 , FH140 , K140 ) 

23* WL140 6. 2831*FL140 

24* WH140 6.2831*FH140 

25 * Y140 INT(WL140*WH140*OUT140,0. ) 

26 * 0UT140 INT(K140*WH140*XIN140- (WL140+WH140)*OUT140-Y140,0. ) 

27* BPFIL ESP(0UT140) 


29* 

XIN 

SIN (6 . 2831*FM*T) 

30* 

AMSIG5 

CSP (XIN, FC,. 95,1.) 

31* 


RSP(XX) 

32* 

STFIL 

CSP(XX,BW,FC,1. ) 

33* 


RSP(X) 

34* 

XI 

KP*X 

35 * 

KVAR 

1. - K1*VF 

36* 

X2 

KVAR*X1*KF0RWD 

37 * 

X3 

X2*SIN (6 . 2831*FC*T) 

38* 

BPFIL 

CSP(X3,.5*FM, . 2*FC,1. 

39* 


RSP (X5) 

40* 

LPFIL 

CSP(X3,FCR, 1. ) 

41* 


RSP(X4) 

42* 

VF1 

KF*(X4-VSET) 

43* 

VFLIM 

• 9/KI 

44* 

VF 

LIM(VF1, -VFLIM, VFLIM) 


Figure 48. TIMSIM Listing, Automatic Gain Control System. 
(Continued) . 
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45* 

DTI 

46* 

DT2 

47* 

T1 

48* 

TTEST 

49* 

TLOGIC 

50* 

NTLOG 

51* TLOGIC 

DT 

52* NTLOG 

DT 

53* 


54* 


55* 


56* 



**FURTHER DIAGNOSTICS 


0.05/FCR 

,50/FC 

TPOINT - DT2 
T-Tl 

FSW (TTEST , TRUE , FALSE .FALSE) 
COM(TLOGIC) 

EQL(DTl) 

EQL(DT2) 

FIN(T.TF) 

HDR (TIME , X3 , X4 , VF , X5 ) 

0UT(T,X3,X4,VF ,X5) 

END 

D EXECUTION FOLLOW*** 


ENTER DATA NOW 

FC KP BW 

10000. 2.00000-03 5000.0 


K1 VSET 

1.00000-01 1.0000 


ENTER DATA NOW 

0. 100. 

X4 FCR 

0.00000 100.00 


3. 4.E-3 

KF TPOINT 

3.0000 4.00000-03 


6.E-3 1.E4 

TF KFORWD 

6.00000-03 10000. 


ENTER DATA NOW 

1000. 

FM 

1000.0 


Figure 48. End. 
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TIME 

X3 

X4 

VF 

0.00000 

0.00000 

0.00000 

-3.0000 

5.00000-04 

-4.07690-05 

2.8968 

5.6903 

1.00000-03 

1.02872-04 

2.6581 

4.9742 

1.50000-03 

-6.78247-05 

3.4256 

7.2767 

2.00000-03 

1.84120-04 

2.9041 

5.7124 

2.50000-03 

-9.51716-05 

3.4685 

7.4054 

3.00000-03 

2.84755-04 

2.9241 

5.7723 

3.50000-03 

-1.15230-04 

3.4720 

7.4159 

4.00000-03 

3.95770-04 

2.9258 

5.7773 

4.05000-03 

3.01825-04 

2.9397 

5.8191 

4.10000-03 

1.84063-04 

2.9858 

5.9573 

4.15000-03 

6.21273-05 

3.0570 

6.1709 

4.20000-03 

-4.57693-05 

3.1427 

6.4279 

4.25000-03 

-1.27702-04 

3.2317 

6.6947 

4.30000-03 

-1.79128-04 

3.3140 

6.9417 

4.35000-03 

-2.00990-04 

3.3825 

7.1474 

4.40000-03 

-1.97042-04 

3.4329 

7.2985 

4.45000-03 

-1.71074-04 

3.4630 

7.3889 

4.50000-03 

-1.25771-04 

3.4723 

7.4168 

4.55000-03 

-6.23912-05 

3.4610 

7.3829 

4.60000-03 

1.81439-05 

3.4298 

7.2894 

4.65000-03 

1.14925-04 

3.3801 

7. 1404 

4.70000-03 

2.24564-04 

3.3142 

6.9426 

4.75000-03 

3.39990-04 

3.2358 

6 » 7074 

4.80000-03 

4.49451-04 

3.1508 

6.4525 

4.85000-03 

5.37047-04 

3.0675 

6.2026 

4.90000-03 

5.85543-04 

2.9958 

5.9874 


4.95000-03 

5.81355-04 

2.9459 

5.8376 

5.00000-03 

5.20045-04 

2.9259 

5.7778 

5.05000-03 

4 09642-04 

2.9399 

5.8195 

5.10000-03 

2.68857-04 

2.9859 

5.9577 

5.15000-03 

1.21864-04 

3.0571 

6.1712 

5.20000-03 

-9.49593-06 

3.1428 

6.4282 

5.25000-03 

-1.10613-04 

3.2317 

6.6950 

5.30000-03 

-1.75682-04 

3.3140 

6.9419 

5.35000-03 

-2.05506-04 

3.3826 

7.1476 

5.40000-03 

-2.04383-04 

3.4329 

7.2987 

5.45000-03 

-1.76745-04 

3.4630 

7.3891 

5.50000-03 

-1.25781-04 

3.4723 

7.4170 

5.55000-03 

-5.30914-05 

3.4610 

7.3830 

5.60000-03 

4.04190-05 

3.4299 

7.2896 

5.65000-03 

1.53555-04 

3.3802 

7 . 1405 

5.70000-03 

2.82489-04 

3.3142 

6.9427 

5.75000-03 

4.19098-04 

3.2358 

6.7075 

5.80000-03 

5.49811-04 

3.1509 

6.4527 

5.85000-03 

6.56145-04 

3.0676 

6.2027 

5.90000-03 

7.17897-04 

2.9958 

5.9875 

5.95000-03 

7 . 18865-04 

2.9459 

5.8377 

6.00000-03 

6.53252-04 

2.9260 

5.7779 

6.05000-03 

5.29630-04 

2.9399 

5.8196 


Figure 49. 

T1MSIM Result*, 

AGC System 


X5 

0.00000 

-3.42534-02 

-1.5398 

-.38317 

-.63602 

-1.71602-02 

-.47952 

3.04278-02 

-.46301 

.34545 

1.1311 

1.7666 

2.1556 

2.2567 

2.0872 

1.7065 

1.1917 

.61588 

3.57353-02 

-.51181 

-1.0048 

-1.4290 

-1.7688 

-2.0009 

-2.0913 

-2.0006 

-1.6962 

-1.1710 
-.46135 
. 34709 
1.1327 
1.7682 
2.1572 
2.2584 
2.0888 
1.7080 
1.1931 
.61719 
7.69046-02 
-.51079 
-1.0039 
-1.4282 
-1.7683 
-2.0004 
-2.0910 
-2.0004 
-1.6960 
-1.1707 
-.46095 
.34764 


L 


D. 


Conclusions 


As stated earlier, TIMSIM is basically a simulation philosophy 
rather than a program or language; it has been presented here within 
the framework of the simulation language MIMIC. Several subprograms 
representative of communication system building blocks have been pre- 
sented and the method of generation of subprograms has been demon- 
strated so that additional subprograms can be generated as desired by 
the user. It is worth noting that the TIMSIM concept is also appli- 
cable to non-communication systems. Appropriate subprograms can be 
developed for mechanical systems, and control systems, etc. 
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APPENDIX A 


Listing of Software Modification to CIRCUS 
to Write Output Data into External File. 




(I 


* 

i 

! 


SUBROUTINE LIN<6A(PL- / 'T*NPNT*TIMF *DU V Y2 *LFL AG * .N'-'P » 

1 DELT * NR E AD ) 

INTEGER HOUR ( 2 > ♦ DAY ( 2 ) 

DIMENSION WORD ( 1 ) ♦ PLOT ( NPNT * 1 ) ♦ TT'Ell) 

DIMENSION NWORD(l) 

COMMON WORD* Ml *\5 »N6 »LPDS*LDS* JT I I L E l 1 * 2 ) 

COvviqn inoplt.ifloti ♦ MPLOTS *KPRNT *DUMY1 »DTC JR 

COMMON /SCRTCn/ T I TIE ( 1 vU) » ANAX ♦ Av I n »DA y *HOuR * I BAD* I MAX* 

I V, IN*KK *L A*NGO*NN* J»NWIPE*P I 1 * I 1 * I 2 * I 3 ♦ I 4* 

2 XPLCT 5 ♦ A 

F 0 J ! V A LFNCE ( WORD (12 6) * N7 ) * ( W 0 R D ( 3 8 ) ♦ N2 ) 

PQJTVAL'NCE ( WORD (33 ) * L I NTVL ) 

: G J I V A L r N C E ( WORD » MW ''PR ) 

'QU l VALENCE ( -'OROd :9) .LYVAX » 

'OJfVALENCE ( WCRO(ll-) *LYMIN ) 
fOUIVALEMCF < WORD ( 1 1 1 ) » NPCBLL ) 

► ~jj I VALENCE ( WORD (112) » LPLT ) 

,\R = 5 1 
\C = 21 

IF ( LFLAG-1 ) 39j*385*39J 


rrsE CONTINUE 
LrLAG = 2 

.vRITF (N7) ( PLO. (1*1)* I = 1 ♦ N '.»/ 1 D E ) 

N READ = NWIPB 
R-WIND N7 

READ ( N2 ) X PLOT S ♦ ( T I T LE ( I ) ♦ I = 1 . NPLO T S ) 

v>\ COiNTINUF 

V. RITE ( N 6 , 3 0 ) 

r***-* f; - a D TI“E AND PLOT VARIABLES FOR MPN I VALUES 
\ \i L - N a P 
NPL = N 3 CELL/6D 
DO 2i~ JX = 1*N?L 

% v n, = Na'L + ( 6 ( NPLOTS+1 ) ) - 1 

RF. A,, (,N2) (PLOT (1*1)* I =NWL *NWK ) 

N'\L = N' ,; < + l 
2K CONTINUE 

\PL = N D LOT S + l 

M •, L = NWP 

DO 230 JX = 1 » NPNT 

TIME ( J X ) = PLOT ( NWL ♦ 1 ) 

■ o 2 2 ) i =1 tN plots 
::<l = Nwl + 1 

PLOTtJX*!) = PLOT ( M VL » 1 ) 

^ 2 \j CONTI 4 0 r_ 

r ’ a L = NWL + 1 

? 2 :: CONTINUE 

CALL CLOCK ( 6HR5TIMC , IDUM) 

00 21 J= 1 « NPLOT S 
REF = PLOT ( 1 » J ) 

AMAX = REF 
AMIN * REF 


OF TIME. 
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LINK6A (Continued) 


C**** FIND MAXIMUM AND MINIMUM OF EACH VARIABLE 
DO 41 1*1* NPNT 

IF < AVAX.GT.PLOTI I .J ) ) GO TO 40 
AMAX = PLOT ( I » J ) 

IMAX = I 

AO IF ( AVIN.LT.PL0T ( I * J ) ) GO TO 41 
AMIN = PLOT ( I ♦ J ) 

I MIN = I 

41 CONTINUE 

IF { ABS ( AMAX-REF ) .OT • ABS ( A^I N-RE.F ) ) GC TO 42 
C***# USE *f I NDER* TO DETERMINE 1C - 9C PERCENT RISE AND Fall 
C**** T I vp S FOR PREDOMINATELY POSITIVE VARIABLE;?. 

PT1 = ,1*(Avin-REF)+REF 

call FINDER ( I'M N .NPNT * 1 » IMI N *T4.T 1 *PLCT ( 1 . J ) * T IME ) 
PT1 = .9*( AM IN -REF )+REF 

CALL TINDER ( I Ml N *NPNT * 1 » I M I N * T 3 * T2 »PLOT ( 1 ♦ J ) * T IME ) 

GO TC 43 

42 CONTINUE 

C**** USE *F I NDFR* TC DETERMINE 10 - 90 PERCENT RloE AND Fall. 
C**** TIMES EOR PREDOMINATELY POSITIVE VARIABLES. 

PT1 = . I*( AMAX-RtE )+REF 

CAcL FINDER ( 1 . I M AX ♦ I M AX ♦ NPNT » T 1 » T4 » PLOT ( 1 ♦ J ) * T I ME ) 

P T 1 = .9# ( AMAX-REF ) +REF 

CALL FINDER ( 1 » IMAX * I M AX .NPNT .T2 » 1 3 »PLOT( 1 » J ) *T IME ) 

43 CONTINUE 

C ** ■** WRITE OUT RISE AND FALL TIMES. 

C * * COMMENT BELOW SUPPRESSES PRINI OF RIoE AND FALL TIMES* 
C WRITE ( N6 * 5 C ) TITLE(J) *T1*T2*T3»T4 

C * * * 

21 CONTINUE 

<;**** IF ONLY RESPONSE TIMES DESIRED. RETURN. 

IF ( I PLOT I .F0.4 ) GO TO 80 
C 

C INSERT PLOTTING ROUTINES HERE. 

CALL PLTDTA ( PLOT ( N'VP . 1 ) »NR . NC » T I ME .PLOT .NPNT * NPLOTS . 

1 TITLE. DFLT.NPNT ) 

C 

call CLOCK <6H*PL0T*. IDUM) 

C * * NEXT TWO COMMENTS SUPPRESS LINE PRINTER PLOTS. 

C CALL PLCTER ( PLOT ( NWP . 1 ) » NR » NC ♦ T I ME .PLOT . NPN T .NPLwTS » 

C 1 TITLE.DELT.NPNT ) 

C 1 NPNT ) 

C * * * 

CALL CLOCK (6HR- p L0T , IDUM) 

80 CONTINUE 
RETURN 

30 F3RV!AT(////52X15HRESPCNSE TIVES//31V. ..7ARAVETER4* 

1 8H1C RISE 4X8H90 RISE4X0H1O FALL • X 3 H 0 FALL/) 

50 FORMAT I 33XA6. 1X1P4E12 .4 ) 

END 




.. 


‘JcROJT I NE PLTCTA ( NSCR ♦ NR * 'NC * T I V.E *PLOT *NPNT ♦NPLOTS » 
T I TLc.OFLT.NO I -n 

ivr NSK'I PLCT ( '‘O t m , vdloTS ) * T I TLE ( M^LCT 5 ) 

l MSC T ( MR * '-1C ) » TI^E(NPiMF) 

x I T F ( 1 9 > \P*T 

M T ( 1 '• ) ■ ‘PLOT c 

RITE (191 T I ' ‘ E 

P 1 T r ( 1 P ) PLOT 

PIT- (19) TIT L r 


APPENDIX B 


Listing of Software for Displaying Data Extracted from CIRCUS. 


NOTE: Many of the Hollerith strings in format statements were 

delimited with quote marks. The printer used to make 
the following listing did not have the quote character; 
a minus sign appears where each quote should have been. 


( 


> 


MAIN PROGRAM 


PARAMETER N = 32 

DIMENSION TIME! 5000 ) » PLATI5000) 
COMPLEX AIN) 

DIMENSION IBUFU000O) 

LOGICAL CPLT* PDPRT ♦ DOMFLG 
CPLT = .FALSE. 

PDPRT = .FALSE. 

DOMFLG = .TRUE. 

READ! 19) NPNT 
WR1TEI6.200) NPNT 
NP = N 

WRITEI6.202) NP 
READ ( 19 ) NPLOTS 
WRITEI6.201) NPLOTS 
CALL RDT FI NPNT NPLOTS. TIME. PL AT ) 
I GAM = A LOG ( N ) / AL0GI2.) + .1 
1 WRITEI6.101) 

WRITEI6.1001) 

WRITE I 6.23 Jl) 


READ! 5.100.ERR=l»nND=999 ) ISWTCH 
IF I ISWTCH .GT. 0) GO TO 2 
I F I ISWTCH .LT . o) GO TO 3 
CALL WRT DTA I NPNT. NPLOTS » T IME .PLAT ) 

GO TO 1 

2 IF! ISWTCH .GT. NPLOTS) GO TO 5C 
WRITEI6.102 ) 

READ 15.100) ISTART 
DO 10 I = l.N 

ITQ = I + ISTART - 1 + (ISWTCH - 1) * NPNT 
10 All) = CMPLXI PLAT! ITQ ) ♦ 0.) 

TP = I I T IMEI ISTART + N - 1) - T I ME I I START ) 
I F I • NOT • PDPRT) WRITEI6.203) TP 
PDPRT = .TRUE. 

CALL FFT I A. IGAM.-l ) 

CALL LFOLD(A.N) 

WRITEI6.103) 

RE AD I 5 . 1 DO ) JSWTCH 
IFIJSWTCH .LT. 100) GO TO 20 
CALL PRNTIA.N) 

JSWTCH = JSWTCH - 1J0 


) 


IN -1 ) ) * N 


I 

I 



102 



MAIN (Continued) 


f 20 IF ( JSW TCH .EQ. 0) GO TO 1 

DELF = 1. / TP 

i I F ( JSW TCH .LT. 10) GO TO 30 

I F ( • NOT • CPLT ) CALL PLOTS ( IBUF ( 1 )» 10000* 2 ) 
CPLT = .TRUE. 

WRITEC6.108 ) 

READ ( 5 » 1 00 ) FMAX 

CALL CPLOTF ( A , N .DELF , FMAX > 

JSWTCH = JSWTCH - lo 
IF(JSWTCH .EQ. 0) GO TO 1 
30 CONTINUE 

WRITE(6» 104) 

READ! 5.100) FLO 
wRITE(6. 105 ) 

READ ( 5 ♦ 1 00 ) FHI 

CALL TTFP (A*N»DELF»FLO»FrlI ) 

GO TO 1 

3 ISWTCH = AbS< IS wTCH) 

I F ( I SW TCH .GT. NPLOTS) GO TO 50 
WRITE(6*102 ) 

READ( 5*100) ISTART 
DO 11 I = 1 » N 

ITQ = I + ISTART - 1 + (ISWTCH - 1) * NPNT 
11 A ( I ) = CMPLX(PLAT ( ITQ) . 0.) 

TP = ( ( T I ME ( 1 ) - TIME(N)) / (N -1)) * N 
I F ( .NOT • PDPRT ) WRITEC6.203) TP 
PDPRT = .TRUE. 

§ wRITE(6*103) 

i READ ( 5 ♦ 100 ) JSWTCH 

| I F ( JSW TCH .LT. luu) GO TO 21 

I CALL PRNT(A.N) 

| JSWTCH = JSWTCH - 100 

| 21 I F ( JSWTCH .EQ. 0) GO TO 1 

f I F ( JSWTCH .LT. 10) GO TO 31 

| I F ( • NOT • CPLT) CALL PLOTS ( IbUF ( 1 ) » 10000 .2 > 

! CPLT = .TRUE. 

f C * * 

5 C PLACE CALCOMP TIME PLOTS HERE 

| WR I TE ( 6 »90u ) 

! C * * 





MAIN (Continued) 


JSWTCH = JSWTCH - 10 
I F ( JSWTCH • EQ • 0) GO TO 1 

31 CONTINUE 

WRITE (6*304) NP 
WR I T E ( 6 » 30 3 ) 

RE AO ( 5 » 1 00 ) NST ART 
NST = NSTART 
WRITE(6»3C5) 

READ (5*100) NST OP 
NSP = NSTOP 
wRITE<6»306) 

RE AD ( 5 » 1 00 ) N JUMP 

CALL TTTP ( A, NST .NSP.NJ'JMP) 

GO TO 1 

50 WRITE (6*107) NPlOTS 


GO TO 1 

100 FORMAT ( ) 

101 F ORM AT ( - ENTER 0 FOR PRINT OF DATA SET FROM CIRCUS-) 
1001 FORMAT (- ENTER PLUS DATA SET NUMBER FOR FREQ FCN-) 
2001 FORMAT (- ENTER MINUS DATA SET NUMBER FOR TIME FCN-) 

102 FORMAT (- ENTER I START- ) 

1 (j 3 FORMAT ( - ENTER 100 FOR PRINT* 010 FOR CALCOMP PLOT*-* 
A - OR U01 FOR TTY PLOT- ) 


104 

105 

107 

108 
200 
201 
202 
203 
303 
3u4 

/ 

305 

306 
900 
999 


FORMAT ( - 
FORMAT ( - 
FORMAT ( // 
FORMAT ( - 
FORMAT ( - 
FORMAT (- 
FORMAT (- 
FORMAT (- 
FORMAT (- 
FORMAT ( - 
4 

FORMAT (- 
FORMAT ( - 
FORMAT ( - 
I F ( CPL T ) 
STOP 
END 


ENTER FLO-) 

ENTER F H I — ) 

— ERROR--LARGEST DATA SET NUMbE k IS -»I2/) 

ENTER THE HIGhEST DESIRED FREQ IN THE SPECTRUM. FMAX-/ ) 
THE NUMbER OF POINTS PER DATA SET = -.15) 

THE NUMBER OF DATA SETS - -*I2/) 

THE NUMBER OF POINTS USED IN A TRANSFORM = -*I5) 

THE PERIOD OF THE TIME FUNCTION = -.1PE11.4.- SEC-/) 
ENTER NSTART-) 

ARRAY SI2E = -»I5»- NSTOP MUST BE EQUAL TO- 
• OR LESS THAN THIS VALUE-/) 

ENTER NSTOP-) 

ENTER NJUMP-) 

CALCOMP TIME PLOT NOT OPERATIONAL-/) 

CALL PLOT (v,. ,0., 999) 
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n n n ^ 


r * * 

r * \* tr v 

C 

c 

c 

r ■* * 


170 ') 


u : & ^ 


* 

* 


( 


RUr.RCJTINF C •’ L OT F ( A , ,\i , DE L F ♦ cv A X ) 

, T _ TH c FRLL r -v' I .\T ; CONTROL STATEMENT V-JST hkECEFD the 
r y ~ r u T F STATE v r ^ j pc* RUNS I'SlNfi CAlCO^P PLOTS. 


<3 u 5 =f UNIT #, TOPS 

C 0 P L r X A ( 1 ) 

I R = F AX / 0 F L ~ + .0-1 
A '■* A X = I R * SELF 
FLO = - FVAX 
r h I = FVAX 
FORVAT ( ) 

XTEST = A c: 0- ( F ul I -F L’J ) 

IF ( XTEST.LT.1.--3 ■) Ou ly 
XTEST = X T “ FT / ( A-iS ( "lu ) + Ai'S(Phl)) 

I F ( x r = 5 , T . L T . 1 . L - * ) r :0 TO 9>$9 

YSP'-EO = 7. . 

= \’/2 

NST = '12 - I R + 1 

NS® = N 2 + IR + 1 

T 1 = 1 . r> 3 b 

00 4,0 I = \ST .MS® 

T 2 = CA <S( A ( I ) ) 

IF( T2.GT.Ti ) T1 = T2 
03KAX = 20.*ALCG1 ( T 1 ) 

CALL SCALE ( O-'-'AX. ’AXSCL ) 

XVXSCL = -AXSCL 
X < S C A L = 1 '.**( -x ’XSCL/2C . ) 


call 

FACTOR ( 

. 4 ) 



CALL 

pl: 

TIC. 

i-2~ 

.,3) 


call 

PLO 

T( 1? 

• * 0 • 

,-3) 


call 

PLC 

T ( — 1 

< • * • 

1 4 , , 3 

) 

DC 3 

I = 

1.2 




CALL 

PLO 

T f — 1 

• % < 

.,2) 


CALL 

pl: 

T ( 1 0 

• * ^ • 

,2 ) 


CALL 

□ lo 

T ( 1 0 

. ,-14. ,2 ) 


call 

PLO 

T C - 1 

.0 1 

,-14, 

♦ 2 ) 

call 

PLC 

T ( - 1 

.0 1 


.2 ) 

call 

®LC 

t ( l : 

. 1 . 

C . v 1 , 

2 ) 

call 

PLC 

T(l~ 

.'1, 

-14.0 

1 .2 ) 

CALL 

®LC 

‘ T { — 1 

. ♦ ~ 

14. 1 

, 7 ) 

00 3 


I a gap: = 

1»? 


DC 

LOOP 

* i-j 

cal I 

3RATE 

$ THF 

DO 1 

. I 

■ • . 

7 C 



Y * 

-14, 

► + V- 

.2*1 



I F ( 

'•*00 1 

11,5) 

.HO 

:. ) 

GC TO 


4 



CPLOTF (Continued) 


CALL PLOT (-1 - .1 »V *3 ) 

GO TO 3 

6 IF (MODI I* 1C) .FO. 0) GC T L 7 
CALL PLCT(-1C.16»Y*3> 

GO TC 3 

7 CALL PLOT ( - 1 • 2 » Y » 3 ) 

5 CALL PLOTI-l • ♦ Y .2 ) 

l; CONTINUE 

Y = 2. 

DO 15 I - 1 ♦ 3 
J = I - 1 

Y = Y-2. 

YY = X''X$CL - 1 ' • *J* YSPk FD/7C • 

15 CALL N , J : ’ / 3Ek(-11.0 9*Y-.1^5*.21*YY»D«.-1) 

CALL SYN':!0L(-11.2A,-9.4,.21*14hAVPLIT'JDE ( 03 ) .90 .0 ♦ 14 ) 

FCE^TR = (FLO + FHI1/2. 
f-jPP-R = F H I - F C c. V T P 
I- (FUPPFR.GT.l.-l GO TO 21^0 
1LXP = ALOGIC ( FUPPcR ) ~ 1 
GO TO 2 1C 1 
2 1-' CONTINJE 

X I F X P = ALCG 1C ( F JPPER ) 

I F X P = XIEXP 
PlfyD = IPXP 

IF( ( ASS(XIEXP-RIEXD) .LT.l .E-2L ) . AND. (XIFXP.GE.9TFXP) ) 

1 IPXP = IFXP - 1 
2101 CONTINUE 

FULFCL = FUPPEP*< 1 ' .**<-IEXP) > 

IFLSCL = F’JLSCL 

I T E w P = 10.*FULSCL 

RITE M P = I T E V P 

SCALE 1 = IFLSCL 

SCALE1 = 1-.*SCALE1/FUL£CL 

TtMFS = 1G.MFLSCL 

SCALE1 = SCALE1/TENIFS 

c 00 LOOP 40 CALIBRATES THE 30TT0 M POSITIVE EORDEP 

C * * 

00 4 J I = 0 . I TEMP 
X = i*scalei 

IF (^00(1*5) .FQ. C) GO TC 32 
CALL PL0T(X»-14.1»3) 

GC TC 36 

32 IF (voD(I.lC) .FQ« 0) GC TO 34 
CALL PLOT ( X t-14. 16 » 3 ) 

GO TO 36 

34 CALL PLOT ( X »-14.2» 3 ) 

36 CALL PLOT ( X t-14. .2 ) 

4 CONTINUE 


u u 




CPLOTF (Continued) 


c * * 

C DC LOCP !4T CALIBRATES THE BOTTOM NEGATIVE BORDER 

C * * 

DO 140 I = 1*1 TEMP 
X = - 1 *SC A L F 1 

I F ( '»*0D ( I » 5 ) • EG • ^ ) GO TO 132 
CALL PLOT( X *-14.1 *3) 

CO TO 136 

132 IF ( V00( I »U ) .FO.O ) GO TO 134 
CALL PLOT ( X *-14.16 *3 ) 

GO TO 136 

134 CALL PLOT ( X ,-14.2 * 3 ) 

136 CALL PLOT ( X *-l 4. *2 ) 

140 CONTINUE 

DO L00 D 240 CALIBRATES THE TOP POSITIVE BORDER 

C * * 

00 240 1 = 0,1 TEMP 
X = i*scalei 

1 F ( MOD ( I » 5 ) • EQ • 0 ) GO TO 23 2 
CALL PLOT ( X , • 1 ,3 ) 

GO TO 236 

232 IF(V»OD( I *1C ) .FC. ) GO TO 234 
CALL PLOT ( X , . 1 6 ,3 ) 

GO TO 236 

234 call plot ( x * . 2 * 3 > 

236 CALL PLOT ( X *0 • *2 ) 

240 CONTINUE 

Q * * 

C DO LOOP 34G CALIBRATES THE TCP NEGATIVE BORDER 

C * * 

DO 340 I = 1*1 TEMP 
X = -I*SCALE1 

I F ( VOO ( I *5 ) .FO.O ) GO TO 332 
CALL PLOT (X,.l,3) 

GO TO 336 

332 I E ( MOD (I»10),FG*0) GO TC 334 
CALL PLOT (X,*16,3) 

GO TO 336 

334 CALL PLCT(X,.2»3) 

336 CALL PL0T(X*G.*2) 

34 J CONTINUE 

C - - 

C DC LOCP 20 CALIBRATES THE RIGHT HAND BORDER 

C » * 

00 20 I * 0,70 
Y = -14. + 0.2*1 
IF <*00 (1,5) .EQ. 2) GO TO 16 
call PLOT ( K.1*Y,3) 

GO TO 18 
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CPLOTF (Continued) 


16 IF (VOD(ItlC) .cQ. J) GO TO 17 
CALL PLOT ( 1 - • 1 6 » Y » 3 ) 


17 

13 

20 


10 • 2 » Y ♦ 3 ) 

1C. .Y.2 ) 


* 0 • ♦ 


-n 


GO TO 18 
CALL PLOT ( 

CALL PLOT( 

CONTINUE 

A K 1 = 1C./FULSCL 
CALL NU v ftFR(- 0 .06,-14. 5..21, 

DO 200 I = l.IFLSCL 
X PO A = -.06 + I *A< 1 
XNEG = -.12 - I *A<1 

CALL NUMBER (XP0S»“14.3, .21»1.*I.0.»-1) 

CALL NUMBER (XNEo* -14 .5* . 2 1 »-l • * I * 0 . ♦ -1 ) 

CALL SYMBOL (-A.5»-14*9» .21 * „ 

1 46 H ( FREQUENCY - FCENTFR ) DIVIDED BY F/dCALF, 

CALL SYMBOL ( -4. 5 ♦ -1 5 . 3 ♦ • 2 1 ♦ 1 OH FC ENTER - ♦ -J • * 1 ' ) 

CALL NU M c F R ( - 2 • 0 ♦ - 1 5 • 3 * • 2 1 * FC E N T R » 0 . ♦ 0 ) 

CALL SYMBOL! j.,-15.3,. 21 .9HFSCALE = »0..9) 

f ALL N U M B F. R ( 2 . 4 ♦ - 1 5 . 3 » . 2 1 . 1 0 . * * I E X P ♦ ' ♦ 0 ) 


C3000 CONTINUE 

I STOP = NAP - NAT 1 

q£MCM = ^sP - N C T 

DELX = 2C./DFN0M 

XS = - ( N2 + 1 - NST ) * DELX 

DO 3 0 I = 1 » I STOP 

I I = I - 1 

T 1 = (CABS(A( I+NST-1 ) ) )*X<SCAL 

I F ( T 1 .LT. l.F-7) GO TO 3« 

y i = (7J./YSPRFD) * 4. * ALCG10CTI) 

IF(Yl .LE. —14. 0 ) GO TO 30 

XI = II * DELX * XS 

I F ( X 1 .LT. -l r .) GO TO 30 

I F ( X 1 .GT. 1'.) GO TO 3 0 

IF(Vl.L r .O*I GO TO 2 *> 

CALL SYMhOL(Xl».16».71»l H **D.»l) 

GO TO 30 
230 CONTINUE 

CALL PL0T(X1»-14..3> 

CALL PL0T(X1.Y1 .2 ) 

CALL PLOT ( X 1 1 - 1 4 • * 2 ) 


3j CONTINUE 

CALL PLOT ( 13. . -20. .-3 ) 

36 FORMAT ( 1H1 *2X « 1 4MPL0 T COMPETED) 

/.'RITE ( 6 * 33 ) 

9099 CONTINUE 

return 


END 


I 


. ,46 1 


\ 
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SUBROUTINE FFT ( A, IGAM • I SN ) 

COMPLEX A ( 1 ) *T1*T2»T EMP 

DOUBLE PRECISION P 1 2 • SO » CO »S I » C I »SN » CS 

p 1 2 = 6*28318 53071795864800 

N = 2 ** IGAM 

NBIT * 36 - IGAM 

N1 = N - 2 

DO 30 I = 1*N1 

IFLIP * 0 

IX * 1 

DO 10 J * 1*IGAM 
I OLD = IX 
IX = IX / 2 
IBIT = IOLD - 2 * IX 
10 IFLIP = 2 * IFLIP + IoIT 
IF (I *LE. IFLIP) GO TC 30 
11=1+1 
12 = IFLIP + 1 
TEMP = A ( I 2 ) 

A( 12 ) = A ( I 1 ) 

A ( I 1 ) = TEMP 
30 CONTINUE 

DO 80 I = 1 ♦ I GAM 
MEL * 2** I 
NEL2 * NEL / 2 
NSET = N / NEL 
SI = DS I N ( P I 2 /NEL ) 

Cl = DCOS ( P I 2/NEL ) 

DO 80 J = 1 » NSET 

I NCR * ( J - 1 ) * NEL 

SO = O.ODO 

CO = l.uDO 

DO 80 II » 1.NEL2 

J1 = II ♦ I NCR 

J2 * J1 ♦ NEL2 

T 1 = A ( J 1 ) 

T2 * A ( J2 ) * CMPLX ( CO y ISN * SO) 

A ( J 1 ) * T 1 + T 2 
A ( J2 ) * T1 - T2 
SN * SO * Cl ♦ CO * SI 

CS * CO * Cl - SO * SI 

CO = CS 
80 SO = SN 

IF (ISN *GT. 0) GO TO 12o 
DO 110 I « 1 » N 
110 A ( I ) » A( I )/N 
120 CONTINUE 
RETURN 
END 


SUBROUTINE LFOLD (A*N) 
COMPLEX A ( 1 ) *TI 
N2=N/2 

DO 10 1*1* N2 
I I = I +N2 
T 1 = A ( I ) 

A ( I )*A (II) 

A ( II ) =T1 
Kt T URN 
END 


n n 




lu 
1 JO 
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SUBROUTINE PRNT(A.N) 

COMPLEX A ( 1 ) 

WRI T E ( 6 » 10 1 < 

DO lu I = 1*N 
Db = 1.E30 
T = CABS ( A ( I ) ) 

IF(T • GT • 0*) Do = 20. * ALOGlO(T) 

WRITE(6»100) I »A ( I ) »Db 

FORMAT (1X.I5.1P2E1 5. 4,5X,0PF8.2 ) 

FORMAT (//»- LINE-, 7X, -REAL-,1 IX ,- I MAG- , 1 2X , -Db-/ ) 

RETURN 

END 


SUBROUTINE RDTF ( NPNT ,NPLOTS, TIME, PLOT ) 

DIMENSION TIME! NPNT ) , PLOT ( NPNT , NPLOTS ) 

READ ( 19) TIME 

READ! 19) PLOT 

RETURN 

END 


l 


SUBROUTINE SCALE ( DBMAX »MAXSCL ) 

* * * THIS SUBROUTINE ESTAoLISHES ORUlNATt SCALING FOR 
THE REMOTE SPECTRUM PLOTTER. 

IF (DBMAX. LE.U. ) GO TO lu 
MAXSCL - 0 

1 MAXSCL = MAXSCL + lo 
DIFF * DBMAX - MAXSCL 
IFIDIFF.GT.O. ) GO TO 1 
GO TO 999 

10 MAXSCL » 0 

11 MAXSCL = MAXSCL - 10 
DIFF * i DBMAX - MAXSCL) 

IFID1FF.LE.U. ) GO TO 11 
MAXSCL * MAXSCL + lo 

999 RETURN 
END 




A 
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SUBROUTINE TTFP ( A , N . DELF »FLO» FH I ) 

C * * THIS SUBROUTINE PROVIDES A TELYTYPE PLOT OF 
C SPECTRUM FROM FLO TO FHI 

COMPLEX A ( 1 ) 

DIMENSION I A ( 50 ) »MM ( 6 ) 

C * * * * 

NST = ( N/2 ) + INTIFLO/DELF + SIGN!. 5. FLO)) 

NST = NST +1 

N5P s. (N/2) + lNT(FHl/DtLF + SlGN(.5»cHl)) 

NSP = NSP + 1 

C * * * * 

DBM AX = -I.E30 

I END = NSP “ NST + 1 

DO 1 I = 1 * I ENO 

DtCTMP = CABS (A (NST + I - 1)) 

IF(DECTMP.LT.1.E-3U) DECT MR = l.E-30 

b = 20. * ALOGlOt DECTMP ) 

I F ( b.GT.DBMAX) DbMAX = B 

1 CONTINUE 
C * * * * 

* * * * 

WR I TE ( 6 » 2 ) I END 

2 F ORMAT ( /5X *8HNS IZt = »I5/) 

C * * * * 

CALL SCALE ( DbMAX »MAXSCL ) 


THE FREQUENCY 


C 

c 

c 


WILL VARY FROM ( MAXSCL- 50 ) DB UP 


* * * 

THE ORDINATE 
TO MAXSCL DB. 

DO 33 I = 1*50 
33 I A ( I ) * 1H 
DO 5 I = 1*6 

5 MM I > * MAXSCL - lo*(6-I) 
MAXF = AoS(FLO) 

I F ( ABS ( Fhl I.GT.MAXF) 

NAMEF = 0 

IF (MAXF.GT.1.E3) NAMEF = 3 
I F ( MAXF .GT • 1 • E6 ) NAMEF = 6 
IF(MAXF.GT.1.E9> NAMEF = 9 
IF (NAMEF.EQ.O ) WRITE (6»2uv^ 
I F ( NAMEF .EQ. 3 ) WRITE(6*203) 

I F ( NAMEF .EQ.6 ) WRITE(6*20b) 
IF ( NAMEF. EQ.V ) WRITE(6*2b9) 


MAXF = A3S(FHl ) 


I 





% 

§ 
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TTFP (Continued) 


200 FORMAT (2X »l4hFREQUENCY (HZ) »9X»8HDtClBELS) 
203 FORMAT (2X.15HFREQUENCY ( XHZ > * 8X ♦ 8HDEC I bELS ) 
206 FORMAT (2X.15HFREQUENCY (MHZ) »8X * 8HDEC I BELS ) 
209 FORMAT ( 2X . 1 5HFREQUENCY ( GHZ ) * 8X . 8HDEC I dELS ) 
WR I TE ( 6 » 7 ) ( MM ( I ) ♦ I = 1.6) 

7 FORMAT (/7X»I4»4(6X*I4) ♦ 5X » 1 4 ) 

WRI TE( 6 .3 ) 

6 FORMAT (9X.1HI ,5(10H + 1)) 

F F ACT = 1. 

I F ( NAMEF • EQ • 3 ) FFACT = l.E-3 

I F ( NAMEF • tQ .6 ) FFACT = l.E-6 

IF (NAMEF. EQ. 9 ) FFACT = l.E-9 

FLO = FLO * FFACT 

FHI = FHI * FFACT 

0ELF1 = DELF * FFACT 

FLOPRT = DELF 1* ( NST -1 -M/2) 

DO 10 I = 1 » I END 

OECTMP = CABS ( A ( NST + I - D) 

I F ( DECTMP* LT • 1 .t— 3u ) DECTMP = l.E-30 
B = 20. * AL0G1 0 ( DECTMP ) 

M = 50 + B - MAXSCL 
J = I - 1 
XJ = J 

FREQ * FLOPRT + XU * DELF1 
I F ( M » L T .U> GO TO 50 
IF(M.EQ.O) GO TO 34 
DO 1515 II = 1 »M 
1515 I A( I I ) = 1H- 

WRITE(6»35) FREQ.IA 
DO 1616 II = 1 »M 
1616 I A C I I ) = 1H 

35 F ORMAT ( F8 • 3 » 2H I.50A1) 

GO TO 36 

34 WR1TE16.37) FREO 
37 FORMAT ( IX .F7.3.2H -) 

GO TO 36 

5u *RITE(6.51> FREQ 
51 FORMAT ( 1X.F7.3.2H I ) 

1000 FORMAT ( ) 

36 CONTINUE 

10 CONTINUE 
WRITEI6.il ) 

11 FORMAT (6X .4HFREQ) 

RETURN 

END 


SUBROUTINE TTTPIA.NST *NSP*NJUMP) 

COMPLEX A ( 1 ) 

0 IMENSION I A ( 50) 
dMAX = 0* 

d M I N = BMAX 

DO 1 I = NST.NSP.NJUMP 
a = R E AL ( A ( I ) ) 

IFtb.LT.dM IN) BMIM = d 
i I F ( d.GT .UMAX ) BMAX = U 

IF ( ( BMAX-0M I N ) .LT • 1 • E-3 0 ) GO TO 999 
OC 33 I - 1*50 
33 I A ( I ) = 1H 
100 wR I T E ( 6 » A ) 

4 FORMAT ( 5X » 1M+ ) 

WRITE (6 *3) BMlN.oMAX _ 

3 FORMAT ( 12X»-AV.PLITUDE- MIN -*E?.A,-» MAX -*E9. 
WR I T t ( 6 * A ) 

WRITE (6 .6) , , . , w _ u . 

6 FORMAT ( 8 X * 2 FH U » 3X » 2FU 1 *3X » 2H • 2 * 3 X * 2 Ft • 3 * 3 X * 2H • A 

1 3X*2H.6*3X»2H.7»3X»2rl.8»3X»2H.9»2X*3Hl.Q) 

WRITE (6*7) 

7 FORMAT ( 1H*6X»1HI»10(5H I)) 

DO 10 I = NST » NSP » N JUMP 

B = REAL ( A( I ) ) 
b = (B-BMIN)/ (BMAX-BMIN ) 

M = INT(0*5G.+C.5> 

IF(M.EQ.O) GO TO 3 A 
I A : M ) = 1H* 
w R ’ T E ( 6 * 3 5 ) I * I A 
I A ( M ) = 1H 

35 FORMAT (2X*I5»3H I.5oAI) 

GO TO 36 

3 A WRITE(6*37) I 
3 7 FORMAT ( 2 X * 15 »3H *) 

36 CONTINUE 

10 CONTINUE 
WRlTEl6tll ) 

11 FORMAT (6X*1HN) 

GO TO 900 

999 WRITE (6*9) 

9 FORMAT ( IX ♦ 12HERR0K FINISH) 

900 RETURN 
END 


A,- VOLTS-) 
* 3X * 2H. 5 » 




SUBROUTINE WRTDTA(NPNT *NPLOTS*TIME*PLOT > 
DIMENSION T I ME ( NPNT ) ♦ PLOT ( NPNT »NPLOT S > 
WR1TE(6»110) 

READ (5*102) ISW 

I c ( ISW .GT. 0 ) GO TO 1U < 

WK ITE ( 6* 101 ) 

WR l TE ( 6 ♦ 1 1 1 ) 

WR I TE ( 6 * 100 ) TIME 


WR I TE ( 6 * 100 ) TIME 
GO TO 999 

IF(ISW .GT. NPLOTS) GO TO 20 
WRITE(6»101 ) 

WR 1 TE ( 6 ♦ 1 12 ) ISW 

wRITE(6»100) (PLOT ( I » ISW) * I = 1 »NPNT ) 

GO TO 999 

WRITE(6»113) NPLOTS 
FORMAT ( 1X.1P6E11.A ) 

FORMAT ( /// ) 

FORMAT ( ) 

FORMAT(- ENTER u FOR TIME LISTING OR DATA 
FORMAT (3UX.-TIME-/ ) 

FORMAT (25X.-DATA SET NUMdER -*I2/) 

FORMAT ( //- ERROR — LARGEST DATA SET NUMotR 
WR I TE ( 6 ♦ 102 ) 

RETURN 

END 


NUMBER- ) 


-♦ 12 / ) 
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APPENDIX C 


FATCAT PROGRAM DESCRIPTION AND LISTING 
1. General Description 

The program is coded in FORTRAN IV and the Univac-1108 version con- 
sists of a main program and 31 subroutines} in addition one of the sub- 
routines (CPLOTF) which generates plots of frequency spectra on a CAL- 
COMP plotter requires calls to 5 other subroutines contained in a plotter 
control package. This plotting routine and associated plotter control 
subroutines are not used in the SIGMA-5 version. 

Several of the subroutines contain multiple entry points; the total 
number of subroutine and function entry names in the Univac-1108 version 
is 47. These names are listed in alphabetical order in Table C 1 and 
those which are not subroutine names are identified. 

In the following sections, the main program and all subroutines 
are briefly described, and each description is followed by a listing of 
the routine as used on the Univac-1108. For listings of the SIGMA-5 
versions of those routines that were modified for that machine, see 
Appendix D. 


NOTE: Many of the Hollerith strings in format statements were delimited 

with quote marks. The printer used to make the following listing 
did not have the quote character; a minus sign appears where each 
quote should have been. 


PHb CEDD!0 
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TABLE C 1 

ALPHABETICAL LISTING OF ALL SUBROUTINE ENTRIES IN FATCAT 
(All names denote subroutines unless marked otherwise) 


1 . 

ADJN 

25. 

LFOLD 


2. 

AMDEMO 

26. 

LIM 


3. 

AMP 

27. 

LSTCGM 


4. 

BCDFPT 

28. 

NUMBER -1 " 


5. 

BWBNDP 

29. 

PDCHK 


6. 

BWBSTP* ** (BWBNDP) 

30. 

PERIOD 


7. 

BWHIP* (BWBNDP) 

31. 

PHDEMO 


8. 

BWLCJWP* (BWBNDP) 

32. 

PLOT 4 " 


9. 

CHBNDP 

33. 

PLOTS + 


10. 

CHBSTP* (CHBNDP) 

34. 

PROCES 


11. 

CHHIP* (CHBNDP) 

35. 

PRTFAC* 

(PERIOD) 

12. 

CHLCWP* (CHBNDP 

36. 

SCALE 


13. 

CPLOTF 

37. 

SIGGEN 


14. 

ELFIND 

38. 

STRDTA 


15. 

FACTOR + 

39. 

SYMBOL 4- 


16. 

FETCH 

40. 

SYNBP 


17. 

FFT 

41. 

SYNHP* 

(SYNRP) 

18. 

FILTER 

42. 

SYNLP* 

(SYNBP) 

19. 

FLATSP 

43. 

TELPLT 


20. 

FMDEMO 

44. 

TIMFCN 


21. 

FRQFCN* (TIMFCN) 

45. 

TTFP 


22. 

FRQMUL 

46. 

WRFF (WRTF) 

23. 

IDLMUL 

47. 

WRTF 


24. 

INPFOR 





* Eh try Point In subroutine named in parenthesis. 

**Not used in SIGMA-5 version. 

+CALCOMP plotter subroutines called by CPLOT. 
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2 . MAIN PROGRAM 

Calls: PLOTS*, FETCH, ELFIND, WRTF, WRFF, TELPLT, TTFP, CPLOTF, PRTFAC , 

STRDTA, PROCES, LSTCOM, INPFOR, PDCHK. 

Commons: blank, CFREQ, CDOM, CD AT A, CCIRKT, CWORD, CFLGS. 

Description : MAIN is the overall controlling program which directs the 

operations of command and data input, interpretation of input, data 
storage, and command execution. Most of the detailed work in all opera- 
tions is carried out by subroutines. 

Program Listing: 


PARAMETER NMAX = 2048 
COMPLEX A(NMAX) 

COMMON iM » I GAM* DEL F » JELT .PD.CARRF 
COMMON /CFREG/ NFk*Fk(6) 

COMMON /CDOM/ DOMFLG 
COMMON /CDATA/ jCT * *DAT A(2«u) 
COMMON /CCIRKT/ NdLK . I T YP ( 3o » 2 > 
COMMON/CWOKO/ WuRD(lu) 
CCMMON/CFLGS/ PDFLG*ARFLG 
LOGIC PDF LG* ARFLG 
C DIMENSION I dUF ( 5000 ) 

C CALL PLOTS! IbUF ( 1 ) »5uou*2 ) 

JCTR = 1 

I TYP ( 1*2 ) = JCTR 
NFR = 0 
NoLK. = 0 
IdLK. = 0 
PDFLG = .FALSE. 


ARFLG = 

•FALSE. 




WRITE (6 

* 7006 ) 




DO 2 I 

= 1*10 




WORD! I ) 

= 6H 




CALL FET CHI WORD* 

L* NdAD) 


IF (NOAD 

.EQ. 0) < 

30 TO 

1 


CALL ELFINDIWORD 

* LTYP) 


GO TO! 

10* 20* 

30* 

20* 

5 G » 

1 

>— 
»— * 
o 

* 

ro 

o 

13^* 

140* 

1 5 o ♦ 

2 

210* 220* 

2 3u. 

24 o » 

25o. 

3 

310* 320* 

33v* 

34v> » 

35v* 

4 

410* 42u* 

43J) 

.LTYP 



60* 70* 80* 90* 10u» 

16v» 17v>* 18u» 19«j» 2 0 o * 

26v> * 27o* 28 0 » 29u » 3vu* j 

36o » 3 7o » 38o* 39 j ♦ 


★CALCOMP plotter routine. 
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MAIN (Continued) 


10 IF ( WORD ( 3 ) .EQ. 1H ) GO TO 12 

N 1 = WORD ( 2 ) 

N2 = WORD ( 3 ) 

GO TO 15 
12 WR I T E ( 6 * 70^4 ) 

READ ( 5 »700o ) N1»N2 
15 CALL WRT F ( A » Nl *N2 ) 

GO TO 1 

2u IF ( WORD ( 3 ) .EU. 1H ) GO TO 22 
FRLO = WORD ( 2 ) 

FRHI = WORD ( 3 ) 

GO TO 25 

22 WRI TE ( 6 »70U5 ) 

RE AD ( 5 ♦ 7 000 ) FRLO.FRH I 
25 IF (LTYP .EG!. 4) GO TO 40 
CALL WRFF( A.FKLO.FRrtI ) 

GO TO 1 

C » * TTY TIME PLOT 

3 u IF ( WORD ( 3 > .EG. 1H ) GO TO 32 
NbT = WORD ( 2 ) 

MSP = WORD ( 3 ) 

N JUMP = WORD ( 4 ) 

IF ( N JUMP .EO. 1H ) N JUMP = 1 
GO TO 35 

32 WR I T E ( 6 » 7007 ) 

READ ( 5 *7000 ) NST *NSP *NJUM.° 

35 IF ( N JUMP • LT • 1) N JUMP * 1 
CALL TELPLT ( A*NST »NSP*NJUMP ) 

GO TO 1 

c * * TTY FREOUENCY PLOT 

40 CALL TTFP ( A »FRLOtFRHl ) 

GO TO 1 

C * * CALCOMP TIME PLOT 

50 WRITE(6*7101 ) I < I TYP ( I * J I . J * 1 * 2 )* 
7101 FORMAT (1X*2( I 3 *3X ) ) 

GO TO 650 

C * * CALCOMP FREQUENCY PLOT 
60 CALL CPLOTF(A) 

GO TO 1 

C * * PRINT PRIME FACTORS 
70 CALL PRTFAC 
uO TO 1 

c * * END OF JOB 
60 GO TO 999 


1 * 5 ) 


I 

j 
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MAIN (Continued) 


C * * BUTTERWORTH BANDPASS 
90 CALL STRDTA ( 3 »0 *0 ) 

NTYP = 3 
GO TO 600 

C * * BUTTERWORTH LOWPASS 
100 CALL STRDT A ( 2 * 0 »0 ) 

NTYP = A 
GO TO 600 

C * * BUTTERWORTH HIGHPASS 
110 CALL ST KOTA (2*0*0) 

NTYP = 6 
GO TO 600 

C * * BUTTERWORTH BANDSTOP 
120 CALL STRDTA! 3,0,0) 

NTYP = 6 
GO TO 600 

C * * CHEBYSHEV BANDPASS 
130 CALL STRDTA ( A » O ,0 ) 

NTYP = 7 
GO TO 6o0 

C * * CHEBYShEV LOWPASS 
1 A 0 CALL STRDTA! 3,0,0 ) 

NTYP = 8 
GO TO 600 

C * * CHEBYSHEV HIGHPASS 
150 CALL STRDTA! 3,0,0) 

NTYP = 9 
GO TO 600 

C * * CHEBYSHtV BANDSTOP 
160 CALL STRDT A ( A » 0 *0 ) 

NTYP * 10 
GO TO 600 

C * * SYNCHRONOUS BANDPASS FILTER 
170 CALL STRDT A ( 3 , 0 ,0 ) 

NTYP = 11 
GO TO 600 

C * * SYNCHRONOUS LOWPASS FILTER 
1B0 CALL STRDT A ( 2 , 0 ,0 ) 

NTYP = 12 
GO TO 600 

C * * SYNCHRONOUS HIGHPASS FILTER 
190 CALL STRDT A ( 2 , 0 ,0 ) 

NTYP = 13 
GO TO 600 
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MAIN (Continued) 


C * * SIGNAL GENLRATOR 
200 CALL STRDTA (6*2*1) 

CARRFQ = WORD ( 2 ) 

NT YP = 1 
GO TO 600 

C * * FREQUENCY MULTIPLIER 
210 NTYP = A 
GO TO 600 
220 GO TO A30 
C * * IDEAL MULTIPLIER 
230 CALL STRDTA(2,1»2) 

NTYP = 16 
GO TO 6u0 
2 AO GO TO A3 0 
C * * FM DEMODULATOR 

250 CALL STRDTA ( 1*0*0) 

NTYP = 19 
GO TO 600 

C * * PHASE DEMODULATOR 
260 CALL STRDTA( 1*0*0) 

NTYP = 20 
GO TO 6u0 
C * * AMPLIFIER 

270 CALL STRDTAC 1*0*0) 

NTYP = 2 
GO TO 600 
C * * LIMITER 

280 CALL STRDT A ( 3 » 0 *0 ) 

NTYP = 17 
GO TO 600 

290 NOUT * WORD ( 2 ) 

IFINOUT *LE • N6LO GO TO 291 
WR I T E ( 6 * 700 1 ) NBLK. 

GO TO 1 

291 I F ( NOUT - IBLK ) 295*295*293 
293 IF ( PDFLG) CALL PDCHK. 

ITMP * IBLA ♦ 1 

DO 292 IdLK * ITMP* NOUT 

IBTYP = ITYP( IBL»C*1 ) 

JCTR * I TYP ( I BL< * 2 ) 

CALL PROCES(IBTYP.A) 

292 CONTINUE 
IBLK. * NOUT 

295 WR I TE ( 6 *7002 ) IbLK. 

GO TO 1 




MAIN (Continued) 


300 CONTINUE 
310 CONTINUE 
320 GO TO 370 
330 GO TO 430 
C * * AM DEMODULATOR 
340 CALL STRDT A ( 1 *0 *0 ) 

NTYP = 22 
GO TO 600 

C * * FLAT SPECTRUM GENERATOR 
350 CALL ST RDT A ( 3 ♦ 0 *0 ) 

CARRFQ = 0. 

NTYP = 23 
GO TO 600 
36u CALL LSTCOM 
GO TO 1 
370 GO TO 650 
380 CALL INPFOR 
GO TO 1 
390 GO TO 650 
400 GO TO 650 
410 GO TO 650 
420 GO TO 650 
430 wRITE(6.7C03) 

GO TO 1 

60U NbLK * NBLK + 1 

I TYP ( NbLK* 1 ) * NTYP 
ITYP!NdLK + 1.2) = JCTR 
GO TO 1 

650 WRITE!6.71v»0) WORD(l) 

GO TO 1 
K-00 FORMAT ( ) 

7001 FORMAT!- * * ERROR * * LARGEST BLOCK NO IS -.12*- * *- > 

7002 FORMAT!- PROCESSING COMPLETE THRU BLOCK -.12) 

7003 FORMAT!- * * UNDEFINED STATEMENT * *-) 

7004 FORMAT!- ENTER LOW. HIGH INDICES-) 

7005 FORMAT!- ENTER LOW. HIGH FREQUENCIES-) 

7006 FORMAT!- START-) 

7007 FORMAT!- ENTER NET ART » NSTOP. N JUMP- ) 

7100 FORMAT!- COMMAND -»A6»- IS NOT YET OPERATIONAL-) 

999 CONTINUE 

C CALL PLOTSCO. .0..999 ) 

STOP 

END 


f 


* 
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3. 


Subroutine ADJN 


Called by: PDCHK, PERIOD, FLATSP 

Calls : none 

Commons : blank 

Entries : none 

Description: ADJN adjusts N (number of data samples) to be a power of 2 

If the current value of N (either that computed to meet the Nyquist cri- 
terion or that entered by the user) is already a power of 2 it is not 
changed; otherwise it is adjusted upward to the next power of 2. 

Program Listing: 


SUBROUTINE ADJN 

COMMON N * I6AK ♦ DELF *DELT »PD 

IGAM = A LOG ( h ) /A LOG ( 2 • ) + »999 

N = 2**1 GAM 

DELT = PD / N 

RETURN 

END 


4. 


Subroutine AMDEMO 


Called by: PROCES 

Calls: FRQFCN, TIMFCN 

Commons : blank, CDATA 

Entries : none 

Description: AMDEMO simulates an ideal amplitude demodulator. Operating 

on the frequency spectrum, the negative frequency components are all set 
to zero to give a spectrum characteristic of a complex time function. 

The positive frequency components are then shifted down in the data array 
by an amount corresponding to the center frequency of the demodulator, 
thus positioning the spectrum at baseband. Transforming to the time 
domain produces a complex time wave form; conversion to a real time wave 
form is effected by replacing each time sample with one whose real part 
is the absolute value of the complex sample, and whose imaginary part is 
set to zero. 

Program Listing: 


SUBROUTINE AMDEMO ( A ) 

COMMON N* IGAM»DELF ♦DtLT »PD»CAKKFQ 
COMMCN/COAT A/ JCT K »DATA ( 2<J<j ) 

COMPLEX A ( 1 ) 

PI2 = 6.2831853 

N2 = N /2 

FO = WORD(JCTR) 

CALL FRQFCN ( A ) 

C * * REMOVE THE NEGATIVE FREQUENCY CompuimEimTS 
DO 10 I = 1.N2 

10 A ( I ) * (0. ♦ u. ) 

C * * MOVE THE MODULATED CARS I EH Tu ZE*u FkEuUEimCY 
I FO = Fu / DELF + .5 
NST ART = N2 ♦ 1 
NSTOP1 s N ~ I F v + 1 
DO 11 I = NSTART *N 

11 A ( I - I F U ) = All) 

DO 12 I = NSTOP 1 *N 

12 A ( I ) = ( 0. . 0 • ) 

C * * RECOVER THE AMPLITUDE INFORMATION 
CALL TIMFCN(A) 

DO 20 I * 1 » N 
TEMP = CABS ( A ( I ) ) 

2u A ( I ) * CMPLXITEMP. o.) 

RETURN 

END 
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5. 


Subroutine AMP 


Called by: PROCES 

Calls: none 

Commons: blank, CDATA 

Entries : none 

Description: AMP simulates an amplifier. Its action is simply to mul 

tiply each data sample by a constant. Since the multiplication is the 
same in both time and frequency domains, AMP accepts the data array in 
either domain. 

Program Listing: 


SUBROUTINE AMP ( A ) 

COMPLEX A ( 1 ) 

COMMON N » IGAM 

COMMON/ CDATA/ JCT* , DAT A ( ) 
XP = DAT A ( JCTR ) / 2o. 

G = 10.**XP 
DO 100 I = 1*N 
100 A ( I ) = G*A ( I ) 

RETURN 

END 


6 . Function BCDFPT 


Called by: FETCH 

Calls : none 

Commons : none 

Entries : none 

Description: BCDFPT accepts binary coded characters representing numeri- 

cal quantities and converts them to a real number which is returned 
through the function name, ftiis routine was adapted from CIRCUS with 
a few minor changes . 

Program Listing: 


n ICr 
I = 
J = 
< = 


o T rc'V-T^ -AT A FaO m TC FLOATING POI'T. 

A \ AA.RA Y C -NT A I ,v I *:G T*'" *' -CO 
r A3i TO CONVERTED* 

= l\:-X :F T-if CHAR ACT 'p -EING CONV^TL-. 

= i \ A L X C R A Z S P r 0 I N G T 0 T H F DIGIT J - 1 • 

= 1 VIFN ,,ECD r : I i':G NJ V 3FP PORTION. 

? uhfn decoding fractional pcrti-n. 

3 WmEN CODING cX p O vr NT • 

DIGIT* r » n LU 3 » DECPT » rCH 
EXPPLC-. DI-n-LC. D r C c L 0 » FXSIGN 

N '• ( . r ' ( 1 ) * < F I C *1 ( 1 ) * I N T ~ 0 ° ( ~ ) * P F S' ; L T ( 3 ) * P I G I T ( 1 3 ) 


INTEGER DIGIT* r . D LUS* DFCPT « rU- 
LOGICAL FxP^LC* PI” LG* D ’ C r L G » ::>G! iN 
it >/-n; I 'N L-(. r '( ! ' *<SI GM I ? ) « INTFOP ( 3) ♦RFSMLT < 3 ) *1 
/, j A I A I T / 1 H ♦ HI * 1 H? , 1^° » 1 * 1 H r - * I HA ♦ 1 M7 » ] MB , 1 uc / 

DATA PLJS. MIVUF. E » DFCPT / H+, 1H-. HE* 1 M| 

FXP-LG = • F a l F r . • 

J I o F L G = .FALSI. 

D c C " L G - • - ■■ L S c • 

FXSIGV - .FALSE. 

DO 11 <=1*3 

< S I G \ ( < ) = 1 
INTEGRA ) = 

CONTINUE 
NPLART = ~ 

< = 1 

D 0 31 I - 1 » f : 

I CM A? = : COM I ) 


FUNCTION BCDFPT (Continued) 


TEST FCR SION* DIGIT ♦ DECIMAL POINT* OR E 

IF ( I CHAR-PLUS ) 13*22*13 
li IF ( ICHAR-VINUS ) 14,24*14 

14 DC 15 J = 1 * 1 0 

IF ( I CHAR-0 I G I T ( J ) ) 15*25*15 

15 CONTINUE 

I- ( ICHAP-DECPT ) 16*26*16 

16 IF ( ICH4R-F ) 21*29*21 


2 3 I F ( DIC'FLG ) 
GO TC 31 


TO 2 8 


24 IF ( DIGFLG ) GO TO 27 
< 5 I GN ( 1 ) = -1 
GO TO 31 


PLUS SIGN 


MINUS SIGN 


DIGIT FRO'* C TO Q 


2 5 I \TEC-R ( < ) = 1 ' J * I N T £ G R ( < l+J-l 
nplart = NPLAPT+.<-l 
digflg = .true. 

GO TO 31 

DEC I v AL POINT 

ONLY ONE DECIMAL POINT PER NU V L' F” IS ALLOWED. 
DECIMAL POINT IS NOT ALLOWED IN EXPONENT. 

26 IF ( DECFLG ) GO TO 21 
IF ( EXPFLG ) GO TO 21 
DECFLG - .TRUr. 

< = 2 
GO TO 31 


FUNCTION BCDFPT (Continued) 


*E* FOR EXPONENT 

SLAN< TIME'S TFN ** EXPONENT NOT ALLOWED. 

<SIGN(3) = -1 

IF ; FXSIC-N ) 00 TO 21 

EXSIGN = .TRUE. 

UQ TO 3 :■ 

IF ( EXPFLG ) GO TO 21 

IF ( .NOT. DIGFLG ) GO TO 21 

EXPFLG = .TRUE . 

< = 3 

N D L A S V = NPLAPT 
CONTI VJE 


TUP NUMBER HAS BEEN SEPARATED INTO INTEGER. FRACTION. 
AND EXPONENT PARTS. COVE I NE THE’" TO FOR'" T ,J E 
number in: floating POINT. 

IF ( EXPFLG ) GO TO 32 
EXPCN = 1. 

GO TO 35 

CALCULAT- EXPONENT. an EXPONENT v Av BE ONLY T’-.’O 
DIGITS LONG AND LESS THAN 38 IN MAGNITUDE. 

IF ( NPL ART-NPL ASV-4 ) 33,33*21 
IEXPON = I NT FOR ( 3 ) *<S I GN ! 3 ) 

IF ( I A 5 S ( I C XPCN ) - 37 ) 3A, 34,21 
EXPCN = 1 0 •** I EXPCN 
NPL AR T = NPLASV 

CALCULATE MANTISSA 

RTSHFT = 1 . .**NPL ART 

RESULT ( 1 ) = FLOAT! I NTEGR ( 1 >*KS I GNll) ) 

RE SULT ( 2 ) = FLOAT! I NTEGR ( 2 ) *KS I GN ! 1 ) ) / RTSHFT 

BCDFPT = ( RESULT! H ♦RESULT 12) )*EXPON 

RETURN 


ILL-GAL CHARACTER OR BAD SYNTAX. 


N = -1 
GO TO A 1 
END 
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7. Subroutine BWBNDP 


Called by: PROCES 

Calls: FILTER 

Commons : CDATA, CFILT 

Entries: BWBSTP, BWLOWP, BWHID 

Description: BWBNDP simulates a Butterworth bandpass filter; auxiliary 

entries prod-.ce simulations of Butterworth band stop, high pass, and low 
pass filters. 

The Butterworth filter produces a maxiramally flat response defined 
by the function 


l H VI • -r=R 

\A + (U) p ) 

where n is the filter order and uu p is a normalized frequency, 
of Equation (l) are given by 


The poles 


S = e j6 
b k 6 


where 

e . ( Vl±sl -1) (f) , v. - l. 2. ... 2". «> 

With the poles known, the transfer function can be written in terms of 
the n poles lying in the left half-plane as 


H(S) * (S - S 1 )(S - s 2 )(s - s 3 ) ... (S' s n ) 

The Butterworth filter models are implemented with Equations (2), (3), 
and (4). Ihe computation of the transfer function. Equation (4), is 
carried out in subroutine FILTER. Subroutine BWBNDP computes the values 
of S k , and sets up variables from which the normalized frequency, o p , 
can be determined to satisfy the definitions 

m * uo/ui for low pass, 

P c 


for band pass. 
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* 


UJ 


p 


a) /ot> for high pass, 
c 


Bo) 

2 

cjj = for band stop, 

p cu - u> 
r o 


where 


(jo * 3 dB corner frequency, and 
c 

Bou * full 3 dB bandwidth. 

All Butterworth filters calculated with this model will exhibit 3 dB 
attenuation at the corner frequency (low pass and high pass), or at one- 
half the bandwidth away from the center frequency (band pass and band 
stop). The computed transfer function of the latter two are symmetrical. 

Program Listing: 


SUBROUTINE BwBNDP(A) 

COMPLEX A( 1 ) *S ( 20 ) 

COMMGN/CDATA/ jCT R .DATA ( 2u« ) 
COMMON/CFILT/ Fu.FCOFF.NR. amp.fflg.s 
LOGICAL FFLG 
FFLG = .TRUE* 

GO TO 2 

ENTRY bWOSTP(A) 

FFLG = .FALSE. 

2 F0 = OATA(JCTR) 

FCOFF = DAT A ( JCTR + 1 ) / 2. 

NR = DAT A ( JCTR+2 ) 

GO TO 6 

ENTRY BWLOrtP(A) 

| FFLG = .TRUE. 

| GO TO A 

i ENTRY BWHIP(A) 

f FFLG = .FALSE. 

I 4 F0 = 0. 

! FCOFF * DAT A ( JCTR ) 

| NR = DAT A ( JCTR+ 1 ) 

V 6 AMP = 1. 

t DO 10 < * 1*NR 

ThETA * 1.5707963 * ( ( 2 • ♦ NR-D/NR) 

10 S(K) * CMPLX (CCSCTHETA) » S I N I TrlETA ) ) 
CALL FILTERIA) 

RETURN 

END 
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8. Subroutine CHBNDP 

Called by: PROCES 

Calls : FILTER 

Commons : CDATA, CFILT 

Entries: CHBSTP, CHLCWP, CHHIP 

Description: CHBNDP simulates a Tchebysheff bandpass filter; auxiliary 

entries produce simulations of Tchebysheff bands top, high pass, and low 
pass filters. 

The implementation of the Tchebysheff (equal ripple) filter model 
is identical to that used for Butterworth filters except for the compu- 
tation of the poles. The poles for the Tchebysheff filter are given by 


S c = <*«, + K 


where 


a, = + tanh a sin 9, 

1/ — 


m = COS 0 

k 

1 . u-i 1 

a * — sinh — 
n e 


jj, k * 1, 2, 3, 


e = ripple width, 0 < e < 1* 

CHBNDP computes the poles and then calls subroutine FILTER which 
actually computes the transfer function and applies it to the frequency 
function. 




Program Listing: 


SUBROUTINE CHbNDP ( A ) 

COMPLEX A ( 1 ) »S ( 20 ) 

COMMON/COATA/ JCTR *DATA (200 ) 

COMMON/CF ILT / FO.FCOFF »NR»AMP.FFLG»S 
LOGICAL FFLG 
F FLG = • TRUE • 

GO TO 2 

ENTRY CHBSTP ( A ) 

FFLG = .FALSE. 

FO = DAT A ( JCTR ) 

FCOFF = DAT A ( JCTR+1 ) / 2. 

NR = DAT A ( JCTR+2 ) 

EPSDd = DATA ( JCTR+3 ) 

GO TO 6 

ENTRY CHLOfoP ( A ) 

FFLG = .TRUE. 

GO TO A 
ENTRY CHHIP(A) 

FFLG = .FALSE. 

FO = 0. 

FCOFF = DATA ( JCTR ) 

NR = DAT A ( JCTR+1 ) 

EPSDcJ = DATA ( JCTR+2 ) 

X = 1. / SORT (EXP ( .23*25851 * EPSDd) - 
ARG = X + SORT ( X ** 2 + 1) 

AE = ALOG(ARG) / NR 
CALL FRQFCN ( A ) 

DO 10 K * 1 *NR 

THETA = 1.5707963 * ((2.*(K + NR)-1) 
SIGK. = TANH(AE) * SIN(THETA) 

OMEGK. = COS (THETA ) 

» S(K) = CMPLX ( S I GX » OMEGM 
FFAC * COSH ( AE ) 

IF (.NOT. FFLG) FFAC « l./FFAC 
FCOFF = FFAC * FCOFF 
AMP = 1 • 

00 15 K = 1 »NR 
j AMP = AMP * CAdS(S(K) ) 

IF ( MOD (NR *2 ) .EO. o) AMP = AMP / E) 
CALL FILTER(A) 

RETURN 
END 


NR)-1)/NR ) 


EXP(.11512925*EPSDB) 
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9. Subroutine CPLOTF 


Called by: MAIN 

Calls : FRQFCN , SCALE, FACTOR*, PLOT*, NUMBER*, SYMBOL* 

Commons : blank 

Entries : none 

Description: CPLOTF is used to produce high quality plots of frequency 

spectra. The routine actually generates a data file suitable for driv- 
ing an off-line CALCOMP plotter. Since the routine embodys both equip- 
ment and procedural considerations, its use is probably limited to the 
Univac-1108 and CALCOMP plotter at Georgia Tech. It is included here 
for completeness. 

Program Listing: 


* * 


C 

c 

r- 

c * * 


SUBROUTINE CPLOTF! A) 
vt_ TH- F r L L C a'ING CONTROL STATEMENT 

FXEC'JT r 3 T a T E r N T FOR RUNS USING C&LC9 vp P L ~ T S • 


vi'ST PRECEDE T H 3 


1799 
1 7 1 f» 


400 


@USE UNIT # * TP FJ 

C0 yv 0N N » I G>A V * DELF 
COMPLEX A ( 1 ) 

CALL FROFCN(A) 

FORMAT ( ) 

.*;R I TF ( 6* 1716) 

FORMAT!- FNTFP FLO A NO 
RFA0(5.1799) FLC.FHI 
X T r c T = APS(FHI-FLO) 
lF(XT! r ST.LT,l.--3 ) GO TC 
XTEST = X TEST / ( A OR ( r _0 ) + 

IF! XTEST.LT.l .e-39) GO TO 
YSPREO = 70* 

M 2 = K/2 

NST = N2 + INT ! FLO/DELF ) 
NSP = N2 + 1 \T ( Fn 1 /OElF ) 

T 1 = l.E-35 
DO 4000 I = NST.NSP 
T 2 = C A 3 S ( A ! I ) ) 
IFIT2.GT.T1) T 1 = T2 

d^vax = 2o.*alggi ! n ) 
call scale ( O bVAX.MAXSCL ) 
XMXSCL = MAXSCL 
XKSCAL s 1 0 •**( -XMXSCL /2 ' 


PH i POR CALCOMP SPECTRUM PLOT,-) 


9999 

A 3 S ( F H I > > 
9999 


♦CALCOMP plotter routines 
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CPLOTF (Continued) 


CALL FAC TOPIC. 4) 

CALL d L0T(0.*-2-.»3) 

CALL PLOT I 1 2 • * 0 • t - 3 ) 

CALL PLOT (-1 - . »-14. .3 ) 

DO 3 I = 1 *2 

CALL PLOT ( — 10 . . ,2 ) 

CALL PLOT ( K • . 0 . »2 ) 

CALL PLOT(l% ♦ - ! 4 • ♦ 2 ) 

CALL PLOT (-1 : .Cl ,-14. ,2 ) 

CALL PLOT ( -1C . 31 , v.. l ,2 ) 

CALL PLOT ( ] 0 • 3 1 » r • 1,2) 

CALL PLOT ( U . 1 ,-14. H ,2 ) 

3 CALL PLOT ( — 1 - • ,-14.'. 1 ,2 ) 

JO 3 • 1 C I AGAIN; = 1,2 

■'0 1 , I = ' » 7 C 

Y — — 14. ~ • 2 * I 

1^ (N*00( I ,5 ) .FO. ) 00 TO 6 

call PLOT |-1 . l, v » -a ) 

GO TO P 

6 IF ( (1,1) .-0. .• ) 00 TO 7 

CALL p L r T ( -1C .16, Y, 3 ) 

GO TO ? 

7 CALL PLCT(-i:..2,Y*?) 

A CALL PLOT ( -1C . ,Y,2 ) 

1 ■ CO\T I 00 E 

Y = 2. 

00 15 I = 1,6 
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II 

CPLOTF (Continued) 


SCALE1 = 1~.*SLALL1/FLLSCL 
TEINi I F5 = 1-.*IFU.CL 
C A L r 1 = C C ALE 1 /T E m I F S 
CO A ) I = » I T E M P 

X = I*SCALC1 

IF ( VOD (1*5) .CO. ') GO TC 32 
CALL PLOT(X .-14.1 .3 ) 

GO TO 36 

32 IF (MOOtl.l'') .EO. -) GC TO 34 
CALL PLOT (X. -I a. 16, 3) 

00 TO 36 

34 CALL PLOT ( X ,-14.2 » 3 ) 

36 CALl. PLCT(X,-14.,2 ) 

4, CONTI MU •: 

00 140 I = 1 » I T f v p 

X = - I * G C L r 1 

1 = C'CO ( I ,6 ) ,:C. '■ ) TO 132 
CALL PLOT (X, -14,1,3) 

GO TO 106 

132 IF (VOD( I ,1 - ) . r 0.' 1 G~ T C 134 
CALL PL0T(X»-14.16,3) 

GO TO 136 

13- CAlL PLOT (X *-14.2*3) 

136 CALL PLOT ( X ,-14. ,2 ) 

14 v CONTIMUF 

DC 24 j I = 0,1 TEMP 
X = I * 5 C A L 1. 1 

I F ( v 0 0 ( I » 6 ) • f 0 • ' ) GO TO 2 32 
CALL PLOT (X,. 1*3) 

oO TO 236 ! 

2 32 IF(MCO< I *1 “ ) .EO. ~ > OC TO 234 • 

CALL PLOT ( X , • 1 6 *3 ) ! 

GO TO 236 

234 CALL PL OT(X, .2*3) 

236 CALL PLOT ( X , - • , 2 ) 

2 40 CONTINUE 

DO 34 J I = 1*1 TEMP • \' 

X = -I*SCALF1 

I F ( vQD ( I *6) • EG • .’ ) r '0 TC 332 

CALL PLOT ( X , • 1 , 3 ) t , 

GO TO 336 

332 !F(MOD< I ,1 J) .FO.O) GO TO 334 
CALL PLOT ( X , • 16 ,3 ) 

GO TO 336 

334 CALL PLOT(X,.2,3) 

336 CALL PlOT(X,C*,2> 

340 CONTINUE 
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CPLOTF (Continued) 


Y r I 

! - ("""Ml*-) 


) ro jr u 


1 


C 'LL P L ' T ( 

V' I O 1 ~ 
if if i •cmi, 
r ' L L p L : T ( 1 

or r • i 

17 C « L L P L - T ( \ 

• -j C •• L L =>!_.; i ( l 

. T r ■ i 


. 1 . v * 7 ) 


) '•? T C 1 7 


. 1 6 , V » ' J 

. " ♦ V * ; ) 
i ) '’ < Z l 


. /r jl: ; cl 


v. . I ' . 4 . 

As 1 = 

C A L N „ ' ’ X ( — ^ ~ *i> « — 14.5* «21» 
0 2 I = 3 ♦IPLSCL 
xnc c = -. 6 + I*A<1 

X'.F' = -.1? - T*A<: 


. ♦ - 1 ) 


C - * L L 

’■■ ! i'* - r *: 

( X [> ~ ^ 9 — 1 9 ’T * • 2 1 *1 • * I * 

. .-1 ) 


: l l 

\„ v r’ r - 

( y \ • ■ * — ] rn Z ♦ •21 ♦ — ' 1 • * I * 

, . ♦ - 1 ) 


CALL 

^ V ‘ ' i 

(-£*«?*- I A # r f • 2 1 * 



U*y\ ( 

\ k : OJ ; 

*:Cv - - r L\ T ^i “ I vl 

-Y r ?: C L r t 

( h 7 ) . ' . * 4 6 i 

call 

3Y‘-'L 

! - u • 5 9 - 1 ?•?♦•?! • 1 a ■ H F C r . 

r ' T r ^ 

:> ) 

CALL 

YJ" 2 r ~ 

(-2**~1^»^**21 tFCF\Tr 

* " • 9 ) 


CALL 

C V ' ' • l ^ * 

„ • .-> Vta. 1— 

( •»-i‘7.3*»21»0H cr 5C--L- 

= .0.»?) 


call 

\ _ * r ^ 

(,<«4* — l c .3*»21* 1 1 1 

X P * .•*,-) 


COLT I 

1 , J ", 




I ST O P 

= Yep 

- -cr l 



L, ■ '■ 

= \ r : 

- *.3T + 2 



:• - lx 

= 2 2#/ 

" f C v 


i 

; ~ 3 

I = 1 

, I c TOP 


1 

Tl = 

< C A - M 

A ( I * N 5 T ) ) ) * y x 5 C A L 


1 

I F ( T 1 

• L T • 

l«i-7) GO T 0 3' 


1 


Yl = ( 7 . . /Y c P <-:.■) ) * 4. * A LOO IP ( T 1 ) 

I F ( Y 1 • L F . -14. ) GO TO 5.) 

XI = I * PL LX - 1- . 

I F ( X 1 .LT. -1-.) 30 TO 3 0 
I F ( X 1 .'T. 1 .) 1 ; w TO 3 . 

I cr ( Y 1 . L ^ • • ) ■ ■ 0 Tu 2 F 

CALL 3 Y '* ' C L ( x 1 , . 1 6 * • 2 1 ♦ 1 -A * * ~ • ♦ i ) 

G C T p 3 
2 3 C 7 N T ! \ . . - 

CALL D l ~ T ( X 1 , - 1 4 . ♦ j ) 

CALL P L 0 T ( X 1 * Y 1 , 2 ) 

CALL PL'T (XI » - l 4 . » 7 ) 

3 COM IN J = 

CALL PL0T<n.»-2 .*-3) 

33 F 0 A * ' A T ( 1 H 1 ♦ 2 X » 1 4 H -a L T CO—LCT- ) 

.• A I T - ( 6 » 3 - ) 

Y 1 -* < C 0\T I ••0- 
-X-'TlXO 


\ 
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10. Subroutine ELFIND 




Called by: INPFOR, MAIN, PDCHK 

Calls : none 

Commons : none 

Entries : none 

Description: ELFIND compares a Hollerith string of up to six characters 

to a number of pre-stored character strings. When a match is found, an 
integer is set to a unique value which indicates the matched string. This 
is the basic operation of identifying the input commands; the integer is 
returned to the calling program and used to direct program flow to pro- 
p£tly execute the command. Tus subroutine was patterned after a similar 
subroutine in CIRCUS, but is essentially a complete rewritten version. 

Program Listing: 


C 

c 

c 

c 


SUBROUTINE ELFIND ( N A v E ♦ L ) 
PARAMETER NMAx = A2 


■RIFS TO m 
: N A M E ) AO A 
r TO THE I 
i'PE. 

iTCH ( N v AX ) 
t I ) . I = 1 • N 
»6HPK I NTF . 
,6HENl) 0> J» 
»6HCHL0/lP* 
6HS IGSFN, 

♦ 6 A P H 0 tF '■* 0 » 
1HN* 

♦ 6HI NPUTF . 

♦ NM AX 
MATCH! I ) 


ATCh T wF <f Y WORD F R 0 v an INPUT DATA 
INST ONE OF THF ALLOWABLE INPUT FORMS. 
NDEX WHICH CORRESPONDS TO THF MATCHED 


VAX ) 

6HTPLOTT ,6HTPL0TF»6hCPL0TT ♦6HC?lOTF. 
6HUW3N0P »6H TWLO ■<’? * SHEWN I P ,<SH3wp STP» 
5 HCHH I P » 6HCHBST? * 5HSYN3P * 5HSYNLP* 

F HF ROMUL * 6H »6Hl DLMUL »6H * 

1HA;-’P» 3HL I M ♦ 5 HriLOC < ♦ SHY F r , 

SH ,6HA^1D :: VO.F.HFlATSP.6HLI'-TCO, 

FHD r lFT r * 6H I NSFPT ♦ 6HREPL AC « 6 W P FOE AT / 

) 11 , 21.11 







11. Subroutine FETCH 

Called by: MAIN, PDCHK 

Calls : BCDFPT 

Commons : CFETCH 

Entries: none 

Description: FETCH is the main input routine, it reads in commands as a 

string of BCD characters, decodes the various elements in the input stream 
and stores them in array WORD. All blank characters are discarded; dif- 
ferent elements are delimited by commas. Hollerith strings are truncated 
to the first six characters and stored in WORD. Numeric characters re- 
presenting data are converted to real numbers by BCDPFT prior to storage. 
FETCH was also adapted from CIRCUS, but several changes were made. In 
particular, the program was modified to eliminate two calls to assembly 
language subroutines. 


S’J~.POJTIN r F-' T C-M •. ■'■n *LL ♦ V« : ; ) 

I N T F~ S - 3 -“.nrST \< ♦ CO ,, ‘ , A , r*crt> T , - ^ j a L * PL L - 5 , R o 1 ? f v , - r J 
I “•TF GF9 r. - - 1 * ■ U F e ? * i 5 1. 1 - r ? * •- 1 C 0 r P T tTFL'- *'\ , 0P r 5 » E 
TITL?(12) 

r ~ \ / c. F F T C h / • ■ 1 F P 2 ( S ) » :* 'J F F 1 ( 5 ) 

f‘, ; \f -\ A P 0 S T , F L A \ < , C O' 1 A • 0 r C P T » r 0 1 ' A L * * ' I M- S * *1 I f I r » ? » D L . ' 


/ H - * 1 H »l-= * 1 h - » 1 H 0 


, 1 - 1 + / 


1 

DATA LPARfu^APtN 
1 / 1 h ( * 1-1)/ 

DATA - / 1 H r / 

L = •• 

NCCLS = 8" 

Nr AD = 1 

T /-,j | c rrjr -_C - 1 .-I_r jmd.jt f-'J^P 5 ? P- T' l°\ c 

TH r V ‘p.p'i T ^ A T A IN LL CONSFCJI IN/-.’ CFLLS OF THF 
AR^AY *:»:?. H0LL er <1Th IS TTJNrAT™ TO 6 CHAP... 

] COM I \ JF 

PtAi) ( 5 * 1 0 1 » E NO = 17) ( bUF c 1 ( I ) ♦ I = 1 > 

I T - (6*1- 5) ( rO'F 1 ( I ) * I = 1 ♦ <? 0 ) 

: v = 

? < = • 


I '■-< 

Y 


NC.C 
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FETCH (Continued) 


DC 3 I = 1*6 
BUFF2 ( I ) = BLANK 

IF ( m-NCCLS ) 4 0 * 1 U 0 » 1 o 0 

EXAMINE EACH COLUMN* REMOVE BLANKS * 
SEPARATORS. 

V = M + 1 

TEST = BUFFl(M) 


TEST 


TEST 

TEST 

TEST 

TEST 

TEST 

N+l 


blank : 

COMMA 

equal 

lparen 

RPAREN 


41 * 4 * 4 1 
42*6*42 
43 . 6*43 
) 44 . 6*44 

) 45 . 4.45 


BUFF2 ( N ) = T E C T 
IF ( K ) 5*5.4 

I F TYPE HAS NOT BE 

determine IF it IS 

A 4-8 PUNCH (TITLE 
POINT NUMBER)* OR 


_LER ITH 


WORD 

WHEN 

WHEN 

WHEN 


N SET ( K = 0 ) ♦ TEST CHARACTER T " 

A DIGIT OR SIGN (NO DECISION), 
CARD), A DECIMAL POINT (FLOATING 
ONE OF THESE. IN WHICH CASE A 
SSUMED • IF K IS SET. IT V ILL 3 E 


OR NONE OF THESE. IN WH 
IS ASSUMED. IF K IS SET 
AN INTEGER 

A FLOATING POINT NUMBER 
A HOLLERITH WORD. 


NCOMMA = 



IF 

( TEST 

- NZ ) 52*51*51 

51 

IF 

( TEST 

- NINE ) 4,4.52 

52 

IF 

( TEST 

- PLUS ) 53,4,53 

53 

IF 

( TEST 

- MINUS ) 54*4,54 

54 

IF 

( TEST 

- APOST ) 55.3C.5 

55 

IF 

( TEST 

- DECPT ) 56.57,5 

56 

IF 

( TEST 

- E ) 561*560,561 

55C 

IF 

( N — 1 

) 561,561,4 

561 

K = 

2 



GO 

TO 4 


57 

K = 

1 



GO 

TO 4 
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FETCH (Continued) 


SELECT *'OD r Or CONVERSION* PASEP UPON <• 

IF ( < - 1 ) 1 * 7 *) 
r 'jFFB = aCDFPT ( .vjFF2*N ) 

IF ( \ ) 1-S.IL6.3 

aCRUL + D = -- 0 F - 3 
00 TO 91 

FNC ( A * 1 . 00 , ■.••00 O ( L + 1 ) ) ( BUFFO ( I ) , T =1 *6 ) 

L = L + 1 

l F ^PT r i f. ; i c. -i " ^ . ITH T)- r 0 '•O'' IVAr.c, r?c I N I T I *-L 17. - 

A*;0 COM I fJU 0 • 1= TrF NC'M.S C^L CCMTAIMC A 

, np-C r ;S r - T‘i r M~XT C • , . jr R I c r » cr T T^r 

N'JVNER CF /.OP OF. CO f VERTED IN LL •* RETURN. 

I c ( V-NC.CL5 ) 2 * 1 ' * 1 - 
IF ( TOOT - CO’^-'A ) 11*1*11 
LL = L 
R E T U R N 

vpvE A TULA’ CAPO INTO T‘ J E T I T LF A9RAY. 

LNC OOF ( 7 2 » 1 0 - 1 * T I T L£ ) ( u.l'FF 1 ( I ) ♦ I = 1»T2) 

50 TO 1 

a (;ARO iv* pi- .-ip r F\ PROCESSED. IF T UL LAM 
,v p m - p l a t'K SY -HOL '-'AS A CO’'"- 1 A (MCCO'A= . 5 * E l A_- 
THF N r XT CAPO. -OTHCSvISE ThFFF 10 IMF DR'-’ATl CM 
I % *£UFF ?* TO BE CONVERT E 0 * AFT' R 0 H I r , — TO 
..ILL RETURN TO T-U CALLING PP^GRA”. 

IF ( NCO^vA ) 6*1*6 

F £T C i r 0 L N L C 0 \C 0 '•’ I T -Vi T SEPARATORS OR A f £9 .*• I T r 
v o R £ THAN 15 0151 TO A NO COUL r NOT CONTINUE. 

.-.RITE (6 *2- ) ( .i 1 'r F 2 ( I 1 * I = 1 » 6 ) 

MAO = , 

RETURN 

STOP 

FCRVAT (SOU) 

FOR v A T C £ A! ) 

F OR’.' AT ( 1X8. A 1 ) 

FORMAT (?5u:.«* FETCH CANNOT CECO-r- 6Ui4h **//) 

E NO 


12. Subroutine FFT 


Called by: TIMFCN 

Calls : none 

Commons : none 

Entries : none 

Description: FFT performs the direct and inverse fast Fourier transform 

This program is substantially the same FFT routine developed under Con- 
tract NASA8-20054 and previously reported*. It has been modified, how- 
ever, to remove the FLD function, available in FORTRAN V, which appeared 
in the original version. These changes appear in the DO 10 loop, and 
the version listed here contains only standard FORTRAN- IV statements. 




\ 


*W a lsh, J. R. and R. D. Wetherington, CCS PfflBrL i nK SpSgVrtl Sfcttfl l s a.* 
Technical Report No. 7, Contract NAS8-20054, Georgia Institute of Tech- 
nology, 29 May 1970. 


% 


142 


> 



Program Listing: 

SUBROUTIME FFT < A* IGAM » I SN ) 

COMPLEX A(l).TltT2*TEMP 

DOUbLE PRECISION P 1 2 *SO *CO ♦ S I »C I »SN »CS 

P 1 2 = 6.2831b53o71795«640Do 

N = 2 ** IGAM 

NBIT = 36 - IGAM 

N1 « N - 2 

DO 30 I * 1 *N1 

IFLIP = 0 

l X = I 

00 10 J = 1 » I GAM 

1 OLD = IX 
IX = IX / 2 

IbIT = I OLD - 2 * IX 
10 IFLIP = 2 * IFLIP + loIT 
IF (I .LE. IFLIP) GO TO 30 

11 = 1 + 1 

12 = IFLIP + 1 
TEMP = A ( I 2 ) 

A ( I 2 ) = A ( 1 1 ) 

A ( I 1 > = TEMP 
30 CONTINUE 

DO 80 I = 1*1 GAM 
MEL = 2** I 
NEL2 = NEL / 2 
NSET = N / NtL 
SI = DS I N ( P I 2/ NEL ) 

Cl = DCOS ( P I 2/NEL ) 

DO 80 J = 1 »NSET 

I NCR * ( J - 1 I * NEL 

SO = O.ODO 

CO = l.ODO 

DO 80 II = 1 »NEL2 

J1 = II + I NCR 

J2 = J1 + NEL2 

T 1 = A ( J 1 ) 

T 2 = A ( J2 ) * CMPLX ( CO 9 ISN * SO) 

A ( J 1 ) = T1 + T2 
A ( J2 ) = T1 - T2 
SN s 30 * Cl + CO * SI 

CS = CO * Cl - SO * SI 

CO = cs 
80 SO = SN 

IF (ISN .GT. 0) GO TO 120 
DO 110 I = 1 »N 
110 All) * A( I )/N 
120 CONTINUE 
RETURN 
E NO 



t 

I 

i 
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13 . Subroutine FILTER 


* 


Called by: BWBNDP, CHBNDP 

Calls : FRQFCN 

Commons: blank, CFILT 


Entries : none 

Description: FILTER operates on the components in the frequency array 

to complete the computations for any type of Butterworth or Tchebysheff 
filter. Given the poles, S k , determined by BWBNDP or CHBNDP, FILTER 
calculates the transfer function 


H(f p ) 


< £ p - V (£ p - V ••• <£ p - V 


where n is the filter order (number of poles) and £ p is a normalized 
complex frequency defined by 

r f - f 

j 2 — j low pass band pass filters, 

f cutoff 

j cutof f. ^ high pags or band stop filters. 

All spectral lines subject to more than 300 dB rejection are set to zero. 
Program Listing: 


SJdROUTINE FILTtR(A) 

COMPLEX A ( 1 ) » S ( 20 ) »HD * L 

COMMON N * I GAM * DEL F tDELT 

COMMON/CF I LT/ Fu.FCOFF ♦ NR * AMP ♦ FFLG *S 

LOGICAL FFLG 

TEST * EXP '3b. / NR ) 

CALL FRQFCN ( A ) 

DO 30 I = 1 *N 
I I * I - 1 - N/2 
F = II * DELF 

FP * S I ON ( 1 • ♦ F ) * (AdS(F) - FO) / FCOrF 
IF (FFLG) 00 TO lb 
IF ( AdS < FP ) .LT. l.E-16) GO TO 2b 
FP * -l./FP 

15 IF (AbS(FP) .OT. TcST) GO TO 2b 
l - CMPLX ( 0. »FP ) 

HD = CMPLX ( 1 • » 0 • ) 

00 20 K * 1 ♦ NR 

20 HD ! HD * (Z - SCIO) 

A ( I ) * AMP * A ( I ) / HD 
GO TO 30 

2b All) * CMPLX! 0. .0. » 

30 CONTINUE 
RETURN 
END 
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14 . Subroutine FIATS P 

Called by : PROCES 

Calls : ADJN 

Commons : blank, CDATA, CDOM 

Entries : none 

Description: FLATSP loads the frequency array with components of uniform 

amplitude thus simulating the spectrum of an impulse function. It is 
useful in examing the transfer functions of filters in detail. 

Program Listing: 


SUBROUTINE FLATSP ( A ) 

COMPLEX A ( 1 ) 

COMMON N t IGAM*DtLF *DELT »PD 
COMMON /CDATA/ JCTR »DAT A ( 0 ) 

COMMON /CDOM/ DOME |_G 
LOGICAL DOMFLG 
AMP = DATA(JCTR) 

DlLF = DATA! JCTR+1 ) 

PD = l./DELF 
N = DATA ( JCTR+2 ) 

CALL ADJN 

DO 1U I = 1*N 

10 All) * CMPLX ( AMP * 0 • ) 

DOMFLG = .FALSE# 

RETURN 

END 


l 



i 
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15. Subroutine FMDEMO 


Called by: PROCES 

Calls: FRQFCN, TIMFCN 

Commons: blank, CDATA, CDOM 

Ehtries : none 

Description: FMDEMO simulates the action of an FM demodulator. Operat- 
ing on the frequency spectrum, the negative frequency components are all 
set to zero to give a spectrum characteristic of a complex time function. 
The positive frequency components are then shifted down in the data array 
by an amount corresponding to the center frequency of the demodulator, 
thus positioning the spectrum at baseband. Transforming to the time o- 
main and taking the complex logarithm of each time sample produces an 
imaginary part equal to the phase angle (modulo 2rr). A tracking loop 
corrects for excursions beyond the + rr range thus reconstructing the 
phase deviation due to the angle modulation* 



Program Listing: 


ZERO FREQUENCY 


SUbROUTINE FMDEMO(A) 

COMMON N* I GAM* DEL F *DElT .PO.CARRFQ 
COMMON /COAT A/ JCT R » DAT A ( 2 Uu ) 

COMPLEX A < 1 ) 

P 1 2 = 6.2831853 
J - 0 • 

N2 = N /2 

FO = WORD ( JCTR ) 

CALL FRQFCN(A) 

REMOVE THE NEGATIVE FREQUENCY COMPONENTS 
DO 10 I = 1 ♦ N2 
A ( 1 ) = ( U # ♦ 0.) 

MOVE THt MODULATED CARRIER TO ZERO FREQUI 

IFu = FO / DELF + .5 

NSTART = N2 + 1 

NSTuPl = N - IFC + 1 

DO 1 1 I = NSTART *N 

A ( I - I F 0 ) = A ( I ) 

DO 12 I = NST0P1.N 
A(I) = ( 0 • » o • ) 

RECOVER THE ANGLE INFORMATION 
CALL TIMFCN(A) 

A ( 1 ) = C LOG ( A ( 1 ) ) 

THETA2 = A I MAG (All ) ) 

DO 2 u I = 2 » N 
AID = C LOG ( A ( I ) ) 

THETA1 =■ A I MAG ( A { I ) ) 

THETAT = THLTA1 * THETA2 
I F ( THET AT .Lt. 0.) GC- TO 21 
GO TO 29 

I F ( AoS ( T HET A 1 ) .LE. 1.57) GO TO 29 
I F ( THET A2 .GE. 0.) GO TO 22 
J = J - 1 
GO TO 29 
! J = J + 1 
> THETA2 * THETA1 


INFORMATION 


. L E . 

0 • ) 


1.57) GO 
GO TO 22 


TEMP = 
A ( l ) = 
A ( 1 ) = 
ZERO T 
CALL F 
A ( N2 + 
RETURN 
END 


= THETA1 + P 1 2 * J 
= CMPLX ( TEMP. d. ) 

= CMPLX ( A IMAo ( A ( 1 ) ) . 
THE D-C COMPONENT 
FRQFCN(A) 

+ 1) = CMPLX (o.. 0 . ) 


* 


I 


16. Subroutine FRQMUL 

Called by: PROCES 

Calls : TIMFCN 

Commons : blank 

Entries : none 

Description: FRQMUL provides the action of a biased half-wave rectifier; 

operating on the time function, it passes only those time samples whose 
amplitude exceed a fixed threshold (currently set at 0.5 volts). The 
resulting signal is rich in harmonics of the carrier frequency. Particu- 
lar multiples can be isolated by filtering. Note that the output of 
FRQMUL is not bandlimited and the user should be aware that aliasing may 
be a problem. 

Program Listing: 


SUBROUTINE FRQMUL ( A ) 
DIMENSION A ( 1 ) 

COMMON N ♦ IGAM 
CALL TIMFCN(A) 

ThRES = .5 
NDBL = 2 * N - 1 
DO 100 I = 1 * NDoL * 2 
A ( I ) = A ( I ) - THRES 
A ( I ♦ 1 ) = 0. 

IF ( A ( I ) .LT • 0. ) All) * 0. 


100 CONTINUE 



17. Subroutine IDLMUL 


Called by: PROCES 

Calls : TIMFCN 

Commons: blank, CDATA, CFREQ 

Entries : none 

Description: IDLMUL is an ideal multiplier which operates in the time 

domain and generates the product of the signal being processed and a 
local oscillator signal. 

Program Listing: 


SUBROUTINE IDLMJL(A) 

COMMON N » ISAM. DELF .CELT 
COMMON/ COAT A/ JCTK .DATA ( 2 «u ) 
COMMON / C F R E 0 / NFR . FR ( 6 ) 

COMPLEX A ( 1 ) 

AMPLO = DAT A ( JCTR ) 

FLO = DATA! JCTR + 1 ) 

CALL TIMFCN(A) 

P 1 2 = 6.2831db3 
wLO = P 1 2 * FLO 
DO 1 I = 1 »N 
11=1-1 
T = I I * DcLT 

All) = All) * AMPLO * SIN I WLO * 

RETURN 

tND 




18. Subroutine INPFOR 


Called by: MAIN 

Calls: ELFIND 

Commons : CWORD 

Entries : none 

Description: INPFOR is a service routine that will list the input format 

and define the parameters of any block input consnand. It is added as a 
convenience to the remote terminal user; it has no affect on the circuit 
or signal being processed. 

Program Listing: 


SUBROUTINE INPFOR 
COMMON /CWORD/ W0kD(1u 

logical flg 

FLG = .TRUE. 

CALL ElFIN0(W0R0(2 > *L ) 


GO TO ( 

1* 

1* 

1 * 

1* 

1* 

1 » 

1* 

1* 

9* 

10* 

A 

11 ♦ 

12* 

13 » 

14* 

15* 

16* 

17. 

18. 

19* 

20. 

6 

21* 

1* 

23 ♦ 

1* 

2b * 

26 * 

27* 

26* 

1* 

1* 

C 

1 » 

1* 

1 * 

34 • 

35* 

1 * 

1* 

1* 

1* 

1* 

D 

1 » 

1* 

1 * 

L 








*R I TE ( 6 * 700 1 ) *UR0(2) 
oO TO 999 
WRI TE ( 6 * 7020 ) 
wRI TE ( 6 » 7 1 v4 ) 

IF (FLG) GC TO 999 
wRITE(6»703b) 

*RI TE(6*711 1 ) 

IF (FLG) GO TO 999 
WRI TE ( 6 ♦ 70o9 ) 
wRI TE(6*7101 ) 

IF (FLG) GO TO 999 
WRI TE (6*701 u) 
wKlTE(6*71u2) 

IF (FLG) GO TO 999 
WRlTt(6»70ll ) 
wRITE(6*71v2 ) 

IF (FLG) GO TO 999 
wRI TE (6*7C12 ) 

WRI TE (6*7101 ) 

IF (FLG) GO TO 999 
wRITE( 6*7G13> 
wRITE(6.7lGl ) 
*RITt(6»71u3) 

IF (FLG) GO TO 999 


150 


i 


( I 

INPFOR (Continued) 


lu 

•; R I T - <6 , 7 'U ) 
■•7 I T r ( 6 ♦ 7 1 ? 1 

• ' R I T c (6*71' 3 ) 



IF ( r L 0 ) 90 T; 

o 9 9 

IF 

•: R I T - ( 6» 7, 'l t - ) 
>.RIT r (6*71 ,2 ) 
v R I T r ( 6 » 7 1 .. ? ) 



I- (-L9) oC TO 

9 9 'i 

16 

R I T r ( 6 ♦ 7 16! 

<\'x IT” (6*71 .. 1 ) 
... R I T 0 ( 6 » 7 1 . .?■ ) 



IF (FLO) 00 TO 

099 

17 

IT rr (‘S»"’'l 7 ) 
.. R I T ( 6 , 7 1 1 ) 



IF (FI.G) <~.r 1~' 

099 


'.•■RI T- ( 8 ) 

• k I T : ( 6 » 7 1 ' ' ? ) 



IF (FLO) 00 TO 

9 q \) 

1 ? 

;• R I T F ( 6 , 7 . ' 1 9 ) 
V. R I T r ( 6 * 7 1 ? ) 



IF (FLO) r-z TO 

'9 9 9 


2 o a R I T L ( 6 ♦ 7 , 2 F 1 
.• RI T- ( 6 ,71 . 7 ) 

I c r -C TO 999 


Ti ,v R I T E ( 6 » 7 r 4 ) 
•"? I T* ( f, .71 . 7 ) 



INPFOR (Continued) 


F "D'«M ( ) 

FCR’-'AT ( IX ♦ 
FORMAT ( /- 
FORMAT ( /- 
FORMAT ( /- 
FORMAT ( /- 
FORMAT ( /- 
FORMAT ( /- 
FORMAT ( /- 
FORMAT ( /- 
FORMAT ( /- 
FORMAT ( /- 
FORMAT ( /- 
FORMAT ( /- 
FORMAT ( /- 
FORMAT ( /- 
FORMAT ( /- 
FORMAT ( /- 
FORMAT ( /- 
A 

FORMAT ( /- 
FORMAT ( /- 
FORMAT (/- 
FORMAT ( - 


A6»- IS MOT A M INPUT COMMAND-) 

BWfcNDP* FI * B'W* NR-) 

EWLOWP* FC. NR - ) 
b VH I P ♦ FC* NR- ) 

6 WB S TP* F 0 » B'W, NR- ) 

CH5NDP, F.., ;.:W, NR* EPSDB-) 

CHLOWP* F'C* NR* C.P5D3- ) 

CHHIP* FC‘* MR* EPSDB-) 

CH5STP* Fv, dW* NR* EPSDB-) 

SYNBP* F", * BW* MR-) 

SYNLP* FC* NR- ) 

SYNH D » F C ♦ MR-) 

SIGGEM* F.;, FMOO * AM* PM* E v » A- ) 
FRQMUL-) 

IDLMUL* A|.n t FLO-) 

F V D E m 0 * F 0 — ) 

PHOEMO* Fwi-) 

AMp» GAIN-/- CAIN = VOLTAGE GAIN* D.,- 
(6 OB = FACTOR OF 2)-) 


71o5 FORMATt- 
71 jS FORMAT ( - 
A 

7107 FORMAT ( - 

7110 FORMAT (- 
A 

B 

711 1 FORMAT <- 
A 

6 - 

7112 FORMAT!/ 
999 WHITEC6. 

RETURN 

END 


L I M * CL* CH. GL-) 

AM[)E v O * F j — ) 

FLATSP* AMP, OELF* N-) 

p- _ CENTER FREQ* HZ-/- BW = BANDWIDTH* HZ-/ 

. MR = NUV.3 C P OF SECTIONS-) 

PC = CORNER F R E G * HZ-/- NR = NUM9FP OF SECTIONS-) 
PPSDR = CHEBYSH C V RIPPLE FACTOR* Do-) 

Fo"s CARRIER FREQ. HZ-/ 

F MO D * MODULATION FREC* HZ-/ 

AM = PERCENTAGE AMPLITUDE MODULATION-/ 
pvi = PEAK PHASE DEVIATION* RADIANS-/ 

FM . PEAK FREQUENCY DEVIATION. HZ-/ 

A = pfAK AMPLITUDE, VOLTS-) 

(NO PARAMETERS)-) 

ALC - PEAK AMPLITUDE OF LO SIGNAL* VOLTS-/ 

FLO = FREQ OF • LO. HZ- ) 

FO = CENTER FREG* HZ - ) 

CL = LOW CLIPPING LEVEL* VOLTS-/ 

CH = HIGH CLIPPING LEVEL* VOLTS-/ 

GL = LIMITER GAIN* VOLT S/VCLTS- ) 

AMP = AMPLITUDE OF SPECTRAL LINES-/ 

0ElF - prctq SEPARATION OF lines. HZ-/ 

N = ARRAY SIZE-) 

- INPFCR cannot decode — * a 6 / ) 

70UC ) 


(I 
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19. Subroutine LFOLD 

Called by: TIMFCN 

Calls : none 

Commons : none 

Ehtries : none 

Description: LFOLD provides the action of folding and unfolding the fre 

quency spectrum to meet the requirements of the FFT. ftie frequency do- 
main representation is always ordered by frequency except when entering 
or leaving the FFT. 

Program Listing: 


SUBROUTINE LFOLD (A.N) 
COMPLEX A (1).T1 
N2=N/2 

DO 10 I = 1 *N2 

‘ I I = I +N2 

i T 1 = A ( I ) 

i A ( I ) = A (II) 

l 10 A ( I I) =T1 

f RETURN 

i END 



i 


r 



20. Subroutine L1M 


Called by : PROCES 

Calls : T1MFCN 

Commons : blank, CDATA 

Entries : none 

Description: LIM operates on the time function to produce ideal limit- 

ing. Signal excursions are clipped to specified upper and lower limit 
levels. Note that the output of LIM is not bandlimited and the user 
should be aware that aliasing may be a problem. 

Program Listing: 


SUdROUTINE LINK A) 

COMMON N 

COMMON /CD AT A/ JCTR .DATA ( 2JU ) 
CLEVL = DATA ( JCTR ) 

CLEVH = DATA! JCTR+1 ) 

GL = DAT A ( JCTR + 2 ) 

COMPLEX A ( 1 ) 

CALL TIMFCN(A) 

DO 1 I * 1»N 
A ( I ) = GL*A( I ) 

IF ( RE AL ( A ( I ) ) .Lt. CLEVL) A(I) 
IF (REAL ( A ( I ) ) .Gt. CLEVH ) A(I) 
RETURN 
END 


CMPLX(CLEVL»0. ) 
CMPLX (CLEVH *0« ) 


♦ 


(! 


21. Subroutine LSTCOM 

Called by: MAIN 

Calls : none 

Commons : CWORD 

Entries: none 

Description: LSTCOM is a service routine that will list all of the valid 

commands that are recognized by FATCAT. It is included as an aid to the 
remote terminal user; calling LSTCOM has no affect on the circuit or sig- 
nal being processed. 

Program Listing: 


3 JFP'O'T I V~ t3TC0 v 

/CWORO/ ( 1 ) 

.•:RI T- ( S * 7 ^ 1 ) 

I = ’ 

IF (ajR'U?) 0. 3M-'LL) 
I = ~ 


if (.vcprm 
if ( vjR:: ( 2 > 
IF (A'CRD(2) 
1= ( .'/CRT ( 2 ) 

I - Ca’0"^(2) 

IF C.OR r M2) 
IF l.\CR''(2) 


eg. bhSoi'^CD go to ?c 
eg. 6h- i lt - 1 1 go to 3 0 
j. 6 ' 0 [''c.)j ) oo r~ <*c 

r .~. 4 - 1 v I s C ) GC T ~ 5 ^ 
-0. 6-tCOVTPC) GO T~ Ff 
“0. 5^\ C !7F) GO Tj 7, 

-'3. 6MLIFTCD GO 8 0 



aR I T c 

( 6 f 7"? 

~ ) 

0 7 

D ( 2 I 


GO TO 

999 





2 ' 

vR I TF 
IF (I 

{ 6 * 7 
• L Cv • 

2 ) 
) 

r~ r\ 

\J 

TO 

99 9 


WRITE 
IF (I 

(6*7*. 

* . vw • 

7 ) 


T* 

9 9 9 

A .j 

•■RITE 
IF (I 

( 6 ♦ 7 

.EG. 

A ) 
. ) 

V.vJ 

TO 

999 

c. 

y J 

RITE 
IF (I 

( 6 > 7 : 

♦ c . * 

r o ) 

) 

3': 

TO 

0 99 

A ^ 

a.'RI Tf 

1= (I 

( 6 ♦ r 

. r 0 . 

6) 

) 

o-o 

T 0 

999 

7 ; 

'.'RITE 

(6.7' 

7) 





GO TO 

9 99 





8 

•.RITE 

( 5 ♦ 7 ' ' 

8 ) 





99? R F T JR\ 

7..1 FCRVATt /1 C X, -FATCAT C ^ 'O S0''"ARY- ) 

7_2 FCR^ATl /I .X »-GC .'FC-: 0 l G G i. ♦ 3 X » - ? I o f l A (_ GcNr.R-T 
A - FLAT5P-»5X t-FL-'T SPEC*'**?. GFV?.\T 9-) 
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LSTCOM (Continued) 


7GG3 FORMAT ( / * luX .-FILTERS-//- BUTTERWORTH — / 

A - 3WENDP-.5X.-BAND PASS-/ 

6 — 3WLOWP — * 5X * — LOW PAS-S—/ 

C - BWHIP-.6X. —HIGH PASS-/ 

0 - BWBSTP-.5X .-3AND STOP-/ 

E /- TCHEBYSHEFF— / 

F - CHBNDP-.5X.-3AND PASS-/ 

G - CHLOWP-.5X .-LOW PASS-/ 

H - CHHIP-.6X.-HIGH PASS-/ 

1 - CH6STP-*5X .-BAND STOP-/ 

J /- SYNCHRONOUSLY TUNED — / 

< - SYNBP- »6X »-BAND PASS-/ 

L - SYNLP-*6X»-L0'a' PASS-/ 

M - SYNHP- ,6X*-HIGH PASS-) 

7004 FORMAT ( /*! LX ♦-DEMODULATORS — // 

A - FMDEM0-.5X »-FM DEMODULATOR-/ 

3 _ amdevo- *5x .-am demodulator-/ 

c - PrlDEMC-*5X » -PHASE DEMODULATOR-) 

70u5 FORMAT ( /.1CX»- MISCELLANEOUS -//» 

A - FRQMUL-.5X. -FREQUENCY MULTIPLIER-/. 

C - IDLMUL-.5X, -IDEAL MULTIPLIER-/ 

D _ AMP-. 8X .-AMPLIFIER-/ 

E - LIM-.8X , -LIMITER-) 

7006 FORMAT (/ »1GX .-CONTROL COMMANDS-// 

A - PR I NTT- . 5X .-PR I NT TIME FUNCTION-/ 

B _ PRINTF-.5X. -PRINT FREQUENCY FUNCTION-/ 

C - TPLOTT-. 5X .-PR I NTEP PLOT CF TIME FUNCTION-/ 

D . TPL0TF-.5X .-PRINTER PLOT OF FREQUENCY FUNCTION-/ 

F - CPLCTT-.5X .-REMOTE PLOT OF TIME FUNCTION-/ 

F "PLOTF— . 5X *-REMOT c PLOT OF FREQUENCY FUNCTION-/ 

0 I 3L0C<-»6X» -PROCESS SIGNAL TO OUTPUT OF BLOCK SPECIFIED-/ 

H - PRIMEF-SX.-LIST PRIME FACTORS OF ALL SOURCE FREQUENCIES-/ 

1 - FND OF JOE-. IX. -TERMINATES RUN-/ 

J - I N PUT FORMATS » BLOCK N Am £ LIST INPUTS FOR NAMED BLOCK-/) 

7007 FORMAT!/. 4X. -SPECIAL COMMANDS TO SPECIFY ARRAY SIZE-// 

A - YES-. 8X. -LISTED ARRAY SIZE ACCEPTABLE-/ 
a - N 0 -. 9 X.-LISTED ARRAY SIZE NOT ACCEPTABLE-/ 

0 - N-.10X.-SET ARRAY SIZE TO SPECIFIER VALUE-/) 

•trs'Q FORMAT!/.- LIST COMMANDS- COMMANDS ARE LISTED BY GIVING-. 

iSJ WORD SETS SEPARATED ST « CO». -. 

•3 -THE FIRST IS -.1H-.-LIST COMMAND .-. 1H-/ . - POSSIBLE . 

C - SECOND WORDS AND THE RESULTING OUTPUTS ARE — / 

D - ALL-. 12X. -LIST ENTIRE COMMAND SET-/ 
c - SOURCE 5 - . BX »-L I ST COMMANDS FOR SOURCE BLOCKS-/ 

F - FILTERS-.8X.-HST COMMANDS FOR FILTERS-/ 

G - DEMODULATORS-. 3X .-L 1ST COMMANDS FOR DEMODULATORS-/ 

H - MISC-.11X.-HST OTHER BLOCK COMMANDS ( AMP , L If .ETC ) -/ 

I - C0NTRDL-.8X .-LIST CONTROL COMMANDS-/ 

J - NSIZE-.1CX.-LIST COMMANDS CONTROLLING ARRAY HZt-/ 
i - l l ST COMMANDS- .2X .-LIST THE ABOVE INFORMATION-/) 

7020 FORMAT!/- LSTCOM CANNOT DECODE -»A6/) 

END 
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22. Subroutine PDCHK 



Called by: MAIN 

Calls: P ERICH), FETCH, ELFIND, ADJN, PRTFAC 

Commons: blank, CWORD , CFLGS 

Ehtries : none 

Description: PDCHK is used in checking the period of the data set and 

ascertaining that the Nyquist criterion is met. This subroutine is 
largely executive in nature; the major calculations are carred out by 
subroutine PERIOD. PDCHK is called when processing is called for, pro- 
vided a new source frequency has been added since the last processing 
call. On its initial call, it determines the period of the data set, 
and the associated Af, At, and array size. On any subsequent calls 
(which happen only if a new source frequency has been introduced), it 
checks to see if all frequencies are still periodic on the established 
period. If they are, processing continues; if not, an error message 
is written and the run terminated. 

Program Listing: 


S'JBOOi.'TI \ r PPC.FK 

CDMVCN *■: * T o n- _ F * r L T ♦P"' 1 

cr N vv"'n /o.'T'pn/ vofoii' ) , 

C-T^CN /erLT-S/ POFLu.ARFLG * 

LOGICAL POf- L-‘ ,i-FLo . 

‘!FLC: = o 
\ 5 F T = N 
CALL PER I 00 
IF(APFLC) C,C TO 5 > 

\SFT = N 

■. R I T r ( 6 ♦ 7 t-. . ) P 0 ♦ 0 E L F 

•;R I T F ( 6 ♦ 7 ' 1 ) M 1 1 G A v * 0 r L T ^ 

call fetc>-'C.-cr?*l *. r ) 

CALL FL^r^T-'FPO.LTVP ) \ 

ifiltyp . c g. 3 ;) 3 ' t ~ «:■" ; 

IF(LTYP . c 0. 31) 30 TO 3 j , 

IFILTYP »E 0. 32) CO T ~ a: 

IFILTYP . r :.. 6) STOP ! 

a'R I T r I 6 ♦ 7 ' 2) 

(>0 TO ? 
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PDCHK (Continued) 


%\> WRITE <6 *7003) 

GO TO 2 0 
4C N = WORD (2) 

NFLC- = 0 
CALL AD JN 

I F ( N . G FI • NS FT ) GO TO 10 
NFLG = 1 
GO TO 900 

5 0 I F ( N5ET .GE. N ) GO TO 990 
W R I T E < 6 ♦ 7 0 C 5 ) N 
CALL PRTFAC 
WRITE (6 *7006) 

•STOP 


1 \J • 


IF ( NFLG 


• r \ o. « 


) GO TO 999 


WRITE (6*7004) N » N S E T 
WRITE (6*7'' l ,3 ) 

GO TO 2 0 
990 N = NSET 
999 PDFLG = .FALSE. 

ARE* LG = . T R U f . • 

700^ FORMAT!/*- PERIOD = -*1PE10.3,- SECONDS »-* -DE LTA-F - -♦ 

7001 FORMATION = - » 1 6 * - » IGAV. = -»I2*- » DELTA-T 

, .1PE10.3/- IS THIS SATISFACTORY-) 

7002 FORMAT!- INPUT MEANINGLESS * ENTER YES* NO* OR f- . VAL* JE- > 

7003 FORMAT!- "NTER N* VALUE-) 

7004 FORMAT!- N = -.16.- UNACCFPT ABLF ** N V UST- 

7004 F am - be-, 16,- TO MEET NYQUIST CRITERION-) 

7005 FORMAT!- SOURCE FREQUENCIES RE Q U * RE SIZE 

1 -TO BE -*I9//— PRIME FACTORS ARE-/) 

7006 FORMAT!//- RUN IS BEING TERMINATED-) 

END 


X 



23. Subroutine PERIOD 


•» 


I 


Called by: PDCHK 

Calls: ADJN 

Commons : b lank , CFREQ 

Entries : PRTFAC 

Description: PERIOD factors each source frequency into prime factors 

and constructs the highest common factor to determine Af and the asso- 
ciated smallest period on which all of the frequencies are periodic. 
From the period and the highest frequency present, the number of sam- 
ples to meet the Nyquist criterion is computed; this number is adjusted 
upward (if necessary) to a power of 2 by ADJN. 

Entry PRTFAC will produce a listing of the prime factors of the 
source frequencies. 

Program Listing: 



PERIOD (Continued) 


’ACTOR ARRAY 


INTO PRIME FACTORS. 


IDEL 


• NE • 


V A X F A C = 0 

30 100 I = l.NFR 

CONVERT ITH FREC TO INTEGER. 

IFR = FR ( I ) 

CLEAR ITH COL OF FACTOR ARRAY IA. 

00 10 J = 1.30 

1 A ( J , I ) = C 
J = 1 
ITFST = 2 

I OF L = 1 

FACTOR ITH FRFQ INTO PRIME FACTORS. 

GO TO AO 

ITEST = ITEST + IDEL 
IDEL = 2 

IF ( IFR .EQ. 1 ) GO TO 90 
ALT' 1 = SORT ( IFR) 

IF ( ITEST .GT. ALI- V > GO TO 60 
IF(M3D( IFR. ITEST) • NE • 0) GO TO 30 
I A ( J ♦ I > = ITEST 

J = J + 1 

I*(j .GE. 30) GO TO 990 
IFR = IFR / ITEST 
GO TO 50 
i I A ( J . I ) = IFR 

1 IF (j , GT • MAX FAC ) VAXFAC = J 
; CONTINUE 

if (Flag) go to 1^2 

lim = VAXFAC + 1 
DO 101 I = l.LIM 

L WRITE(6»711) I » ( I A ( I »J)*J*1*NFR) 

GO TO 999 

*■ ALL FR'OS FACTORED. PRI^F FACTORS I* 
* OF IA. FIND COMMON P AC TORS AND PLAC ! 

2 I DC T = 1 

DO 110 I = l.NFR 
0 ICT ( I ) = 1 


VAXFAC 


FACTORS 


FTRSI NFR 
IN I OUT. 


rOLb 


I DC T = 1 

DO 110 I = l.NFR 
I CT ( I ) = 1 
I MAX = 0 

DO 130 I = l.NFR 
IF(IA(1.I) .GT. IMAX) IMAX 
ADVANCE ALL ARRAY POINTERS 
DO 200 I = 1 *NF R 
J = ICT ( I ) 

I F (I A ( J* I ) »EQ. o ) GO TO A 
I p ( I A ( J . I ) .OF. IMAX) GO Ti 
ICT ( I ) ■ J + 1 
GO TO 140 
CONTINUE 


* I A( 1, I ) 

TO A FACTOR 


I MAx 


160 


fh 


V****>. 


FERIOD (Continued) 


r * * if all poi;.T'-~ -act 09 5 are not i*'ax» at.vance ivax. 

21 - OZ 22 - I = 1 L,- 

J = ICT ( I ) 

I P ( I .A ( J * I ) . % ■ . r : A V ) 0 0 TO 3 u C 

?:• , C - N T I v J E 

r 4 , * ALL ‘-'ACTOR’S A'- I'*AY. T9iv^"’ TP IC.'T* ADVANCE ALL 
X * # s-T._ T c- 5 - ryr^, AM- ? F T r 5 T Fpr.' CP 'POL FACTOR* 

I pi;T ( I DC T 1 - ! *' a X 

I X T = I OCT + 1 
>■> ' I = l*’)F~ 

-> ? i c t ( i ) = i c t t n + 1 

TO 2 1 - 

I -1 AX = I A ( J, I) 

O T 'T 13 3 

I XT = I O C T - 1 

~ l f = 1* 

, » I - \ ACT; 11 l\ I 0 J T , C K L F = 1 . "THvIR '■ IS'- D - L F = 

- v 1 \r j ~c ILL STI'.'f S I *■' I0"T# 

I f (I : C t . t - . ) TO Air 

0*> A1 ' I = 1 * T XT 

U ] ' ’ , - l - = 1_ F *■ I r 1 J T ( I 1 

L\- o-'i = 1 . /PF-LF 

F “ A X = r . 

;vj ^2 : 1 = i*nff 

42 I - ( - A? ( I ) *CT. F f -’AX) A VAX = FRU) 

q * * S A v p l I N 0 PERIOD T (. '/EFT NYL-I.-.-I C R I I F •'< I ON • 

% = 2. * PC * ~ ' A X + .3 

r # * ARRAY SIZE ANP SA-P PERIOD OILT) TO SATISFY - P T. 

call ADJN 

r, 0 t C R 9 9 
99: R I T - r ( A , 7 ) I 

STOP 

7 FORMAT ( / / 1 H * - F R IM. 1 E NC Y -* 12*-HA5 v 0RE THAN 79 FACTORS-/) 
711 FOR VAT ( 1 h *M IX » 16 >) 

-•;o return 


24. Subroutine PHD DIO 


% 


(\ 


Called by : PROCES 

Calls : FRQFCN, TIMFCN 

Commons : b lank , CD ATA 


Entries : none 

Description: PHDQIO simulates the action of an ideal phase demodulator. 

Operating on the frequency spectrum, the negative frequency components 
are all set to zero to give a spectrum characteristic o£ * ** 

function. The positive frequency components are then shifted down i 
the data array by an amount corresponding to the center 
demodulator, thus positioning the spectrum at baseband. Transforming 
to the time domain and taking the complex logarithm of each time 
produces an imaginary part equal to the phase angle (m^ulo 2n) . A track 
ing loop corrects for excursions beyong the + rr range thus reconstruct 
ing the phase deviation due to the angle modulation. 


\ 
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Program listing: 



SUBROUTINE PhDEMO ( A ) 

COMMON N » IGAM.DELF » DELT .PO.CARRFQ 
COMMON/CDAT A/ JCTR .DATA ( 2JH ) 

COMPLEX A ( 1 ) 

P 1 2 = 6.2831653 
J — 0 • 

N2 = N / 2 

FO = aORUUCTR) 

CALL FRQFCN1A) 

C * * REMOVE TmE NEGATIVc. FReuUENCY CGMPONeNTS 
DO 10 I = 1 » N2 
luAU) = <U.»o.) 

C * * MOVE THE MCDULATcD CARRIER TO ZERO FREQUENCY 
I Fu = Fo / DELF + .5 
NSTART = N2 + 1 
NSTOP1 = N - IFO + 1 
DO 11 I = NSTART. N 

11 A ( I - IFO) = A ( I ) 

DO 12 I = NST0P1.N 

12 A ( I ) = (u. > 0. ) 

c * * RECOVER The AfooLE INFORMATION 
CALL TIMFCN(A) 

All) = CLOGIh l 1 ) ) 

THE T A2 = A I MAG ( H ( 1 ) ) 

DO 2<J I = 2 ♦ N 

A ( I ) = CLOG ( A ( 1)) ! 

THETA1 = A I MAG ( A ( I ) ) 

THETAT = THETA1 * THE T A2 ' 

IF ( THETAT .LE. 0.) GO TO 21 
GO TO 29 

21 IF ( A6S ( THETA1 ) «Lt. 1.57) GO TO 29 
I F ( THET A2 .GE. 0.) GO TO 22 

J * J - 1 
GO TO 29 

22 J * J + 1 

29 THETA2 s THETA1 

TeMP * THETA! + PI2 * J 
2j A( I ) * CMPLXl TEMP. w, ) 

All) * CMPLXl AIMAG1A ! 1 )) * ^.) 

C * * ZERO THE D-C COMPONENT 
CALL FR0FCN1A) 

A(N2+l) a lC..O.) I 

RETURN j 

END 



I! 


25. Subroutine PROCES 
Called by: MAIN 

rails* SIGGEN AMP BWBNDP , BWLCWP, BWHIP, BWBSTP, CHBNDP, CHLCWP, CHHIP, 
OTBSW', SYNBP , SYNLP , STOP. FRQMUL, IDIMUI, LIM, FMDEMO, PBDEMO, 

AMDEMO, FLATSP. 

Commons : none 

Entries: none 

Description: PROCES is simply a switching routine that calls the proper 

block model for processing the signal. 


* 


I! 


Program Listing: 


SUBROUTINE PROCtS ( I * A ) 

COMPLEX Ail) 

(jO TO ( 1» 2, 3, 4* 1? * 6 » 7» a* 9, 10, 

1 11, 12,13, 14, 15, 16,17, 18,19, 20, 

2 21, tl , 23), I 

1 CALL SIGGEN(A) 

RETURN 

2 CALL AMP ( A ) 

RETURN 

3 CALL bWdNDP(A) 

RETURN 

4 CALL 6v»L0WP(A) 

RETURN 

3 CALL BwhlP(A) 

RETURN 

6 CALL BWBSTP(A) 

RETURN 

7 CALL CHoNDP(A) 

RETURN 

8 CALL CHLOWP(A) 

RETURN 

9 CALL CHHIP(A) 

RETURN 

lu CALL CrlBSTP ( A ) ' 

RETURN 4 

11 CALL SYNBP(A) ? 

RcTURN | 

12 CALL SYNLP(A) 

l RETURN ; 

i 13 CALL SYNHP(A) 

! RETURN 

' 14 CALL FRGMUL(A) 

15 RETURN 

16 CALL IOLMUL(A) 

RETURN 

17 CALL LIM(A) 

la RET URN 

19 CALL FMOcMU(A) 

RETURN 

f 2 J CALL PHDEMu(A) 

*■ 21 RETURN | 

$ 22 CALL AMOEMO(A) 

RETURN 

23 CALL FLATSP(A) , , 

RETURN i 

END 


% 

\ 
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n n 


* 


{ 


26. Subroutine SCALE 

Called by: CPLOTF, TTFP 

Calls : none 

Commons : none 

Entries : none 

Description: SCALE establishes the (floating) ordinate scaling for spec 

trum plotter routines. 

Program Listing: 


* # 


1 


10 
1 1 


99 9 


subroutine scale ( obmax .maxscl ) 

* This SUoRCUTINt ESTAoLlSrlES ORDINAT 
THF REMOTE SPECTRUM PLOTTER. 


IF ( OoMAX.Ll.u. ) GO TO 1^ 
MAXSCL = 0 

MAXSCL = MAXSCL + 10 
DIFF = DoMAX - MAXSCL 
IF(OIFF.GT.G.) GO TO 1 
GO TO 999 
MAXSCL = 0 

MAXSCL = MAXSCL - lo 
OIFF = (OBMAX - MAXSCL) 
IF (OIFF.LE.O. ) oO TO 11 
MAXSCL = MAXSCL + 1 ^ 
RETURN 
END 


E 


SCALING for 
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27. Subroutine SIGGEN 


Called by: PROCES 

Calls : none 

Commons: blank, CDATA, CDOM 

Entries : none 

Description: SIGGEN generates and stores in the data array samples of 

the time function for the specified carrier frequency, modulation fre- 
quency, modulation types, and modulation indices. Modulation types in* 
elude AM, FM, and PM; AM can be used in combination with either of the 
other two. At least one of the indices for FM and PM must be zero. 
SIGGEN flags the data array as containing a time function. 

Program Listing: 

c j _* ■q.'MjT I ,v' 5 I 'r-; \ ( ' 1 
C^vvn.v/CDATA/ JCTR*OATA ( 2 > 

C ■V.IV'VJ /COC''-' / PC’r L r 
C 0'-’ v C \ ■:* I G ■- “ * ~ L f- » r LT 
C ? v P L i X A ( 1 ) 


logical 

_ 

^ * v j 

” L’ 




= 

• Ti 

R‘J r • 



r - 

^ , 

iT- 

•• ( J 

r. T p ) 



r VC:) 

= 


AT A 

( JCTP 

+ 

PC J A * 

■/ , 

= 

OAT 

A( JC 

T 

R H 

P<Ph! 

v 

= 

[' m \ 

T A ( j 

r r x 


PV 

= 

r\ A 

TA ( J 


T ri 

A PS 

= 


r a 

( JCT 

K 

+ 

P 12 

= 

6 • 


!r:’ 



.VC = 

P I 2 


FC 



~ 

D 

12 

* 

F'TI 



; T A 

.V-* 

= 

cr 

T " 

/ 

1 

1 F ( 2 

<o 

H r > 

V • 

r-T. 

• 

1 

I - (p 

< F 


V • 

r,T. 

• 

1 

r 0 1 

I 

= 

1 f 

V 



I I = 

I 

- 

1 




T = 

I I 

# 

0 L 

LT 



C 0 F K 

= 

A 

MDC 

> * ( 


1* 

M 's kj 

3 

VC 


T + 

K 

<p 

A ( I ) 

= 


■•’PL x ( ( c : 

■ F 


T ~ 7 > 



•0 

,-• F T A 

rV = P/CT-V / F V"- 




' ~ 2 

I = 1* M 




1 ! = 

1 - 1 




T = 

I I * OF L T 




C OC = 

s V p c * ( 1, +• -FT AAV * 

c ( 



A\G 

s WC * T ♦ ».FTAF-.‘ * ' I '• ( “ 

■ # T 


2 

A (I J 

= C M PLX ( ( CC£r * FI \ < ' ' 1 j ) 

) * • 



T 

0 PRR 



5 .. 

•; x I T 

F ( 6 ♦ 1 M » 

-::TH 

$ 

M 

r -;qv 

AT(lr ,-TIvf FC\ r «Rr-< - 

* 

If. 


Ax 

- C.PF C I r I r ' - ) 



* ! ) ) 
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28. Subroutine STRDTA 


Called by: MAIN 

Calls : none 

Commons : CDATA, CFREQ, CFLGS, CCIRKT. CWORD 

Description: STRDTA stores the input parameters of all circuit blocks 

in permanent storage. The incoming data is removed from array WORD 
(placed there by FETCH) and stored in the basic data storage array DATA. 
Storage is dynamic, with the data location specified by JCTR. The type 
of block and the associated value of JCTR is stored external to this 
routine in the two column array ITYP. 

Program Listing: 


SUBROUTINE STRUT a(K1*K.2»K3) 

COMMON/ CD AT A/ JCTR *DATA < 20u ) 

COMMON /CFREQ/ NFR *FR ( 6 ) 

COMMON /CFLGS/ PDFLG ♦ ARFLG 
COMMON /CCIRK.T/ NoLK. ♦ I T VP ( 30 *2 ) 

COMMON /CWORD/ WORD(IC) 

LOGICAL pdflg 
JCTR = I TYP ( NBLK + 1*2) 

DO 10 I = 1»<1 

10 DAT A ( JCTR ♦ I -1) = W0RD(I + 1) 

JCTR = JCTK + K 1 
I F ( <2 .EQ. 0) GO TO 999 
PDFLG * .TRUE. 

DO 20 I = l*K2 
TEMP = WORD ( K.3 ♦ I ) 

IF (TEMP .LT. 1.) GO TO 20 
NFR = NFR + 1 
I F ( NFR .GT. 6) GO TO 50 
FR(NFR) * TEMP 
20 CONTINUE 
GO TO 999 
50 WRITE(6»70vu) 

STOP 

7000 FORMAT ( /- INPUT FREQUENCIES EXCcED SIX-/) 
999 RETURN 
END 
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29. Subroutine SYNBP 

Called by: PROCES 

Calls : FRQFCN 

Commons: blank, CDATA 

Entries: SYNLP, SYNHP 

Description: SYNBP provides a model of a synchromously tuned bandpass 

filter. Entries SYNLP and SYNHP provide low pass and high pass models. 
The filter transfer function is computed to be 


(r^r) 


where 


n * number of sections. 


f/f' for low pass, 
co 

- f' /f for high pass, 
co 

(f - f )/f' for bandpass, 
. o co 


f' = f /$(n), 
co co 


and f is the nominal corner frequency in Hz. For this idealized filter, 
the insertion loss at the center of the pass band is zero, and the atten- 
uation at the corner frequency is 3 dB. The computed transfer function 
for the bandpass filter is symmetrical about center frequency. 


7 


I 


Program Listing: 


* 


SUbROUT I NE SYNBP(A) 

COMPLEX A( 1 ) *0*01 
COMMON N * I GAM » DEL F »DELT 
COMMON /CDATA/ JCTR»DATA( 2oO> 

LOGICAL FFLG 
FC = DATA(JCTR) 

FCOFF = DAT A < JCTR+1 ) / 2. 

NR = DAT A( JCTR+2 ) 

FFLG = .TRUE. 

GO TO 20 
ENTRY SYNLP(A) 

FFLG = .TRUE. 

GO TO 10 
ENTRY SYNHP(A) 

FFLG = .FALSE. 

10 FO = 0. 

FCOFF = JATA(JCTR) 

NR = OATA(JCTR + 1) 

2 o TEST - EXPOS. / NR) 

CALL FRGFCN(A) 

X = EXP( .69314718 / NR) 

PSI = SORT ( X - 1.) 

FCPR = FCOFF / PSI 

00 50 I = l.N 

II ■ I - 1 - N/2 
F = II * DcLF 

FP = S I GN ( 1 . » F ) * (AbS(F) - FO) / FCPR 
IF (FFLG) GO TO 30 

IF (AbS(FP) .LT. l.E-16) GO TO 40 
FP = -l./FP 

3^ IF (A6S(FP) .GT. TEgT) GO TO 40 

01 = CMPLX ( 1 . » FP ) 

D * D 1 **NR ! 

A (I) * A (I) / D i 

GO TO 50 j 

40 A ( I ) = CMPLX ( 0 • *0. J 1 


50 CONTINUE 
RETURN 
END 





30. Subroutine TELPLT 




Called by: MAIN 

Calls: TIMFCN 

Commons : blank 

Entries : none 

Description: TELPLT is a routine which provides a remote teletype plot 

of the time waveform data. The routine provides the capability of plot 
ting selected portions of the array containing the data. This portion 
of the data array to be plotted is controlled by specification of the 
input parameters NST, NSP, NJUMP which specify the starting point in 
the data array, the stopping point, and the number of data array points 
skipped between plotted points. 




i 

i 


i 
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Program Listing: 


dUbROUTINE TELPLT ( A * NST *NbP»NJUMP) 


COMPLEX A ( 1 ) 

DIMENSION I A ( 50 ) 

COMMON N ♦ I GAM* DEL F *DELT *PD 
CALL TIMFCN(A) 
bMAX = REAL (A (NST ) ) 
bMlN = bMAX 

DO 1 I = NST ♦ NSP ♦ N JUMP 
d = REAL ( A( I ) ) 

IF(o.LT.bMIN) oMIN = d 
1 IF(o.GT.dMAX) bMAX = b 

IF( (dMAX-bMIN) •LT.l.E-3'-') GO TO 999 
DO 33 I = 1 » 50 
33 I A ( I ) = 1H 

low a R I T E ( 6 ♦ 4 ) 

4 FORMAT ( 5X * 1 H+ ) 

aR I T E ( 6 t 3 ) BMINtoMAX w .. lTt % 

3 FORMAT ( 12X*-AMPlITUl)L- MIN -»L9«4*-f MAX VOLT S- 

toR I T E ( 6 » 4 ) 


6 FORMATUxlzH 0 . 3X ♦ 2 H . 1 . 3X . 2H . 2 . dX . 2h . 3 » 3X ♦ 2H . 4 , 3X *2H . 5 » 

1 3X»2h.6»3X»2H*7»3X»2n«8»3X»2rt.9»2X»3nl.G> 


ihRITE(6*7) 

7 FORMAT ( lh,3X»lHI »lu( 5H I)) 

00 10 I = NST * NSP * NUOMP 
d =• REAL ( A( I ) ) 
u = ( b-oM I N ) / ( bMAX— bM I N ) 


M = INT(3*5u»+«5) 
IF(M.EQ.O) CC TO 34 


I A ( M ) = 1H* 

WRITE(6»35) I * I A 
I A ( M ) = 1H 

35 FORMAT (2X* I5.3H I * 5 ^ A1 ) 
oO TO 36 

34 wR I TE ( 6 * 37 ) I 
37 FORMAT (2X» I5t3H *) 

36 CONTINUE 

10 CONTINUE 

* R I T E ( 6 » 1 1 ) 

11 FORMAT (6X*1HN) 

GO TO 90U 

999 WRITE (6*9) 

9 FORMAT ( 1X*12HERK0K FINISH) 

90 o RETURN 


END 
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31. Subroutine TIMFCN 

Called by: AMDEMO, FMDEMO, FRQMUL, IDLMUL, LIM, PHDEMO, TELPLT , WRTF 

Calls: LFOLD, FFT 

Commons: blank, CDOM 

Entries : FRQFCN 

Description: TIMFCN checks the type of function stored in the data ar- 

ray (indicated by DOMFLG) and transforms when necessary. A call to 
TIMFCN will assure that a time function is in the data array, while a 
call to FRQFCN will assure that a frequency function is in the data 
array. 

Program Listing: 


SJoROuTINE TIMFCN (A) 
COMPLEX All) 

COMMON N » I oAV: 

COMMON /CDOM/ DOMFLG 
LOGICAL DOMFLG 
IFluCMFLG) GO TO 999 
call LFOLD(A.N) 

CALL FFT ( A , I GAM. 1) 

DGMFlG - .TRUE* 

GO TO 999 
ENTRY FROFCN(A) 

I F ( .NOT • DOMFLG) GO TO 999 
CALL FFT ( A, lGAMt-1 ) 

CALL LFuLD(A.N) 

UOMFLG = .FALSt* 

999 RtToKN 
c.ND 
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32 . Subroutine TTFP 


Called by : MAIN 

Calls : FRQFCN, SCALE 

Commons : blank 

Entries: none 

Description: TTFP is a routine which provides a remote teletype plot of 

the frequency function. The frequency spectrum is plotted between the 
limits specified as input parameters. These input parameters are FLO 
and FHI which specify the low and high frequency limits, in Hz, of the 
spectrum to be plotted. 

Program Listing: 


C 

C 


* * 


SUBROUTINE TTFPl A*FLO*FhI ) 

THIS SUdROUTINE PROVIDES A TELYTVPE 
SPECTRUM FROM FLO TO FHI AND PRINTS 
COMPLEX A ( 1 ) 

COMMON N * I GAM , DclF 
DIMENSION I A ( 50 ) »MM ( 6 ) 

CALL FRurCN(A) 


PLOI OF I HE 
THE CARRIER 


c * * * * 



NST = 

( N/2 ) 

+ 

NST = 

NST + 

1 

NoP - 

(N/2 ) 

+ 

NSP = 

NSP + 

1 


INT (FLO/DtLF S 1 6N 1 • 5 » FLO ) ) 
I NT ( FH I / DcLF + S I ON ( • 3 ♦ Fh I ) ) 


FREQUENCY 

frequency 


174 


TTFP (Continued) 


ubMAX = -1 • E 30 

I END = NSP - NST + 1 

DO 1 I = 1*1 END 

DECTMP = CAoS ( A ( NST + I - 1)) 

1F< JECTMP.LT. l.t-3o) DECTMP = 

b = 2 0. * ALGG1 o(DECTMP ) 

IF(b.GT.OriMAX ) JoMAX = b 

COiNT I Nut 


l.E-30 


WRITE (6*2) I END 
FORMAT ( /5X.8HNS IZE = 


I b/ ) 


CALL SCALE (DoMAX »MAXoCL ) 


THE ORDINATE WILL VARY FROM ( MAXbCL~5u ) 
TO MAXSCL Do. 

DO 33 I = 1.30 

> I A < 1 > = Ih 

DO 5 I = 1.6 

> Mi', ( I ) = MAXSCL - lu*(6-I) 

MAXF = AdS(FLO) 

IF(AdS(FHl l.GT.MAXF) M^XF = Ai 
NAMEF = 0 

IF (iMAXF.GT .1.L3) NAMEF = 3 
IF (MAXF.GT . 1 . t6 ) NAME : = 6 
I F ( M A X F • G T . 1 . E 9 ) NAME.- = 9 
I F ( NAMEF • EG * 0 ) WRITE I6«2«u) 

IF ( N AM t F • E G • 3 ) WRITt (6.2^3 ) 

I F ( NAMEF . EG » 6 ) WRITE (6.2^6) 

IF (NAMEF. Eo. 9) WRITE(6.2u9) 

J F0RMAT(2X,14HFRhjUENCY (HZ>»9 
3 FORMAT (2X*13HFRcOJENCY KHZ). 
b FORMAT (2X. 15HFREOOENCY (MHZ). 

9 FORMAT (2X, ISHFRlGJENCY (GHZ). 

WRITE (6 .7) (MM(I).l = 1.6) 

7 FORMAT ( / 7X» IA,4(6X.IA ) . 5>X * I 4 ) 

*R I TE ( 6 »d ) 

b FORVAT(9X.lHl.3(luH ♦ 1 


.*( 6-1 ) 


M&XF = AoS ( FH I ) 


(HZ ) »9X»6hUtCIot.Lb) 
KHZ) . dX . SHDEC I dElS ) 
(MHZ ) » dX . 8 HD EC 1 bELS ) 
(GHZ ) .ttX.bHDECIdtLo) 
. 6 ) 


4 


in 


TTFP (Continued) 


FFACT * 1. 

IF ( NAMtF.EU.3 ) FFACT - 
IF ( NAMEF .EC.6 ) FFaCT * 
I F ( NAMEF .EG » 9 ) FFACT s 
FLO * FLO * FFACT 
FHI = FHI * FFACT 
DELF1 * DELF * FFACT 
FLOPRT = DtLF 1* ( NST -1 
l»0 10 I = 1 ♦ I END 
UECTMP = CAbS ( A ( NST + 


l.t-3 

l.t-6 

l.E-9 


-N/2) 


- 1 >> 


I F (OECTMP.LT. 1 .E-3^ ) oECTMP = l.E-30 
e = 20. * ALOG1 o ( DECTKP ) 

M = 50 + 6 - MAXSCL 
J = I - 1 
XJ = J 

F REG * FLOPRT ♦ XJ * DELF1 
IF (M.LT .0 ) GO TO 5 o 
IF(M.EQ.^» GO TO 34 
DO 1515 II = 1 »M 
I A( I I ) = 1H- 

mR I TE ( 6 • 3 5 ) FREu.IA 
DO 1616 II = 1 »Fi 
I A ( 1 I ) = 1H 

FORMAT ( F8 • 3 »2n I»5oAl) 

GO TO 36 

*R l TE ( 6 » 37 ) FREU 
FORMAT ( IX. F7.3.2H -) 

GO TO 36 

inRI TE ( 6 *51 ) p RLQ 
FORMAT ( 1X»F7.3*2H I) 

FORMAT ( ) 

CONTINUE 

CONTINUE 

wRITE(6tll) 

FORMAT ( 6X*4HFREvj) 

RETJRN 

END 


l 


33. Subroutine WRTF 


Called by : MAIN 

Calls: TIMFCN, FRQFCN 

Commons: blank 

Entries : WRFF 

Description: WRTF generates a printed listing of the time function be- 

tween selected limits. Entry WRFF generates a printed listing of the 
frequency function between selected limits. 

Program Listing: 


■SuoRCuTlNE a'kTF ( A * iNi 1 * N2 ) 

COMPLEX All) 

COMMON N * I GAM ♦ uclF 
CALL TIMFCN(A) 

WRITE (6*704) 

DO 10 I = M.N2 
lo *RITt(6*7C0 ) I »A(I 1 

00 TO 999 

tNTRY RAFF!A»FLu»FhI ) 

CALL FROFCN(A) 

NCTR = N/2 + 1 

FLOT = FLO - DELF/1-. 

FhIT = Fnl + DcLF/l^« 

*RITt(6.7u2 ) 

DC 2u I = 1 » N 
F = (I - NCTR ) *0ELF 
IF (F .LT. FLOT) oO TO 2o 
IF (F • 0 T • FhIT) 00 TO 999 
THETA - 1 «c3^ 

T - C AbS ( A ( I ) ) 

IF (T «L T • l.E-30) 00 TO lb 

ThETA » 67.29378*ATAN2(AIMAG(A( I ) ) .KlALIAI I ) ) ) 

13 Do * 1*E3C 

IF (T .uT. v.) Do = 2 u.*AlOG1^( 1 ) 

■/»KlTE!6*7ol ) I »F»m( I ) »T *L'o»ThLTA 
2^ oONTINDE 
999 WRlTt(6.7C3) 

700 FORMATUh , l b.2! IPE1 -.3.2X )) 

701 FORMAT ( 1H » I 5 » 1 Ac. 1 2 • 4 * 3 ( 2 X * oPF d # 3 )» 2 X » F 7 • 2 * 2 X » F6 • 1 ) 
70 2 FORMAT!/- L I Nt-* bX * -FR£w-» 9X * -R t AL- *6X * -l MA(j-»6X • 

1 -MAG— » 6X *-Db — » bX ♦ -PhA Sc — / ) 

703 Pfo-TVAT I / ) 

7v/4 FORMAT!/- L I Nr.-» OX » -REAL - * 3X * - ! MAG-/ ) 

REToRN 

END 
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APPENDIX D 


Listing of SIGMA-5 Version of those FATCAT Routines 
which were Modified to Adapt them to the SIGMA- 5 - 


NOTE: Many of the Hollerith strings in format 

statements were delimited with quote marks. 

The printer used to make the following list- 
ing did not have the quote character; a minus 
iign appears where each quote should have been. 


pkrtrding page 


ELAN'K' Mo? 


FTT'nj) 
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» 


MAIN PROGRAM 


C 


COMPLEX A ( 1 024 ) 

COMMON N * I GAM »D£LF *l)ELT »PD .CARRFQ 
COMMON /CFREQ/ NFR »FR ( 6 ) 

COMMON /COOM/ OOMFLG 
COMMON /CDATA/ JCTR .DAT A < 2oO ) 
COMMON /CCIRXT/ NdLK ♦ I T YP ( 30 * 2 ) 
COMMON/CwORD/ wORD(12) 
COMMON/CFLGS/ PuFlG .ARFLo 
LOGICAL PDFLG, ARt-Lo 
JCTR = 1 

I TYPl 1*2) « JCTR 
NFR = 0 
NdLX = vj 
IdLlC = 0 


= .FALSE* 


ARFLG = .FALSE. 
WR I T E ( 6 » 7 006 ) 

1 JO 2 I * 1 ♦ 12 

2 WORD( I ) = Ah 
CALL FETCh(WORD» 


L • NbAJ) 


IF 

(N6AD 

• 0) 

oO T C 

i 1 






CALL ELFINJI 

WORD. 

LTYF 

1 > 





90* 

GO 

T 0 ( 10. 

20. 

3o . 

20. 

5 o . 

60. 

7u » 

80. 

no. 

120. 

liv^. 

1 A 0 . 

15o. 

160. 

1 7o » 

180. 

1 90 $ 

1 

210. 

220 . 

2 3o * 

2 AO t 

2 5 o . 

260 . 

2 7o . 

280 » 

c 90 t 

1 

310. 

320. 

330. 

3 AO . 

350. 

360. 

370. 

380. 

3 90 • 

► 

410. 

A20. 

<*3o ) ■ 

.ltyp 






IF 

( WORD ( A ) 

• EQ. 

1H ) 

GO TO 12 





N 1 

* WORD ( 3 ) 









N2 

* WORD ( A ) 









GO 

TO 15 










10 


12 WRITE(6»700A) 


RcAD ( 5 » 700o ) N1.N2 
15 CALL WRTFIA.N1.N2) 

GO TO 1 

2yj IF ( WORD ( A ) #tO. In > uO TO 22 
FRLO * WORD ( 3 ) 

FRH I * WORO(A) 

GO TO 25 

22 WRITE(6.70Q5) 

READ ( 5 » 70CO ) FRLO.FRH I 
25 IF (LTYP .£0. A) GO TO Au 
CALL WRFF ( A.FRLO.FRhI ) 

GO TO 1 

* * TTY TIME PLOT 
30 IF (WORD ( A ) .EG. lrt ) GO TO 32 
NST * WORD ( 3 ) 

NSP * WORD (A) 

NJOMP * W0RD(5) 

IF ( N JUMP .EQ. 1H > N JUMP * 1 
GO TO 35 


100 . 
200 * 
3 00 . 
400. 
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( 


MAIN (Continued) 


32 MR I TE ( 6 » 7007 ) 

READ ( 5 *7Q0U ) NST *NSP .NJUMP 
35 IF ( N JUMP .IT. 1) NJUMP = 1 
CALL TELPLT(A,NST.NSP.NJUMP) 

GO TO 1 

C * * TTY FREQUENCY PLOT 

4 u CALL TTFP(A*HRLO»ri\fiI) 
oO TO 1 

C * * CALCOMP TIME PLOT 

50 mRI T t ( 6. 71ul ) ((ITYP(I.J). J = 1.2). I = 1.5) 

7 1 u 1 FORMAT ( 1X.2( I 3.3X) ) 

6o 00 TO 650 

C * * PRINT PR I Me f-ACTORS 
7 0 CALL PRTFAC 

00 TO 1 

C * * cNO OF JOo 
oo 00 TO 999 

C * * DUTTERAORTn oANDPASS 
9*j CALL STRDTA(3»G»Q) 

NT YP = 3 
GO TO 6^o 

C * * ooT TERmGRTn lQMPaSS 
loo CALL ST R JT « ( 2 » 0 » 0 ) 

NTYP * A 
00 TO 600 

C * * cO T TER MQRT H HIOuPaSS 
110 CALL STRDT A ( 2 , 0 .0 ) 

NTYP = 5 
GO TO 6 v o 

C * * dOT TER aORTh dANDSTOP 
120 CALL STRUT A ( 3 . o .0 ) 

NTYP * 6 
00 TO 60o 

C * * CrtEoYSHeV uANDPASS 
1 3o CALL STRDTA(A.U.o) 

NTYP = 7 

00 TO 600 

C * •> CHEdYSHEV LOwPASS 
140 CALL STRDTA(3.0.0) 

NTYP = 8 
00 TO 600 

C * * ChEdYSHcV nIOhPAdS 
15o CALL STRDTA (3.0.0) 

NTYP = 9 
00 TO 6*0 

C * * CHtoYSMtV dANDSTOP 
160 CALL STRDTA ( 4 . 0 .0 ) 

NTYP » 10 

00 TO 600 


\ 
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MAIN (Continued) 


C * * SYNCHRONOUS oANDPASS FILTER 
170 CALL STRUT A ( 3 » 0*0 ) 

NTYP = 11 
GO TO 600 

C * * SYNCHRONOUS LOWPASS FILTER 
160 CALL STRUT A ( 2 * 0 *0 ) 

NTYP = 12 
GO TO 6uO 

C * * SYNCHRONOUS HlGHPASS FILTER 
19u CALL STRUT A { 2 * o *0 ) 

NTYP = 13 
GO T 0 600 

C * * SIGNAL GENcRATOR 
200 CALL STRDTA (6*2*1) 

CARRFO = WORD ( 3 ) 

NTYP = 1 
GO TO 600 

C * * FREQUENCY MULTIPLIER 
NTYP = A 
GO TO 600 
220 GO TO 430 

L * * IDEAL MULTIPLIER 
230 CALL STRDTA (2*1*2) 

NTYP = 16 
GO TO 600 
240 GO TO 430 

C * * FM DEMODULATOR 

250 CALL STRDTA( 1,0.0 ) 

NTYP * 19 
GO TO 600 

C * * PHASE DEMODULATOR 
260 CALL STRUTAC l.OtO ) 

NTYP * 20 
GO TO 600 

C * * AMPLIFIER 

270 CALL STRDTA ( 1 * 0 *0 ) 

NTYP » 2 
GO TO 600 

C * * LIMITER 

2dO CALL STRDTA ( 3*0*0) 

NTYP = 17 
GO TO 600 

290 NOUT * WORD (3) 

IFtNOUT .Lt. NOLA) GO TO 291 
WR I T t ( 6 * 70o 1 ) NdLA 
GO TO 1 

291 I F ( NOUT - IbLK.) 295*295*293 
293 I F ( PDFLG ) CALL PDCHa 

I TMP » I OLA ♦ 1 

DO 292 1 6LA * ITMP* NOUT 
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MAIN (Continued) 


IoTYP = »TYP( I 3L< » 1 ) 

JCTR = I T YP ( I BLK *2 ) 

CALL PRCCES ( I t?TYP * A ) 

292 CONTINUE 

IoLK * NOUT 

295 *R I TE ( 6 t 7uu2 ) IutR 
oO TO 1 
3 Uu CON T I NU t 
3 1 o CONTINUE 
320 00 TO 370 
330 GO TO 430 
C * * AM OEMOOULATOR 

340 CALL STROTm ( 1 ,G ,0 > 

NTYP = 22 
00 TO 600 

C * * ELAT SPECTRUM GtN.RATOK 
350 CALL STROTA < 3.0 »U ) 

CARKFQ = u • 

NTYP = 23 
GO TO 6^0 
360 CALL LSTCOM 

00 TO 1 
370 00 TO 65C 
360 CALL INPFOR 

GO TO 1 
390 GO TO 650 
400 GO TO 650 
410 GO TO 650 
420 oO TO 650 
430 WRITE(6*70o3) 

GO TO 1 

6v,G NdLiC = NdLK ♦ 1 

I TYP( NttLA.l ) = NTYP 

1 T YP ( NdL K. + 1*2) = OCTR 
GO TO 1 

650 RRITE(6»710v;) a/ORC ( 1 ) » IaORo ( 2 ) 

GO TO 1 
7000 FORMAT ( ) 

C 7000 F ORM AT ( 3 G • 0 ) 

70ul FORMAT!- * * ERROR * * LAKotST dLUCK. NO IS -»12»- * *-) 

7002 FORMAT!- PROCESSING oOMPlcTE ThRU OLUCk -»I2> 

7003 FORMAT!- * * uNjtHNtD sTATt-MLNT * *-) 

7004 FORMAT!- ENTtK lOM * HlsH INDICES-) 

7005 FORMAT!- ENTtK LOa» HIGH KSEu J c.NC I ts- ) 

7006 FORMAT!- START-) 

7007 FORMAT!- ENTtR NsTART* NSTOP* NJUMP-) 

7100 FORMAT!- COMMAND -»A4.A2»- IS NOT YE I CPE.R A I I ONAL- ) 

999 CONTINUE 
STOP 
END 
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* 


c 

c 


SUBROUTINE FETCH(WORD.LL*NBAD) 

INTEGER APCST. BLANK ♦CONN-. A, DECPT .EQUAL .PLUS .RP AREN *TEST 
INTEGER BUFFI .BUFE 2 .RUEF 3 .BCDFPT.TITLE.WORD.E 

DIMENSION T1TLEU2) 

COMMON/C fetch/ 3U C ”2 ( 6 ) ♦ BUFFI ( 80 ) 

DIMENSION WCRD(l) 

DAT A^APCST • BLANK *C0MMA .DECPT .EQUAL 1 H0^ » 1H+ * 

! / 1H- .1H . 1H . ♦ 1 H • *1H= *1H- .1H9 *1H0 , 1H+ / 

DATA LPARPN.RPAREN 
1 / 1H( . 1H) / 

DATA E / 1HF / 

L = 0 

NCOLS = 80 
NBAD = 1 

FETCH IS A FR^E-FIELO INPUT SUBPROGRAM WHICH RETURNS 
THE INpSt DATS IM LL CONSECUTIVE CELLS OF THE 
IrRaI SoRd! HOLLERITH IS TRUNCATED TO 6 CHARS. 

1 r?« ,N '(?.ioci.end = lO-T. I BUFFim.IM.BC 

WRITE <6. IOCS) ( BljFFK I ) .1 = 1*30 ) 

20 * - 0 

2 K * 0 
N = 3 

NCOMMA = 0 
DO 3 I = 1*6 

3 6UFF2 ( I ) = 3LANK 


) 


C 

C 


A 

* 

AO 


A 1 
A2 
A3 
AA 
A5 


IF 


( M-NCCLS ) AO » 1 » 1 oO 

exavinf each column, remove 

SEPARATORS. 

M = M+l 

TEST = BUFF 1 ( M ) 


blanks* and TEST FOR 


IF 
IF 
IF 
IF 
IF 

N = N+l 
BUFF2IN) 
IF ( K ) 


( 

( 

( 

( 

( 


TEST 

TEST 

TEST 

TEST 

TEST 


BLANK 

COMMA 

EQUAL 

lparen 

RPAREN 


Al.A.Al 
A2.6.A2 
A3. 6*43 
) AA .6 *A6 
) A5*A*A5 


= TEST 
5 » 5 ♦ A 
IF TYPE HAS NOT 
DETERMINE IF IT 
A A-8 PUNCH 
POINT NUM8EP) ♦ 
HOLLERITH WORD 
C WHEN 

1 WHEN 

2 WHEN 


3-EH SET ( K=0 ) ♦ TEST CHARACTER TO 
IS A DIGIT OR SIGN (NO DECISION). 
(TITLE CARD). A DECIMAL POINT (FLOATING 
OR NONE OF THFSE. IN WHICH CASE A 
IS ASSUMED. IF K IS SET. IT WILL BE 
AN INTEGER 

A FLOATING POINT NUMBER 
A HOLLERITH WORD. 
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'■ ^ n o n rs r. o o n 


* 


FETCH (Continued) 


5 NCOMMA = 1 

IF ( TEST - NZ ) 52*51*51 

51 IF ( TEST - NINE ) 4*4*52 

52 IF ( TEST - PLUS I 53*4*53 

53 IF ( TEST - MINUS ) 54*4*54 

54 IF ( TEST - APOST ) 55*30*55 

55 IF ( TEST - OECPT ) 56*57*56 

56 IF ( TEST - E ) 561*560*561 

560 IF ( N-l > 561*561*4 

561 K = 2 
GO TO 4 

57 K = 1 
GO TO 4 

* SELECT MODE OF CONVERSION* 6 AS EC UPON K • 

6 IF < K — 1 ) 7.7,9 

7 6UFP4 = GCDFPT ( BUFF2.N ) 

IF ( N ) 1-6*106*6 

B WORD(L+l> = ; 3UFF3 
GO TO 91 

9 ENCODE ( 4 « 1 002 ♦ WORD ( L + l ) ) ( B U - c 2 ( I ) *1=1*4) 
ENC0DE(2»l ri C? .WORD (L + 2H ( BUFF 2 E I ) * I =5 * 6 ) 

L = L + 1 


91 

* 


lu 

11 

# 

30 


IF 

IF 

LL 


C 

100 

C * 

c 

106 


IF NOT FINISHED WITH ThE CARD IMAGE* RF I N I T I AL I l r 
AND CONTINUE* IF THE NCOLS COLUMN CONTAINED A 
COMMA, PROCESS THE NEXT CARD. OTHERWISE* SET TH r 
NUMBER OF WORDS CONVERTED IN LL AND RETURN. 

( M-NCOLS ) 2*10*10 
( TEST - COMMA ) 11*1*11 
= L 
RETURN 

MOVE A TITL* CARD INTO THE TITLE APRAY. 

ENCODE! 48, IE 01 * TITLE ) ( E'JF-1 ( I ) * I = 1*48) 

00 TC A CARD IMAGE HAS BEEN PROCESSED. IF ThE LAST 
NON-BLANK SYMBOL WAS A COMMA ( NCOMMA sQ), READ 
THE NEXT CARD. OTHERW T 5E THERE IS INFORMATION 
IN *BUFF2* TO PE CONVERTED* AFTER WHICH, FETCH 
WILL RETURN TO THE CALLING PROGRAM. 

( NCCMMA ) 6*1*6 

p c t C H FOUND CONCOMITANT SEPARATORS OR A NUMBER WITH 
MORE THAN 15 DIGITS AND COULD NOT CONTINUE. 


IF 


WRITE ( 6 ♦ 2 C C r ) 
NB AD = 0 

RETURN 
107 STOP 

1001 F ORMA T ( 8 0 A 1 ) 

1002 F ORMAT ( 6 A1 ) 
1005 FORMAT ( 1X80A1 ) 
2000 FORMAT (25H0** 

END 


( BUFF 2 ( I ) ♦ I s 1 * 6 ) 


FETCH CANNOT DECODE 6A1.4H **//) 
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» 


SUBROUTINE INPFOR 
COMMON /CwORD/ wORD(12) 


INTEGER WORD 
LOGICAL FLG 
FLG * .TRUE. 


call 

ELF I NO ( WORD( 3 ) 

♦ L ) 



1, 

1, 



60 TO 

( 1. 

1* 

1 • 

1* 

1, 

1 » 

9 » 

A 

11. 

12. 

13. 

14, 

13, 

16 . 

17, 

18. 

19 » 

d 

21. 

1. 

23. 

1. 

25, 

26. 

27, 

28, 

1 * 

C 

1. 

1. 

1. 

34, 

35, 

1 , 

1. 

1, 

1 * 

D 

1. 

1. 5 

• OU ) » 

L 







I *RITE(6»700l> WORD ( 1 ) ♦ WORD ( 2 ) 
60 TO 999 

20 wR I TE ( 6 , 7020 ) 
wR I TE ( 6 » 7104 ) 

IF (FLG) GO TO 999 
35 WR I TE ( 6 , 7033 ) 

WRITE( 6*7111) 

IF (FLG) GO TO 999 
9 WRITE(6»7009) 

WRI TE(6»710l) 

IF (FLG) GO TO 999 
10 WR I TE ( 6 ♦ 70 1 0 ) 

WRITE(6.71o2 ) 

IF (FLG) GO TO 999 

II WRI TE ( 6 » 701 1 ) 
wR I TE ( 6 » 7 1 02 ) 

IF (FLG) GO TO 999 

12 WRITE(6»7012) 

WRITE(6.71ul> 

IF (FLG) GO TO 999 

13 WRITE(6,7013) 

WRITE(6»7101 ) 

WRITE<6»7103) 

IF (FLG) GO TO 999 

14 WRITE(6*7014> 

WRI TE ( 6*7102 ) 

WR I TE ( 6 • 7 lu3 ) 

IF (FLG) GO TO 999 

15 WRITE(6»70l5) 

WRITE(6*7102) 
wRI TE ( 6 *71o3 ) 

IF (FLG) GO TC 999 

16 WRITE(6.7016) 

WRITE(6»71C1) 

WRITE(6»7103) 

IF (FLG) GO TO 999 


10 . 

20 . 

1* 

1 . 
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INPFOR (Continued) 


17 WRITt(6.7Cl7) 
wRIT£(6.71ul ) 

IF (FLG) GO TO 999 

18 WRIT£<6»7018) 

'wRITE(6*71s>2) 

IF (FLG) GO TO 999 

19 *RITE(6.7019> 

WRI T E ( 6 * 71 02 ) 

IF (FLG) GO TO 999 

25 WRITE(6*7025) 

WRITE(6*71o7) 

IF (FLG) GO TO 999 
34 WR I TE ( 6 » 7u34 ) 

*RITE(6*7lo7) 

IF (FLG) GO TO 99> 

26 WRITE(6*7026) 
wRITE(6.71u7) 

IF (FlG) GO TO 99^ 

21 wRI TE ( 6*7321 ) 

WRITE(6»7105) 

IF (FLG) GO TO 999 
23 WRI TE(6.7023 ) 

WRI TE ( 6 » 7106 ) 

IF (FLG) GO TO 999 

27 WR I T E ( 6 » 7G2 7 ) 

IF (FLG) GO TO 999 

28 WRITE(6»7028) 

WRITE(6.711u> 

GO TO 999 

500 IF ( WORD ( 3 ) .EQ. 3HALL ) GO TO 510 
WRI TE ( 6*7112 ) WORD ( 3 * 

GO TO 999 
510 FLG = .FALSE* 


GO TO 2d 
7000 FORMAT () 

70u 1 FORMAT < IX *A4*A2 .- 
7wu9 FORMAT!/- owoNDP. 
70 1 d FORMAT (/- BWLQwP* 

7011 FORMAT ( /- own IP* 

7012 FORMAT!/- owoSTP* 

7013 FORMAT ( /- CHoNDP. 

7014 FORMAT (/- CHLOwP. 

7015 FORMAT!/- CHHIP* 


Ic NOT AN INPUT COMMAND-) 
F o * oW » NR- ) 

FC» NR- ) 
r C » NR-) 

F j » 6 W * NR- ) 

FO* OW. NR* EPSCB-) 

FC* NR* EPSDd- ) 

FC* NR* EPS DO-) 


INPFOR (Continued) 


7016 FORMAT </- 

7017 FORMAT (/- 

7018 FORMAT (/- 

7019 FORMAT </- 

7020 FORMAT (/- 

7021 FORMAT </- 
7023 FORMAT (/- 

7025 FORMAT </- 

7026 FORMAT (/- 

7027 FORMAT </- 
A 

7028 FORMAT (/- 

7034 FORMAT (/- 

7035 FORMAT ( /- 

7101 FORMATT- 
A 

7102 FORMATt- 

7103 FORMAT ( - 
7 1 u4 F ORMAT ( - 

A - 

8 
C 
0 
E 

7105 FORMAT ( - 

7106 FORMAT <- 
A 

7 lu7 FORMAT <- 

7110 FORMAT <- 
A 

B 

7111 FORMATt- 
A 

B 

7112 FORMAT ( / 
999 WR I TE ( 6 » 


CH8STP » FO ♦ Bw* NR * EPSDB-) 

SYNBP * FO* bW* NR- ) 

SYNLP* FC* NR- ) 

SYNHP » FC. NR- ) 

SIGGEN* FO* FMOL) * AM* PM* FM» A- ) 

FRQMUL-) 

1DLMUL* ALO* FLO-) 

FMDEMC. F O- ) 

PHDEMO* Fu-) 

AMP. GAIN-/- GAIN * VOLT AGt GAIN* Dd- 
(6 DB * FACTOR OF 2)-) 

LIM, CL* Cri* GL-) 

AMDEMO. FO- ) 
flatsp* AMP* OELF. N— ) 

FO = CENTER FREQ. HZ-/- BW = dANOWlDTH, HZ-/ 
NR = NUMBER OF SECTIONS-) 

FC = CORNER FREU* HZ-/- NR 


= NUMBER OF SECTIONS-) 


EPS08 = CrttbYShtV RIPPLE FACTOR* Ob- ) 


Fc = 
FMOO 
AM = 
PM = 
FM = 


CARRIER FREQ* HZ-/ 

= MODULATION FREu* HZ-/ 

PERCENTAGE AMPLITUDE MODULATION-/ 

PEAK. PHASE DEVIATION. RADIANS-/ 

PEAK. FREQUENCY DEVIATION* HZ-/ 

A = PEAK AMPLITUDE* VOLTS-) 

(NO PARAMETERS)-) 

*LO = PEAK AMPLITUDE OF LO SIGNAL* VOLTW 
FLO = FREQ OF LO* HZ- ) 

FO - CENTER FREQ. HZ- ) 

CL = LOW CLIPPING LEVEL* VOLTS-/ 

CH = HIGH CLIPPING LtVEL* VOLTS-/ 

GL = LIMITER GAIN. VOLTS/VOLTo- ) 

AMP * AMPLITUDE OF SPECTRAL LINES-/ 

DE.LF * FREG SEPARATION OF LlNtS* HZ-/ 

N = ARRAY SIZE-) 

- INPFOR CANNOT DECODE -»A6/) 

7000) 


RETURN 

END 
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SUBROUTINE ELF I NO ( NAME *L ) 

DIMENSION NAME ( 1 ) * MATCH(84> 

DATA ( MATCH ( I ) » I = 1 *84 ) 

A/4HPRI N.2HTT .4HPRIN.2HTF .4HTPL0.2HTT .4HTP10.2HTF. 
B 4HCPL0.2HTT.4HCPLC.2hTF . 4HPRIM »2H tF .4HENC0 *2 HF J . 
C 4HdWON»2HDP.4Hd*'LC.2h*P »4HowhI »2HP .4HbWdS »2HTP ♦ 
D 4rlCHdN»2rlOP»4ttCnL0»2HvgP ,4HCMHI .2hP »4nCHdS *2 hTP * 
E 4HSYNB » 2HP »4HSYNL»2HP .4HSYNH.2HP »4rtS IUG » 2 HEN » 
F 4HFRQM.2HUL .4h .2H »4H I DLM * 2HJL »4H »2 h » 

0 4HFMDE »2HMC * 4HPH JE * 2hMO * 4HAMP »2H »4HL I M »2H ♦ 

H 4HbL0C ♦ 2HK. »4HYES »2h .4HNO *2H »4HN »2H * 

1 4H »2h ,4 HAMl)E * 2HVO »4HF L AT ♦ 2HSP *4HL I ST *2HCO ♦ 

J 4HCIRC ,2HUI * 4H I NPU ♦ 2 hT F .4KDELE .2HTE»4HlNcF»2HRT » 



4HREPL .2 

HAC 

» 4HRE P E » 

2 HAT / 


NMAX = 84 





DO 11 I = 

1* 

NMAX. 2 



IF ( NAME ( 

1 ) 

- MAT Cl-! ( 

I ) ) 11.5.11 

5 

I F ( NAMc ( 

2 ) 

- MAT Cn ( 

I -► 1)) 11.21.11 

li 

CONT INUc 





I = NMAX 

+ 1 



21 

L = ( I ♦ 

1 » 

/ 2 



RETURN 

END 
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SUBROUTINE LSTCOM 
COMMON /CwORD/ WORD<12) 
INTEGER WORD 
WRITE(6*70ol) 

I = 1 



IF 

( WORD! 3 ) 

• 

EQ. 

3HALL) ( 

GO TO 2 

!0 


I = 

0 








IF 

( WORD ( 3 ) 

• 

EQ. 

4HS0UR) 

GO 

TO 

20 


IF 

( WORD ( 3 ) 

• 

EQ . 

ahfilt ) 

GO 

TO 

30 


IF 

( WCRD< 3 ) 

• 

EQ . 

ArtDEMO) 

GO 

TO 

40 


IF 

( WORD ( 3 ) 

• 

EQ. 

4HMISC ) 

GO 

TO 

50 


IF 

( WORD ( 3 ) 

• 

EQ . 

4Ht_0NT ) 

GO 

TO 

60 


IF 

( WORD ( 3 ) 

• 

EG. 

4HNSIZ ) 

GO 

TO 

70 


IF 

( WORD ( 3 ) 

• 

EQ. 

AHLIST ) 

GO 

TO 

80 


WRI 

TE(6.702 

0 ) 

WORD ( 3 ) 





GO 

TO 999 







20 

WRI 

TE< 6.7002 ) 







IF 

(I .EQ. 

0) 

GO 

TO 999 




30 

WRI 

;TE(6.700 

3) 







IF 

(I » E Q • 

0) 

GO 

TO 999 




AO 

WRI 

[ TE ( 6 ♦ 7004 ) 







IF 

(I .EQ. 

0) 

GO 

TO 999 




50 

WR] 

1 TE ( 6 » 7005 ) 






IF 

(I .EQ. 

0) 

GO 

TO 999 





60 WRlTt(6.70U6) 

IF (I .tQ. 0) GO TO 999 


70 WRITE (6.7007 ) 

GO TO 999 
80 WR I TE ( 6 * 7Go8 ) 

999 RETURN 

7001 FORMAT (/15X.-FATCAT COMMAND SUMMARY-) 

7002 FORMAT (/ 10* .-SOURCES-//- S I GGtN- . 5X *-S I GNAL GENERATOR-/ 
A - FLATSP-.5X.-FLAT SPECTRUN GENERATOR-) 

7003 FORMAT (/.l OX. -FILTERS-//- BUTT ERWCRTH — / 

A - 6W8NDP-.5X.-BAND PASS-/ 

b - BWLOWP-.5X .-LOW PASS-/ 

C - BWHIP-.6X .-HIGH PASS-/ 

D - 6WBSTP-.5X.-BANO STOP-/ 

E /- TCHEdYSHEFF — / 

F - CHdNDP-»5X *-oAND PASS-/ 

G - CHLOWP-.5X.-LOW PASS-/ 

H - CHH I P- »6X .-HIGH PASS-/ 

I - CH6STP-.5X .-dAND STOP-/ 

J /- SYNCHRONOUSLY TUNED— / 

K - SYNBP-.6X.-BAND PASS-/ 

L - SYNLP- »6X » -LOW PASS-/ 

M - SYNHF-.6X.-HIGH PASS-) 
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LSTCOM (Continued) 


7004 FORMAT ( /»1CX. -DEMODULATORS — // 

A - FMDEMO-»5X *-FM DEMODULATOR-/ 
p _ AMDEMO-.5X.-AM DEMODULATOR-/ 

C - PHDEM0-.5X. -PHASE DEMODULATOR-) 

7005 FORMAT ( / • 1 OX .- MISCELLANEOUS -//♦ 

A - FRQMUL-.5X .-FREQUENCY MULTIPLIER-/. 

C - IDLMUL-.5X. -IDEAL MULTIPLIER-/ 

0 _ AMP-. 3X .-AMPLIFIER-/ 
r _ LIM-.8X. -LIMITER-) 

7006 FORMAT ( / .1 OX .-CONTROL COMMANDS-// 

A - PRINTT-.5X. -PRINT TIME FUNCTION-/ 

1 - PRINTF-.5X. -PRINT FREQUENCY -UNCTION-/ 

r _ TPLOTT-.^X .-PRINTER PLOT CF TT'S FUNCTION-/ 

D - TPLOTF-.5X .-PRINTER ^LCT OF FREQUENCY FUNCTION-/ 

E - CPLOTT-.5X .-REMOTE PLOT OF TI"E FUNCTION-/ 

F - CPLOTF-.5X. -REMOTE PLOT OF FREQUENCY ^UNCTION-/ 

G - 5LOCK-.6X, -PROCESS SIGNAL TO 0‘JT D UT OF BLOCK SPECIFIED-/ 
ri - PR IMEF-5 X * -L I ST PRIME FACTORS CF ALL SOURCE FREQUENCIES-/ 
I - END CF JOB-. IX. -TERMINATES RUN-/ 

J - INPUT FORMATS * BLOC KNAMf LIST INPUTS FOR NAMED BLOCK-/ ) 

70^7 FORMAT (/ «4X .-SPECIAL COMMANDS TO SPECIFY ARRAY SIZE-// 

A - YFS-.8X. -LISTED ARRAY SIZE ACCEPTABLE-/ 

> _ N0-.9X.-LISTED ARRAY SIZE NOT ACCEPTABLE-/ 

C - N-. l'^X .-SET ARRAY SIZE TO SPECIFIER VALUE-/) 

7008 FORMAT!/.- LIST COMMANDS- COMMANDS ARE LISTED BY GIVING-* 

A - TWO ALPHANUMERIC-/- WORD SETS SEPARATED BY A CC V MA • -♦ 

B -THF FIRST IS -.1H-.-LIST COMMAND »-* 1H- /» - POSSIBLE-. 

C - SECOND WORDS AMD THE RESULTING OUTPUTS ARE — / 

D - ALL- » 1 2X *-L I ST ENTIRE COMMAND SET-/ 

E - SOURCES- *8X .-LIST COMMANDS FOR SDL'RCC BLOCKS-/ 

F - FILTERS-. 8X.-LIST COMMANDS FOR F I LTtr RS-/ 

G - DEMODULATORS-. 3X. -LIST COMMANDS FOR DEMODULATORS-/ 

H - MISC-.11X.-LIST OTHER BLOCK COMMANDS IAmp.LII .ETC)-/ 

I - CONTROL-. 8X, -LIST CONTROL COMMANDS-/ 

J - NSIZE-.1QX.-LIST COMMANDS CONTROLLING ARRAY SIZE-/ 

K - LIST COMMANDS-. 2X» -LIST THE ABOVE INFORMATION-/) 

7020 FORMAT!/- LSTCOM CANNOT DECODE -»A6/> 

END 


SUBROUTINE PDCHK. 

COMMON N.IC:A'NDELF,DELT*PD 

COMMON /CWCRD/ WORD! 12) 

COMMON /CFLGS/ PDFLG * ARFLG 

LOGICAL POFLC- »ARFLG 

NFLG = 0 

NSET * N 

CALL PER I CD 

IF( ARFLG) GO TO 

NSET = N 

v* R I T E ! 6 * 7 0 0 C ) PO.OELF 
1 v WRITE (6*70. 1 ) N.IGAM.OELT 
2-J CALL FETCH ( WORD. L.NsAD) 
CALL ELF I ND t WORD »LTYP ) 
IF(LTYP • FO • 30) GO TO 300 
I F ( L T YP • E 0 • 31) GO TO 3 0 
IF(LTYP .FQ. 32) GO TO 40 
IF(LTYP *EG. B) STOP 
WRITE (6 *7002 ) 

GO TO 20 

30 a R I T F ( 6 » 7 0 0 3 ) 

GO TO 20 
40 N = WORD ( 3 ) 

NFLG = 0 
CALL A 3 JN 

I r ( N .GE. NSET) GO TO 10 
NFLG = 1 
GO TO 900 

50 I F ( NSE T .OF. N ) GO TO 990 
WRITE (6*7005) N 
CALL PRTFAC 
WRITE ( 6*70 j 6 ) 

STOP 

900 I F ( NFLG .EQ. 0) GO TO 999 
WRITE (6 *70. 4) N.NSET 
WR 1 TE ( 6 * 7003 ) 


GO TO 20 
?9j N = NSET 
jqq POFLG = *FALSE • 

ARFLG * *TRUF. 

p rr 1 1p M 

30 r. FORMAT </ *- PFR ICO * -.1PE10.3.- FEC0N05 . -.-PELT A-F • 

A E 1 0 • 3 > 

r,-* FORMAT ( - N * - * 1 6 * - » I GAM = -.12*- » OELTA-T = - 
...t FOR t J pE1 . <3/ . , s THIS SATISFACTORY-. 

002 FORMAT!- INPUT MEANINGLESS » ENTER YES. NO. OR N .VALUE- 

003 FORMAT ( - ENTER N* VALUE-) , 

ijj4 FORMAT!- N * — *16*— UNACCEPTABLE ** N ’UST — 

ij«j4 FORMAT I N T0 MEET NYCUIST CPITFRION-) 

0C5 FORMAT!- SOURCE FREQUENCIES REQUIRE THE array SIZE - 

1 -TO 3E -* 19//- PRIME FACTORS APF-/) 

006 FORMAT!//- PUN IS BEING TERMINATED-) 

END 
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> 


1U 


2 \j 

b u 

7 0 U U 
9 99 


SUBROUTINE STRDTA(Kl»K2*X>3) 

COM VON/COAT A/ JCTR .OAT A ( 2us . ) 

COXMON/CFREO/ NFR *FK ( 6 ) 

COMMON /CFLGS/ PDFLG*ARFLG 
COMMON /CCIRsT/ NdLK. * ITYPl3v»2) 

COMMON /CwORD/ *CRD(12> 

LOGICAL PDFLU 

JCTK = ITYPtNdLN + 1*2) 

DO 10 I = l*Rl 

uaTm(JCTR + I -1) = 'aOKLHI + 2) 

JoTR = JCTR + <1 
I F ( K2 .LO, o) GO TO 999 
PDFLG = .TRUE. 

L)C 2 0 I = 1.K2 

TEMP = W0R0IK.3 + 1 + 1) 

IF (TEMP .LT. 1.) GO 10 2o 
NFR = NFR + 1 
I F ( NF R .GT. 6) GO TO bl> 

FR(NFR) = TEMP 
CONT INOL 
GO TO 9*9 
l«RITE(6*70w/u) 

STOP 

FORMAT!/- INPUT FREQUENCIES EXCttU Six-/) 

RETURN 

END 


j 
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LSTCOM (Continued) 


7004 FORMAT (/ » 1 0 X » -DEMOO JL A TORS — // 

A - FMDEMO- , 5x * -FM DEMODULATOR-/ 

' 5 ~ AMDEMO- » 5X » -AM DEMODULATOR-/ 

c __ ~ Ph, DEM0-*5X »-PHASE DEMODULATOR-) 

7w ° /O kmaT(/,1.x.- MISCELLANEOUS -//, 

A ~ FRQMUL- » 5X » -FREQUENCY MULTIPLIER-/. 

C - IDL v UL-*5x»-IDEAL MULTIPLIER-/ 

D - A m p - , g x »— AMPLIFIER-/ 
r - L !M-,8X .-LIMITER-) 

•ormat (/.UX. -CONTROL COMMANDS-// 

- PRINTT-t5X. -PRINT T IMF rUNCTT^/v./ 

- PRINTF-,5X, -PRINT FREQUENCY “UNCTION-/ 

‘ Tp L'DTT-,‘X,-PRlMTr3 PLOT OF Tivr FUNCTION-/ 

"»FX, -PRINTER PLOT OF frequency function-/ 

^PL'-.TT-, jX QT C PLOT OF TT-'E F! i\r T I QN- / 

I plot Of fffVOf.OCY PUNCTION-/ 

■3L0C <- *6X , -PRCC- p. c SIGNAL TO OUTDi'T OF Bi^Ck cprrrcr-' 

: ^Ir^rST Cf ilL ' S k’, 

u --*1X f- T nR*-. iNATf S p(j\-/ * 

"RMi??/ T / y 0R eJI5? t?L0C<>, ' AMF LIST X NPtJTS FOR N A v - h , L 0C<-/ ) 
^R-AT(/,4X,-SPlCiAL COMMANDS to SPrciFY A R R A v RI 7 “ ’// * X 

- YFS-.8X, -LIFTED APPAY SI7E ACC'PTA.BL^-/ Z " 

i ~ [\°~* 9y »-LISTcD array SIZE NOT ACCEPTA^., 

7ooa format.,.- lift :: -voo- $ . 


7 0 C 6 Fi 
A 
•1 
r 

n 

E 

F 

G 

H 

I 

J 

70,; 7 F { 

A 


70 ?; 


”TMF°F?3c? Af Tc V,r!, , rr " / “ ,V0R:) SETsGepARATED RY A c^ma . 
* T ^ F F IR ~ T ? -»1H-,-lIST COMMAND,-, 1H-/ PCSM&L"- 
S.:CONJ ’AC" OS AN 0 THE RFSULTIMG OUTPUTS ARE--/ *" * 

- ALL-, 12X ,-L I ST ENTIRE COMMAND SET-/ 

’ F?i U ?Ffc"’a X *’ L ! ST C0V ‘' ANDS F0K SOURCE nLOCXS-/ 
f lLTERS-,fiX,-LlST COMMANDS FOR FIIT-RS-/ 

DEMODULATORS-, 3X, -LIST COMMANDS p V demodulator f - / 

MI SC-, 1 IX ,-L 1ST OTHER BLOCK COM'‘A'DC (Amd.i ti r r f , . 
CONTROL-, ex, -LIST CONTROL Commands-/ L ‘ TCi-/ 

f*fi? fr T ,1?X, “ LIST COMMANDS CONTROLLING apqay <1?"-/ 

LIST COMMANDS-, ?X, -LIST THE A-^ V - I^n; W J4 I ' 
LSTCC-f CA-INOT dfcccf 1 ,WI 

END “ 


D - 
E - 
F - 
G — 
H - 
I - 
J - 
K - 
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